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ADAPTATION OF VENTURI FLUMES TO Re 
MEASUREMENTS IN CONDUITS. | 


By HAROLD K. ‘PALMER’, M. Am. Soc. C. E., AND FRED jak Bows, 


Assoc. M. Am. Soc. a E. 


sian » fe, Synopsis 


A weir can be Mdared as a control section through which’ water flows 


at critical depth. The sharp angle at the face and the fact that the flow is 
‘convex where it passes the critical section introduce energy losses between 
the point of measurement and the control section, the amount of these losses 
depending upon the setting of the weir in the channel. The ordinary. weir 
formula is an empirical equation that is accurate only as long as the funda- 
mental conditions upon which it was developed can be duplicated. These 
conditions can rarely be coraplied with in a confined channel, such as a sewer 


or an irrigation canal. 

This paper presents the theoretical hydraulic principles involved, and the 
results of special tests made, in the adaptation and construction of various 
Venturi flumes for measuring flow in conduits of uniform cross-section, where 


weirs have proved unsatisfactory. 
Since the uncertainties or variations in weir coefficients are due to eyes 


minable energy losses, it is reasonable to suppose that if these losses can be 


eliminated or reduced to a negligible amount by the use of some other device, 
the uncertainties in the rating curve will be eliminated. Such conditions are 
found in the so-called Venturi flume which, in this paper, includes any stream- 
d channel partly full, having 


lined device placed in an open channel, or a close 
critical depth with parallel 


a sufficient constriction to cause water to flow at 
filaments. 


hape of throat may be used, and the flow can be determined from a 


Any s 
] formulas. 


single depth me 
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The writers show how the rating curve may be EE and give Ree for 
ruse With circular conduits. : 
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INTRODUCTION 


’ 


Weirs have long been considered the standard devices for measuring run- 
ning water, but their use involves empirical formulas, and they must be built 
under restricted conditions, unattainable in many classes of conduits. The 
ponding effect of the water up stream from the weir causes deposits of sus- 

- pended matter which often alter the hydraulic conditions so that the empirical 
formula fails to give the correct flow. In the case of sewers, sludge deposited 
in this manner will decompose in time and become an added source of trouble. 

* Another objection to the use of weirs in a closed conduit is the relatively large 
loss of head. 

Many factors affect the proper installation of weirs as shown by Schoder 
Bh and Turner® in an able discussion of precise weir measurements. The United 
¥ States Bureau of Reclamation has found from experiments‘ that the ordinary 
me sharp-crested Cipolletti or rectangular type of weir must have certain very 
ek: definite requirements for a fair degree of accuracy; for instance, the distances 
: from the bottom and sides of the channel to the edges of the weir must not 
be less than twice the depth of water over the weir, and the channel above 
the weir for a distance of 20 or 30 ft must have a cross-sectional area at least 
six times that of the over-flowing sheet of water at the weir crest. When 
the depth of the water passing over a weir exceeds three-tenths of its 
width, the engineers of the Reclamation Bureau have found that the standard 
formulas indicate quantities that are too small. The error ranges from zero 
at three-tenths to 30% when the depth equals its width. These restrictions 
limit the capacity of the weir to about 40% of the capacity of a closed conduit. 

Occasionally, the true Venturi flume has been used to advantage because 
it gives a minimum of ponding and the frictional losses, being small, result 
in little loss of head. However, it has the disadvantage of requiring the 
measurement of the area of cross-section as well as the velocity of the water 
flowing in the throat. This has only been possible in the past by making 
two simultaneous depth measurements. 

Use of the Parshall meter,’ which is one form of a Venturi flume, obviates 
some of the difficulties previously mentioned, but it is not readily installed 
in conduits already constructed because of the required 3-in. drop in the invert 
grade at the throat. It could scarcely be adapted to an existing, standard 
sewer manhole, say, 4 ft in diameter. In the design of a new system of con- 


5 “Precise Weir Measurements,” by EH. W. Schoder, M. Am. 
Kenneth B. Turner, Hsq., Transactions, Am. Soc. C. B., Vol. 93 “(ibee) pes *S 94 mis vpse 


4 “Measurement of Irrigation Water’, U. S. Dept. 
mnivd Bitten sees: ept. of Interior, Bureau be Reclamation, 


5“ The Improved Neutne Flume”, by os Parsh 
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ey a me! Bee mH L 
ree eas grade sacrificed, Parshall meters could be (an 
ee neon 
or more of record ae. ious betaine esky inspection when a wee 
bor om equired, because the orifice between the channel and | 
. the wet-well often becomes clogged and the formation of sludge after a day 
Bor two in the wet-well tends to raise the float, thereby causing an error in. ote bah 
the stage recorder reading. Rae 

Ordinary Venturi tubes may also be used for liquid flow measurement, but 
larger structures are required for their installation. Such a tube must be 
entirely submerged and in sewers this is likely to prevent the proper handling 
of light floating material. Its efficiency is often impaired by clogging of the 
pressure recording pipes and, including the recorders, the initial installation = 
_ is more costly than the weir and stage recorder units. a ; 
~~ Tn 1920, Hinds® suggested placing a sufficient constriction in a channel to. 
-. eause the water to flow at critical depth, on the assumption that in this case 

the energy head is fixed, whereas. the depth is uncertain. He presented experi- 
ments made by the U. S. Bureau of Reclamation,’ on a flume similar in many 
respects to one type considered by the writers and showed an error of less — 
" than 5% in measurement of flow, assuming no transition losses. 
. The following analysis of the adaptation of the Venturi flume to flow 
4 measurements in conduits of regular cross-section is applicable not only to 
= clear water, but also to any liquid carrying suspended matter. Among the 
E latter are included silt-laden irrigation water, storm water, and sewage. 

The investigations herein reported were made in connection with a sanitary 
sewerage system and as sewage presents most of the difficulties encountered 
in the other classes of water, the word, “sewage,” will be used subsequently 
to denote water that carries settleable solids. 

In passing, it may be stated that an accurate measurement of sewage 
« flows in various parts of a sewerage system is quite important in order to be 
able to predict with reasonable accuracy the need of future works, relief 
sewers, additional pump capacity where pumping is required, ete. In 
é joint sewerage works an accurate measurement of sewage is necessary to 
__ apportion costs to respective participants. Furthermore, as a check on pump- 

ing efficiency and power requirements, knowledge of the actual sewage flow 
is necessary. An accurate gauge of the flow is required for efficient and eco- 
nomical operation, where pre-chlorination® is practiced since the chlorine feed 
depends on both the strength and the quantity of the sewage. One difficulty 
encountered in an accurate measurement of sewage flow in new sewerage 
works is the fact that a sewer is usually constructed of a greater capacity - 
than is required at the time it is first placed in operation. The flow in the 
line is then relatively shallow. Furthermore, in a new system of any- large 


e “Venturi Flume Data Throws Light Upon Control Weir’, by Julian Hinds, M. Am. 


} Soe. C. B., Engineering News-Record (1920), p. 12238. 
7“~The Improved Venturi Flume”. by R. L. Parshall, Assoc. M. Am. Soc. C. B., dis- 
ae Soe. C. B., Vol. 85 (1926), 


cussion by Julian Hinds, M. Am. Soc. C. E., Transactions, Am. 


p. 864. 
8 “Control of Sewage Condition by Chlorination”, by F. D. Bowlus, and A. P. Banta, 
Assoe. Members, Am. Soc. C. B., Water Works and Sewerage, November, 1932, Vol. 79, 
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fat! ‘fees oul Brae 
ficulties see eed ae Ba see weirs is othe’ tole’ denibea at Ef 
> measurement ‘of sewage indicate that the most satisfactory meter would | i} 
corporate the use of a’stream-lined flume, with minimum head loss, pre- 
ing deposits of sludge above it; without sharp édwes to catch rags ere 
er floating refuse; and with sufficient capacity to measure the entire designed 
7 flow and yet be sensitive enough to record low flows in the first years of its 
use. The flume must be adaptable to the use of a simple float mechanism for 
_ -water-stage recording, and must be designed so as to be readily installed in 
an existing conduit with minimum interference to the sewage flow. This 
paper deals with a measuring flume designed so that the flow may pass without 
objectionable ponding up stream; which offers no obstruction to floating or 
suspended solids; and which affords an opportunity to measure the flow from 
a continuous record of the water-surface elevation at a single point in the 
channel up stream from the flume. 
ay The formulas developed for the flume are not empirical but are based 
a rf entirely on theory and can be applied to conduits of any regular size or shape. 
ee _ Trregular channels would require regulation for a short distance up stream 
Soe by lining’ of some kind. Uniformity of shape rather than factors affecting 
ii friction are important. The method has been developed especially for sewers, 
wie but is equally applicable to other channels, irrigation canals, storm drains, ete. 


Notation 


Py ‘. The symbols introduced in this paper are defined as follows: 


a bottom width. 

a subscript denoting “a critical depth.” 
“function of.” A 
acceleration due to gravity. ; 
Kutter’s coefficient of roughness. 
a subscript denoting “velocity.” i 
area. 


surface width; width of water surface. 
Chezy’s coefficient. 

diameter. 

head; Hy, = velocity head. 


a constant in Kutter’s formula — A QAM R. 

rate of flow, or discharge. 

hydraulie radius. 

hydraulic slope; slope of a conduit, expressed as a percentage. : 

average velocity in a section. | 

“difference in;? AH, — difference in velocity heads in throat 
and conduit. 

energy head; «. = energy head at. critical depth.' 

“function of.” 
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conduit is the height of the water surface above 
the velocity head. In passing from one section of a conduit to 
nal a: id -lossés are measured by the decrease in the total 
changes in depths and velocities at the two places. 


urement of the total energy head is the logical basis for the de 


ment of any theory in the design of a Venturi flume of the general _ 


illustrated in Fig. uN : 


eee eee | 
Stage 
Recorder 


Critical Depth 
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(a) SECTION AT ay: 
CENTER OF THROAT (6). LONGITUDINAL SECTION OF METER IN 24-IN. PIPE . 


(c) PLAN OF FLOAT SUPPORT 


Under certain conditions a calculation of the flow may be made by simply 


measuring the depth of the water flowing in a conduit of regular cross-section — 


at any convenient point immediately above the throat of the measuring flume. 
This is only possible when the throat is designed to cause the water to flow 
‘through it at critical depth. For any given quautity of water, and size and 
shape of conduit, Fig. 2 shows that the energy head has a minimum value 
when the water is flowing at eritical depth in the throat. For a given channel 
this critical depth is a function of the quantity of water, and attempts have 
been made by Woodburn” to use a broad-crested weir as a measuring device 
recording the critical depth. However, in this case, it was found impossible 
to locate the section of critical depth accurately as it changed position with 
different quantities of flow. Near the section of critical depth in the throat 
of the Venturi flume, the depth is uncertain, but the energy head is fixed; 
therefore, if the throat is designed so as to cause the water to flow at critical 
depth the depth in_ the conduit of the regular cross-section above will 
be the average energy head less the velocity head at that point, except for 


frictional losses, which will be discussed subsequently. 


9 “Mests of Broad-Crested Weirs’’, by James ve Woodburn, Assoc. M. Am. Soc, C.7H.; 


“Transactions, Am. Soc. C. E., Vol. 96 (1982), p. 387. 
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Energy Head, € in feet Rate of Flow, Q, in cu ft per sec 

Fie, 2.—EHnerGy Hap CURVES IN A VENTURI FLUME. Fic. 3.—DERIVATION OF 

d . RATING CURVE, VEN- 
TURI FLUME. 


For any given quantity of water, according to Bernoulli’s theorem, 


A G0 Oe ee ASSES ine 8 ee (1) 


eo 2 | 
ie -=4+($) 2 oe he ae ae (2) 
= Ay 2g 
a. in which Q = rate of discharge; and A = area of cross-section. For a 
( rectangular cross-section: ; 
‘ ASS bd: 3. . dupa Jesiesda Shee (8a) 
i : 
; . Q 
: bd ~ 
and, } 
2 
-=a+(2 dciéhen Sateuston ar ee | 
bd J 29 } 
in which b = the bottom width of a channel; and e and d are the vari- 


ables. In Fig. 2, Curves F and T represent Equation (4) for the conduit and 
throat, respectively, Q being the same in both cases. Neglecting frictional 
losses there would be no loss of energy in or through the Venturi flume, and 
if Q@ cu ft per sec flowed at a depth, A, in the conduit, it would flow at a 
depth, A’, in the throat. Theoretically, Q could be determined by measuring 
these two depths, A and A’. The necessity in the past for the two measure- 


~ 


% showing the RES between Q and d could be drawn, but in practice 2 
_ this requires an unnecessary amount of labor. 4 oe 
oe ‘For a given size and shape of throat, the value of «> is a definite Sanoned Chas 
_ of Q, or to write the inverse case: : i 


a aC or cS See gepiehp sd) 


When e and Q are known, Hy, in nie conduit can be found by trial, and sub- 
tracting this H, from e gives d (Fig. 3). 
 - For rectangular throats, Equation (5) takes a simple form which will be 
evaluated from the general formula, but for other throat shapes, this becomes _ 
a transcendental function too complicated to be considered, and it is easier — 
to compute a few points on the curve and prepare the rating curve graphically. . 
Tt has been found from practice in computing these points that the best pro- 
~ cedure is to assume various critical depths in the throat and compute Q and e 


2 
2 by the following formulas: 
2 
4 
a 


At critical depth, se — 0; differentiating Equation (2): 


ey. . de Q? 1 dA ( 
: Ah Sh eee ESAT tee Ee eee pea ee 8 6) 
pe dd g es A$ Es dd : s ‘ 
a F . . . : ty dA ¥ 
Z Equating this formula to zero, and substituting B for its equivalent®, ae %, 
Equation (6) becomes: a 
/ ; apes nt | 4 
, Pewed Sagtent Canes gteie ae Bas (7) : 
and, therefore, f 
corns ri 
A 
ny Fag Y fee 1 Sa On he RE SOR eR a C8) 
o B 


Equation (8) is the accepted formula for the quantity flowing at a critical 
depth in a channel of any shape. Since Q@ = AV, by substituting the value 


of Q in Equation (8): 


iE AMET PET ey tcc ee (9) 
29 2B 
and Bernoulli’s equation becomes, 
: 2 TE ae ale enn eee (10) 
| a 


1 See SS 
fe A of Open Channels”, by B. A. Bakhmeteff, M. Am. Soc. C. E., vp. 31, 


Equation (19). 


va 


2 ‘ a Substituting the value of de obtained from Equation (13) in Equa- 
pe Te Re Rater of sengtab yogis tM SCI Bee De ee 
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which is the form taken by Equation (5) for a rectangular throat. For all 
_ other shapes it is necessary to use Equations (8) and (10). Devaste 
- Having drawn the quantity-energy head curve for a given throat (see — 
. Fig. 3), it is necessary to compute and draw the quantity-velocity head curve 
for the conduit section. If the conduit is circular this may be done by 
means of Fig. 4, which gives for various depths the velocity and velocity head 
for 1 cu ft per sec on a logarithmic scale, and, for other quantities, by merely’ 
adding log Q graphically. The intersection of the velocity abscissa with the 
“guide” line falls on the ordinate of the velocity head. The procedure is best 
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y ~ sects the “guide line” at Hy = 0.225. In the second case let Q = 
per sec, with water flowing 0.27 ft deep in a 36-in. pipe. Set the Fa 
_ tthe distance 1.0 to the left to 0.6 and measure this distance to the left from 
-— imtersection of the 36-in. pipe with the d = 0.27 horizontal line, finding» 


4 & ‘note the assumed depth. Enter Fig. 4 with these values of Q and d and read : 
the velocity head. If this differs materially from the assumed value, a second 


' quantity-depth curve is the required rating curve. 


a= 035 os giving @ =3.8 ‘ft per sec. The ar de line for v= Av 


“So “= 1.86 mae per sec. The vertical line for v = 1.86 intersects the guide Tine. 
To find the velocity head for any, value of Q, first assume an pe ttees ie i 
value ‘of the velocity head, subtract it from the quantity-energy head, and es 


approximation should be made. To draw the velocity-head curve it is obvi- 
ously better to begin with a small value of Q because the error in the assumed 
value of the velocity head is correspondingly small and one or two approxi- Win 
mations are usually sufficient. It will be found that only a few points are 
required to draw in the quantity-velocity head curve. Ordinates on the 
quantity-depth curve are merely the differences between the ordinates of 
the quantity-energy head curve and the quantity-velocity head curve. This 


As an illustration of the method of computing a rating curve, consider a 
Venturi flume with trapezodial throat cross-section as in Fig. 1(a), 12 in. 
wide at the base with side slopes of 2 on 1 and a flat slab base 3 in. above 
the invert of a 24-in. pipe. All measurements are in feet and quantities in — 
cubic feet per second. As previously recommended, various critical depths 4 
will be assumed and Q and « computed. In this case the formulas for the 
‘eross-sectional area and water-surface width in the throat are: 


A= ( + “) td oe. Racine alrite sean (15) q 


aad Pe ic tates) wai eeecnen ne 


and, 


The computations for energy heads are arranged in Table 1. In comput- 
ing the velocity heads (see Table 2), it must be remembered that the crest 
of the throat is 0.25 ft-above the invert; hence, the depth in the conduit must 
be increased by that amount. The velocity head in the conduit can be 
obtained only by trials preferably starting with the small values OLE Onan 
Table 2 when Q = 2 cu ft per sec, the depth in the conduit will be somewhat 
less than e + *0.25 = 908. For 2 cu ft per sec, at a depth of 0.9 ft in a 24-in, . 
pipe, Fig. 4 gives Hy > 0.03. Assuming Hy = 0.03, the depth becomes 0.88 
with a resultant value of H, = 0. 035. This is not a sufficient change to effect 
the result so 0.035 can be assumed as the final value. Since the crest of the 


nay be 0 obtained by. ecddina thd o curve drawn throng ou 


as 


_ TABLE 1—Quaxtiny-Exenoy 1 Han ) Ovnve 


“By weit and pinniniates Knee Fig. 1). 
Table 1 is extended to a value, Q = 11.0. The energy head was so regu- 
lar that it was extended by eye to Q = 12, and Table 2 is computed to this 


limit although it is subject to some suspicion for Q > i. 0. Thasexs curves are 
all shown in Fig. 3. 


TABLE 2.—Quantitry-VeLocity Heap Curve. 


Q = 2 4 6 8 10 12 
0.658 0.978 1.225 1.438 1.630 1.812 
0.908 1.23 1.48 1.69 1.88 2.06 
0.03 0.07 0.09 0.14 0.19 0.25 - 
0.88 1.16 1.39 1.55 1.69 1.81 
0.035 0.064 0.100 0.145 0.198 0.256 
0.623 0.914 1.125 1.293 1.432 1.556 


An exponential formula may be written for some small sections of the 
final rating curve, but it is not advisable. The inaccuracies of such a 
formula can be appreciated by equating the value of « in Equation (2) for 
both the throat and conduit sections; thus, 


(aria Re yet 
a+ (2) + =a+(2) Societe sess ees (17) 

c “9 A “9g 

which reduces to, 
Q= 29 ar ee .. (18) 
A—4s 

Writing A = f(d) and Ac = ¢$ (de): 4 

Q=1J29 CH= WIMOOG Pee .. (19) 
f (d) — o (d.) 


Since f(d) is not a simple function in a circular pipe, it is evident that 
the use of an exponential formula is generally inadvisable because of the 
many different combinations depending on the relative size of the pipe and 
throat, and the height of the crest of the throat above the invert. 
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| the size and setting. This change in the constant can be interpreted as 
variation in the energy loss, and it follows that if such variations can 
_ limited by the use of other devices, the change in so-called constants w 
likewise be restricted. In the case of the weir such control is difficult an 
am at times impossible; with a Venturi flume it may be accomplished easily. 
_ In the adaptation of the Venturi flume, to the writers’ ideas, considerable _ 
4 attention was given to the design of the transitions to minimize loss ot) 
energy. In the hydraulic design of flume transitions Hinds" found, an 
_ making twenty-nine tests on ten flume inlets, that only three had losses 
~~ amounting to more than 0.1A Hy. Ina number of cases no measurable loss - 
-oceurred and the average was about 0.04 A Hy. The 4-ft diameter of a stand- _ 
ard sewer manhole limits the length of the transition that can be built” 
readily in an existing sewer. In the sewerage system of the Los Angeles © 
County Sanitation Districts, a ratio of three longitudinally to one trans- : 
versely has been adopted as a compromise. To determine the energy losses 
through Venturi flumes constructed by the Districts, a differential meter 
acting on the principle of a Pitot meter was devised. This is shown in Fig. 5. 
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. Fic. 5.__DIFFERENTIAL ENERGY Fic. 6.—ENERGY CURVES, ConrroL SEecTioN WHIR. 


METER. 


Two bent tubes are set at two points in the stream with orifices at the same 


elevation and pointing up stream, one in the throat and one above the inlet 


transition. These two tubes are connected by means of air tubes to a mano- 


meter filled with colored water; the displacement of the water column in the 
manometer after all water js blown out of the bent tubes measures the differ- 
ence in energy head at the two points. To allow for any compression of the 


Flume and Siphon Transitions’, by Julian Hinds, M. Am. 


; “ i ign of 
Se ae orhada Arn: B., Vol. 92 (1928), p. 1428. 


Soe, C. B., Transactions, Am. Soc. C 


¢ ha ghee not sensitive to such small changes this energy loss is practically 
negligible. Should an appreciable loss be found in a given installation for 
% high velocities it can be added to the quantity-energy head curve. 
In sewage measurement, the mean daily flow is usually required. This 


ins we the lohstotintiod an been cleared away and a little in a ! 
forced into the tubes. Any type of meter in which water is drawn into the — 
oral i be clogged by sewage. This simple, head-differential pean has — 


Biederat tests with the dhitcrintaltutees distes (Fig. 5) infeeeea over 
entire cross- “section have shown that no energy ir greater than 0.005 ft 


may be calculated with sufficient accuracy by using a planimeter on a water- 


‘stage record to get the mean depth. If it is suspected that there are energy - 
‘losses through the Venturi flume at higher flows the differential energy meter 


can be used at various stages, including the maximum, to correct the curve; 
and if the flow is increasing from year to year, the corrections may be deter- 


amined beforehand, as the maximum flow of to-day will be the mean flow some 


-months or some years hence. 
Another possible source of error seems to be that in-an open conduit the 


“surface velocity is about 25% greater than the mean, which should cause 


the velocity head at the surface (where the float is installed for measure- 


ments), to vary as the ratio of 1.00° to 1.25*, or 56% greater than the mean. 


If this condition existed, 1.56 H, should be subtracted from the quantity- 
energy head instead of 1.0 Hy, in order to obtain the proper depth in making 
a rating curve. This condition was tested at one meter formed by placing a 
flat slab, 6 in. thick (at the invert), on the bottom of a 54in. sewer pipe 


without side contractions. In the regular pipe channel, at the upper end of 


the transition section and at mid-stream, a Pitot tube was turned up stream 
and then at right angles to the current. The difference between these two 
measurements was the true velocity head, which, added to the measured depth, 
gave the total-energy head. Applying this energy head to the previously com- 
puted quantity-energy head curve gave the quantity. Readings taken at 
various depths between the center and sides of the stream showed substantially 


TABLE 3.—Retation Berween THE MEASURED AND THE THEORETICAL VELOCITY 
Heaps at Mip-Stream. 


‘Quantity, Q, in cubic Calculated velocity, V, Theoretical head, He, Measured head, H» 
feet per second in feet per second ip feet infeet | 
8.2 2.18 7 0.074 
10.25 2.50 0.097 0.006 


17.5 3.15 0.153 0.142 


« 


icate ‘that within the nae ‘of the Serene the veloci a 
pee it may be assumed to be the mean velocity head. 
% begs investigation os the thane tants correction for Seon flo 


be eet For measurements in clear water, using a float-well for more aN 
"accurate measurements, it should be determined experimentally. 
_ }..,,. When the velocity above the meter is less than 2 ft per sec the saa oes 
) _surface is usually smooth, but above this velocity waves begin to form, — 
probably caused by the difference in velocity heads in the region adjacent * 
to the sides of the conduit. The waves are too small to affect the ball-float, 
but would affect very precise measurements made with a hook-gauge. Where 
_they are usable float wet-wells would obviate this trouble, but in sewers the 3 aL 
floats are affected by the accumulation of decomposing sludge. ae 
One question that may be raised is the effect on the energy head of Kees a 
jc ale velocities at various points in the cross- -section of the throat. It was 
found that the energy head was the same at all points in a transverse plane. . 
As another check, small floats were dropped simultaneously into the upper 
- end of the throat in pairs, one being near the side and the other nearer) 
J the center. Every time, both floats reached the lower end together. As 
a “another test a string was dropped into the throat transversely to the flow 
line and reached the lower end straight.. The explanation of this phenomenon 
is that at, and near, the critical velocity in a short flume the side friction is 


» 

s _ not an appreciable factor. : 
; Efforts were made to find a coefficient to be applied to Hy» in the conduit. i 
The flow was determined by measuring the energy head in the throat with a & 
Pitot tube. Using this value of Q, the measured depth was subtracted from a 


the energy head to give C Hv. Within the accuracy of measurements, flows 
“of 10.3 cu ft per sec, at a velocity of 2.4 ft per sec, and 17.5 cu ft per sec, 
at a velocity of 3 ft per sec, in the conduit both gave C = 1. These were ) 

the largest flows and velocities available. , 


5 : EXPERIMENTAL VERIFICATION 


In the experiments made by the U. S. Bureau of Revlantation,” a wooden 
throat with a bottom width of 2 ft and side slopes of 1 on 15 was installed 
y in an earthen canal with a bottom width of 8 ft and similar side slopes. 
5 The bottom of the throat was 8 in. above the base of the canal and the down- 
| oat sloped downward to a point 8 in. below the 
similar to the Parshall meter. The gauge-point 
here the bottom width was approximately 


y stream | apron end of the thr 
é elevation of the canal base, 


was placed in the upper transition W 


ace ection for Hydraulic Flow”, by M. P. O’Brien, Assoc. M. Am. 
Ce bunt Bee daa ae Jun. Am. Soc. C. E., Engineering News- Record, Vol. 113, 


Soc 
1934, p. 214. 
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riters have itinas energy, curves for the flume (Fig. 6) at t 
or throat section, the gauge-point, and in the main canal for +50 and 
ft per sec, respectively. ; 
; _ ‘Had the gauge-point been moved back into the main canal, the measure do 
depth would have been greater by about 0.1 ft in the case when Q = 50 cu ft 
per sec, but the main argument-in favor of moving up stream is that the 
drop in the water surface indicates convexity, especially as it increases rapidly 
with the tapering of the channel, and it is in the convex water surfaces that — 
se the relation of the water surface to energy head is not correctly given by the 
honey Bernoulli theorem. Had the gauging-point been placed in the regular channel — 
above as suggested, and the sides and bottom lined for a few feet to maintain 
a uniform cross-section, the stream filaments would have been parallel and ~ 
the calculated flow would have followed more closely the volumetric measure- — 
ments. Friction due to the longer distance between the point of measurement — 
and the throat could have been measured by means of the differential energy- 
head meter (Fig. 5) and the energy head increased by this amount. — 

With this change in the point of measurement the recorded tests would 
have shown closer agreement with the calculated rating curve. Since the 
quantity of water flowing through this Venturi flume tested by the Bureau 
Sg of Reclamation was carefully measured volumetrically the data are considered 
most reliable. From the known dimensions a quantity-depth or rating curve = 
was readily constructed by the writers using the method previously outlined. 
As stated by Hinds the difference between the quantities taken from this 
curve for given depths and the actual measurements in no case exceeded 5 per 
cent. Since the writers have pointed out that greater accuracy could have been ; 
attained by the selection of a better place for measuring the depths, and by 
measurement of frictional losses, this test is sufficient proof that all forms of 
the Venturi flume will give accurate measurements when it is possible to 
measure the depth in the conduit of a regular cross-section immediately above 
the throat and where the throat is properly designed to insure that the water 
reaches critical depth. 

In the absence of an experimental laboratory the writers have had to 
make use of active sanitary sewers in the development and adaptation of the 
Venturi flume designs. Standard weirs were set in adjacent manholes, but — 
comparative results were poor except where the flow was small in comparison 
with the capacity of the sewer, since the weir coefficients were uncertain. 
Therefore, these measurements were not considered sufficiently accurate to 
verify those obtained from the Venturi flumes. 

The experiments made by the U. S. Bureau of Reclamation, previously 
referred to, were considered a positive check on the accuracy of the method 
outlined herein. Furthermore, a field test, comparing measurements of a 
temporarily ‘installed Venturi flume with a permanent Parshall meter (made 
through the co-operation of the Cities of Los Angeles and Beverly Hills), 
supplied added proof of the accuracy of the adapted Venturi flume. 
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eter and an indicating and recording 

agh the flume and the rate of flow, in cubic feet per second. Water con- ; 
ously ee into the float-well prevented the accumulation of sludge — aS 
ich otherwise would affect the float results. The wooden throat of the 
mporary Venturi flume was installed in an ordinary standard manhole 
bout 1000 ft above the permanent gauging station, where the grade of the © 
3 1-in. sewer was 0.70 per cent. The bottom of the throat with a width of 10 in. 
owas placed 2 in. above the invert. Although the boards forming the side 
4 slopes were intended to be set at a slope of 2 on 1, the irregularities of the — Bs 
pipe caused them actually to be placed on slopes of 2 on 1 to 12 on 1. 
= Tin transitions for stream-lining the upper end were attached to the 
_ boards, but could not be set in exact position due to the roughness of the chan- | 
_ nel through the manholes. Some caulking had to be done with oakum. The — 
s entire apparatus was quickly installed in the sewer during the low early 
morning flow and readily withdrawn a few hours later when the flow was 
at its peak of about 6 cu ft per sec. The elevation of the bottom of the 
3 throat was referred to a straight-edge placed across the manhole rim in 
- the street pavement and measurements to the water surface above the flume free 
were made with a steel tape. 
a _ Readings were taken simultaneously at the Venturi flume and the Par- 
- shall meter, a correction in time being made between the two stations as 
; determined by passing floats occasionally from one to the other. A correc- 
_ tion was also required for the small flow entering the sewer from a side branch a 
- at a point between the two testing stations. Occasional depth measurements ee 
‘over a V-notch gave this correction with sufficient accuracy. 

Table 4 gives the results of tests when the flow was within the capacity of ' 
the Venturi flume. In Table 4(a) the average discrepancy between the two 
methods of measurement was 4.2%, with a probable error of + 0.4 per cent. 

During the period of high flow (8:28 to 8:58) it was found that the 

critical depth in the throat exceeded the depth of the Venturi flume. These 
- yesults are shown in Table 4. In this second group the discrepancy was 
~ +6% + 0.7%, showing that even when poorly adapted it was still compara- 
~ tively correct. The 4% error in the first group could be corrected by allowing 
- a head loss of 0.25 Hy, or by subtracting 0.75 Hy, instead of Hy, from the 
- energy head. Doubtless the discrepancy could be reduced materially by care- 
fully installing a permanent concrete Venturi flume instead of the temporary 
wooden one which was used. Any obstruction placed in this sewer laid on a 
slope greater than the critical causes the water to jump immediately to the 
conjugate depth on the energy-head curve (Fig. 1) and unless allowance is 
made for this effect, the capacity of the sewer would be curtailed seriously. 
Therefore, it was necessary to use a large throat with a consequent rather 
high velocity of approach. It has been found that the slower the velocity 
of approach the greater the accuracy. This is limited by the capacity of the 
conduit and the possibility of forming deposits in the approach channel which 


would affect the rating curve. 
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~ The proper design of a Venturi flume requires that parallel flow occur 


in the channel” above the flume and in the throat. The necessary critical 


velocity is obtained only when a drop in the energy head occurs just below 
the throat. A small jump in the water surface at this point (see Fig. 1) 
is positive evidence that critical velocity is occurring. In addition, the 
throat section must have sufficient length or the water will not be flowing in 
parallel filaments at the point of critical depth. The first meter stations 
installed by the Sanitation Districts of Los Angeles (Calif.) County had 
short throats which gave good results only on the small flows. Lengthening 
the throat to 3 ft insured parallel flow through more of the length, and tests 
with a Pitot tube showed no apparent change in energy head except at the 


ends. Data are lacking to formulate a rule as to the proper length, but experi- 


ence indicates that the throat should be at least as long as the diameter of 
the pipe. A comparison of flow through two Venturi flumes in consecutive 
manholes, with 36-in. and 6-in. throat lengths, respectively, showed a ratio 
of 1.0 to 0.8, the longer throat indicating the larger flow. 

The important features of the Venturi flume are its adaptation to stream- 
flow measurement in all shapes of conduits not flowing under pressure and 
its ready installation in lines to which access is much restricted. The device 
can always be installed in an existing line and often without serious interrup- 
tion to the flow. 

The size of the flume throat depends upon the size and grade of the con- 
duit and the range of flow it is desired to measure. The ideal flume throat 
would have such a size and shape that the ratio of the cross-section of the 
water in the throat to that in the conduit would be the same for all quantities ; 
but such a throat would be impossible to design and build because the ratio 
of the critical depth to the depth in the conduit is not a constant. 

For a flow that is only a small percentage of the designed capacity of a 
conduit, a throat, V-shaped, or narrow at-the bottom, with sloping sides, is 
satisfactory. However, such side walls are difficult to hold in place. A rec- 
tangular narrow throat may contract the flow so much as to give inadequate 
capacities for larger flows. 

For circular pipes it has generally been found advisable to install a bottom 
slab in addition to the two sides, a feature which helps to support the side 
walls. ‘The slab acts as a broad-crested weir when flows are only a small 
percentage of the designed capacity. Since a thin bottom slab properly 
stream-lined by approach transitions causes no deposition of solids except 
possibly at the very lowest velocities, such a slab is recommended, especially 
on the lighter grades, to offset any possible chance of back-water from an 
unforeseen obstruction below that might otherwise prevent critical velocity 
from being attained. In general, the flatter the slope, the more necessary is 
this precaution and the greater the slab thickness needed. 

A throat of rectangular cross-section with proper transitions may be adapt- 
able in some cases, as in the instance of a semi-elliptical section, or a rec- 
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a3 “Hydraulics of Open Channels”, by B. A. Bakhmeteff, M. Am. Soc. C. B., 1932, p. 28. 
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the placing of forms for building the side walls, a simple device t 
eate critical velocity consists of placing on the invert a pre-cast flat slab 
with transitions. Deposition can only occur above this type at extremely low 
velocities. In similar cases involving large flows critical velocity may be ; 
pee by simply inserting vertical side walls with the proper transitions. 
The important factor in the construction and installation of any form of . 

‘ Venturi device is that a constriction of some sort be placed in the channel — 

? to. produce critical velocity with the least loss of energy, and that the shape, 

_ size, and dimensions of the device are important only in so far as they meet 

the specific problem at hand in a practical manner. 

. The correct size and shape of throat for use in a given conduit is that one 
for which the energy head is greater than the normal energy head in the 
free flowing conduit at all values of Q, and with a bottom slab thin enough to 
prevent deposition of sludge at low flow. A practical method of making this 
selection is first to prepare quantity-energy head curves for several sizes of 
Venturi flume throats (Fig. 7) on tracing cloth. Table 5 can be used for 


Depth, d, in feet 


0 5 10 15 20 25 30 35 40 45 50 10 15 20 
Rate of Fiow, Q,in cu ft per sec ort of robs Q, in cu ft per sec 


Fic. 7.—QUANTITY-ENERGY HEAD Curve FOR VARIOUS Fic. 8.—ENERGY HprapD 
THROAT SIZHS WITH SipE SLopps oF 2 on 1. CURVE, VARIOUS THROAT 
SIZEs. 


constructing these curves for trapezoidal throats having side slopes of 2 on 1. 
Next, on the same scale, a quantity-energy head curve is drawn for the con- ~~ 
duit under consideration, assuming unobstructed flow. Any formula for 
flow in open channels may be used for first constructing a quantity-depth 


curve, but the following procedure is suggested. Bakhmeteff™ expresses the 
Chezy formula, as follows, 


Q = RP RR rRsat . ee eel ey 


in which, 


K =, AC wh Rb Sing. an eee 


: 
4 “Hydraulics of Open Channels”, by B. A. Bakhmeteff, M. Am. Soc. C. B., 19382, p. 13. 
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= 0.018 and S = 0.01 in Kutter’s formula. Drawn to a logarithmic % 
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Rate of | ENERGY Hukp, in Furr, ror Borrom | Rate of | ENerey Hap, 1n Freez, ror Borrom — 24 : 


flow, Q, ‘ IDTHS, IN INCHES i a 
Be in we flow, Q Witus, in INCHES 
e - feet per 
¥ second 
om ¢. 0.2 0.232 | 0.198 | 0.176 | 0.156 | 0.121 10. 1.941 | 1.820 | 1.719 | 1.623 | 1.388 
a 0.5 0.401 | 0.348 | 0.311 | 0.280 | 0.219 12.00 2 AL, 3 | 1.887 | 1.786 
2h. 1.0 0.594 | 0.526 | 0.474 | 0.481 | 0.342 14.0 2.276 | 2.149 | 2.041 | 1.936 | 1.679 
Ad 2.0 0.865 | 0.779 | 0.713 | 0.654 | 0.529 16.0 2.422 | 2.294 | 2.183 | 2.075 
f.. 3.0 1.070 | 0.975 | 0.897 | 0.829 | 0.680 | 18.0 2.558 | 2.428 | 2.315 | 2.205 | 1.930 
4.0 £2939 1.1367] 12 9 OFS) O-8i0) 1, (2050 Aa gts.) fete 2.439 | 2.3827 | 2.045 
ra 5.0 .384 | 1.279 | 1.191 TOS 3) O2O96 Wah. O cay kites.) lee te | cated ate 2.604 | 2.308 
ee 6.0 1.515 406 | 1.314 | 1.228) 1: Si el ome Oh ee tS Shahi Teel | TRS S 2.852 | 2.544 
J 120 .634 | 1.522 | 1.426 | 1.337 | 1.130 AS a A Se | LAR oP aE 9 3. 2.759 
8.0 dash Goss t bso), 4.4358) L220 40 TT. oe etre edt yl pene 
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it is merely necessary to add one-half the logarithm of the slope (expressed 
as a percentage) to obtain the quantity for any depth in the pipe. Velocity 
heads are obtained from Fig. 4 in the manner previously described. Adding 
these velocity heads to the depths shown on the quantity-depth curve gives 
the quantity-energy head curve. Fig. 8 shows a typical diagram. = 
Referring to Fig. 9, the capacity, Q, of a pipe is given by Equation (20) 
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when S = the slope, expressed as a percentage; K = ATX _ a/ BR for any 
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given depth; and C = Chezy’s coefficient. The curves are drawn for a1% slope. 
For values of S less than 1% the flow for any depth will be to the left of 
the curve and for more than 1% to the right. To find the value of ASS 
measure from the slope guide to the right 10-line for S<1% and to the left : 
for S>1 per cent. For example: Let S = 0.25%; D = 8-in. pipe; and, — s 
Q = 0.1 cu ft per sec. In this case S <1per cent. Set the dividers for the 
% distance between the inclined slope line in percentage line and the right 
~ _1.0-line of the diagram at S = 0.25. Transfer this distance to the Q = 0.1- 

line, measuring to the right and intersect the 8-in. curve at d = 0.20 ft. 

Tf Q is given in million gallons daily, set the dividers for the distance between 
a the inclined line and the million-gallon-daily guide line on the extreme right. | 
For S =-0.25% and Q = 0.1 mgd, d = 0.25 ft in an 8-in. pipe. 

On the other hand, let S = 2.25; D = 60-in. pipe; and Q = 15 cu ft 

per sec. In this case S > 1 per cent. Set the dividers for the distance 

from the inclined slope line to the 0.1-line to the left at S = 2.25. Transfer 

this distance to the Q = 15-line, measuring to the left and intersect the 

D = 60-in. line at d = 0.66 ft. 

By superimposing the quantity-energy head curves for the throat upon 

fi that for the pipe and then moving the throat tracing up or down, it is possible 
i to decide quickly which throat size should be used and how high it should 
f be set above the invert. The difference in the two energy heads shows the 


. oY Si © 


X 
\ 


SAWCT 
A 


SsiSPAla as 
15 -- 20, 
Ss 


6.0 80 10 
Values 


Bed Le Seaale 
pa ae 
=e 

heer 
Baleares 
ag 
SA Zs 
— 


06 O08 1.0 


30 =640 


of 


3.0 4.0 


LS. <2.0 


Fie. 9.—DiscHarcr or Preps with WATmR FLow1ne av Various Deprags WHEN Korrar’s n = 0.013. MARE. tes 
Se EEF 5 BS 


i Seh, a 


S=0.110% 


- Depth, d, in feet Above Invert 


0 
0 5 10 15 3 20 
Rate of Flow, Q,in cu ft per sec 
g Fic. 10.—ComPARATIVE ENERGY HEAD IN TRAPBZOIDAL AND RECTANGULAR FLUMES. 
; and it is easily seen that the trapezoidal section can be used as a meter for 


flows as great as 18 cu ft per sec and that at 19 cu ft per sec, which is the 
maximum capacity of the pipe (there being no loss in energy through 
the flume), the carrying capacity of the pipe is not curtailed. The quantity- 
energy head curve of the rectangular flume (Fig. 10) being always above 
the energy head for the pipe, shows loss of energy through the flume at maxi- 
mum designed capacity. This type of Venturi flume would interfere with 
the use of the sewer at full capacity. It is more accurate to compute the rat- 


as peraanent ‘astallation 3 in sewers of the ‘Lae ner acy County's Sanitarian ts 
D istricts the Venturi flumes were readily constructed of concrete iis. in 
; sewer manhole without interruption to sewage flow. A layout sketch of the © 
meter was prepared for a single. installation in a circular pipe and, using “ 
the plan as a pattern, forms for other sizes were obtained by proportion. When a 
5 i possible the entire form for the flume was assembled in the sewer manhole ~ 
Ns Be \ ‘invert and concrete was poured into the forms from one side until the bottom 
slab was filled in order. to assure complete displacement of the sewage. Use 
of a fast-setting cement had distinct advantages in such installations. 
Placing the center of the throat at about the lower side of a manhole gave - 
proper space for installing the stage-recorder float near the upper side at a 
point where the flow in the pipe was uniform. The float consists of an 8-in. 
copper ball fastened to a 2-ft length of 4-in. pipe, hinged loosely up stream 
to a pipe framework (see Fig. 1). The latter is attached to the top of a 
plank placed cross-wise to the sewer and supported by the shelves of the man- 
hole. Without some such apparatus to hold the ball, it tends to float down 
. stream and disturb the record. The stagé recorder is suspended on a spe- 
Sw cially built steel seat immediately under the manhole cover where it is readily 
accessible. Repeated checks of indicated stage records against actual meas- 
ured depths show no appreciable error during the daily range of velocities. 
Braided, insulated, copper wire connects the float to the.recorder. Consider- 
able difficulty was experienced in keeping the stage recorder clocks running 
for a week at a time in the damp, gas-laden, sewer air, but the purchase of 
sealed clocks has practically solved this problem. 
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THE STRESS FUNCTION AND PHOTO- 
ELASTICITY APPLIED TO DAMS — 


By JOHN H. A. BRAHTZ,* Esa. 


§ : Synopsis : 
--—-s- The object of Part I of this paper is to familiarize engineers with the 
use of the Airy stress function for the solution of problems in plane. stress 
and plane strain when ordinary engineering methods fail to give even approxi- 
mate results. : 
The object of Part II is to familiarize the engineer with the photo-elastie 
phenomenon and its application to civil engineering structures and to com- 
pare results obtained in this way with the theoretical results obtained in 
Part I. 
A brief popular outline of the theory is given accompanied by a descrip- 
tion of the apparatus developed at the California Institute of Technology, 
at Pasadena, Calif. Methods of evaluating stresses are explained and applied 
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wey 
ie 


to experiments on Morris Dam. 


¥ INTRODUCTION 


Ordinarily, the theory of the Airy function found in textbooks? is based 
on the assumption that body forces (weight and inertia) can be neglected. 
. This is entirely inadequate in civil engineering structures in which the 
i stresses due to weight often are greater than those caused by boundary forces. 
Furthermore, in texts in which the body forces are included the definitions of 
_stresses are generally such that the form of the function is not invariant to 
a change in co-ordinates. The stress definitions given herein, including both 
boundary and body forces, are such that the stress function will be of the 


same form in rectangular and 
er wifl be closed in December, 1935, Proceedings. 


polar co-ordinates.° 


Nory.—Discussion on this pap 
1 With the U. S. Bureau of Reclamation, Denver, Colo. 
2 See, for example, “Applied Hlasticity”’, by J. Prescott. 
3 “Notes on the Airy Stress Function”, by John H. A. Brahtz, Bulletin, Am. Math. Soc., 


June, 1934. 


“stresses: ie is aay, conceded that-a. slight variation Ce elas 

odulus has little effect on the final stress distribution, except, of cour: 

‘singular points or at points of high stress concentration. With these 

assumptions the stress function may then be applied to a slice of a gravity 

dam. The question of uplift is not considered in this treatment, but it is 
assumed: (1) That there is sufficient resultant average compression at all: 
points of the dam to overcome any internal pore pressures that may exist; 
and (2) that the pores are so small that the average stress distribution ata 

- point may still be found as in isotropic material. 

_ This paper, together with theoretical and experimental work done by — 
others, shows that for purposes of analysis the triangular gravity dam on an 
elastic foundation may be divided conveniently into. three regions: (1) The 
‘upper two-thirds of the dam proper; (2) the lower one-third, including 

the base region of the foundation; and (8) the foundation proper. Finally, 

- special investigation must be made in the regions close to the heel and toe 
as to whether they are sharp corners or fillets, and near the crown. Conse- 
quently, in the present application, stress functions have been derived for 
four cases and a twofold purpose is served: (a) To show the methods of 
deriving stress functions; and (6) to obtain specific results applicable to the 
gravity dam on an elastic foundation. In most cases the derivations have 


: : been omitted, due to the limitation of space. The original manuscript is on 
aay file in Engineering Societies’ Library, in New York, N. Y., and at the Cali- 

; fornia Institute of Technology, at Pasadena, Calif. 

is Application I—The stress functions, stresses, and deflections valid in 
4 the upper part of triangular dams are derived for hydrostatic, and are given { 
* for body, forces with computed results plotted in the case of Morris Dam. 
; In addition, the stress functions are given for a number of special loadings. i 


Application II—The stress functions applicable in the foundation are 
given with stresses and deflections for concentrated and distributed loads, the 
computed results being plotted for a study of Grand Coulee Dam. 

Application III—The “corner-function” applicable at sharp re-entrant 
corners is derived. Two methods of procedure are given for the determina- 
tion of stresses in the lower part of gravity dams. The results are plotted ~ 
for a study of Grand Coulee Dam. 

Application IV.—An approximate method is derived for the determination . 
of stresses at re-entrant sharp and rounded corners. Examples are computed 
for Morris Dam. 

It should be emphasized that even if the gravity dam is used for illustra- 

_ tive examples of application, the aim is not to advance new design criteria. 
The writer hopes, however, that the methods described will help to obtain a 
closer estimate of the stresses that actually occur. It is worthy of note ‘that 
a state of stress computed by the Airy stress function is in equilibrium and 
compatible with Hooke’s generalized law. 


steman, of the California Institute of Technology, P 

i valuable suggestions and co-operation in the preparati 
per. The photo-elastic experiments on Morris Dam were made pos 
the financial support of the Pasadena Water Department. 
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- PART I—THE THEORY OF THE AIRY STRESS FUNCTION — 
4 In this section the stress function is defined in rectangular and moles 
co-ordinates. Convenient forms of boundary conditions are treated. 


__ A two-dimensional elastic system under plane stress or plane strain is in 


brium if the stress components are defined as follows: : 
By rectangular ¢o-brdinates: . eat. 


_— ° >F ; " ; 
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in which gz and oy = components of normal stress parallel to the X-axis and 
_ the Y-axis, respectively; tz,y = shear stress in the direction of the X-axis 
oe or the Y-axis; F = a stress function; and g = the total body force per unit 


4 volume. 
fe In the polar co-ordinates: 
fr Te 
af 2 
5 ¢ = ae bc ORs — gr cos (6 — B)....++- Or Coa) ; 
Ha 7" 08; r or : va 
J oF rs 
OF — g— Rraleneleh s''s t olleleké © -s\v. ele .(26 » 
5 ; 8 a gr cos ( B) (2b) 2 
and, 
3 nae 2 (2 SERS MO NE Cc) 
a or \roé 
4 If, now, F is so restricted that it is a solution to the differential equation: 
4 >t 4 
Sef Pe BS as I a bt BG) 
ox! dz? ay? sé! 
or, 

alps (DOr hog 2 oie 28 ) exighiwe, 20. Pao (3b) 
or? r or r? 06? 
; the stress defined in Equations (1) and (2) will also be compatible with the 
; tress or plane strain have thus 


generalized Hooke’s law. Problems in plane s 
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PsA @ + ae @) + vor aep + Aah ee Same ie) 


i064 Tits 


Pca Equations (3). rae Aa, ahs and A, are arbitrary constants; de fo fs z 


and f, are arbitrary analytic functions of the complex variable, z = 2 gfe ays 


Ped t ys andi= ai. pa 


The second part of the problem (to satiety the noueelaes soi distin) is” 
usually difficult, and must be solved for each individual case. It is important ; 
to realize that as soon as the stress function is known the stresses are also. 


(ees at all points by Equations @.). or Equations (2). 
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Bounpary ConpiITions 


be etareinic to Fig. 1 the boundary conditions to be satisfied by F (xz, y) may 


be formulated in terms of well-known engineering concepts, as follows: Let 


¥ 
BE 


Fie, 1. 


the area enclosed by the hatched boundary curve represent a plane elastic 
body of unit thickness, held in equilibrium by a system of forces, P, (con- 
centrated or continuously distributed), acting on the boundary, ABC, and 
body forces, g, per unit volume; ge and Jy = components of g along the co-or- 


dinate axes; and H and V = the forces along the Y-axis and Y-axis per unit 
length of the boundary. 
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Se. ae faci 
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ecting all forces on the co-or 


a point, B, of the boundary with co-ordinates, Le Yn, the con 


ditions 
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gy Ge Se dy = faye gg tnitiorth COTE) pateiaey 


Sg MOR ee dy — cdx) (ys — y) + (oy'dx — tdy) (az —*x) = dM... - (5c) a 
i 1 which dM indicates the moment of Hds and V ds about Point B, posi- 
_ tive in direction (counter-clockwise about Point O, Fig. 1). The body forces _ 
a of the element are of higher order of magnitude and need not be considered. 
By substituting the stresses as defined by Equations (1) and integrating the 
equations from Point A, Fig. 1, along the boundary in the positive direction 


to Point B, the following equations are obtained: 


B B 
: (=) —(2)-x+0 fewto f Udy ser (6a) 
oy /B OUe fos A A 
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(2) -(=) aa¥te [edt om f y de =.:40b) 
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E in which X and Y are the projections along the X -axis and Y-axis; and 
_ M = the moment about B of all boundary forces between A and B. By a 
~ suitable choice of co-ordinate axes Equations (6) can be simplified somewhat. 
From the definition of the stresses in Equations (1), ete., it follows that a 
constant and terms that are linear in x and. y, added to the stress function, 


contribute no stress. Therefore, by choosing the origin at Point A, Fig. 1, 

ih : oF or 

it is always possible to arrange matters so that F, = (—]) =(—) =92. 
Ox / 4 oy J 4 


The boundary equations then become: 


B 
(=) = X,+ a f x dy + i ee (Ta) 
w/z « 0 2 


: B 
Be (=) = — Yg+ i 92 UB + wf YL osscrecees (7b) 
J or / B 2 ° 


a a Ae 


Eee are to be sea out from the origin Godlee 0, "Tig 1) in 
tive (X-Y) direction. Hence, if Point B is located in the negat’ 
‘ection from the origin, the signs of Xz, Ys, and M zy as defined in conn 


n with Equations (6), must be reversed, exactly as in taking moments fa 5 


hrust, and shear at a section of a beam in bending. This formulation of | 


the boundary conditions will be found: especially useful when approximate 


solutions are desired; that is, in case the boundary conditions are to be 


"satisfied only at a few points. a cee 
a In the cases of exact solutions (that is, if the boundary Lee iieana can be 


"satisfied at all point of the boundary), only two of the three sets of formulas 
Equations (5), (6), and (7) become especially simple and instructive. In 


ats be satisfied automatically. If no body forces exist (that is, if gz = gy = 0) 


Equations (5), (6), and (7) become especially simple and instructive. In 
this case the following simple interpretations of the stress function, F'>, and 


its partial derivatives, “en and oF ; are evident ee Equations (7 }e 
M7 


‘The value of Fp (the stress function for the given boundary forces) at any 
point, B, of the boundary is equal to the moment of all forces acting on the 


boundary between the origin and Point B. The value of oF r at Point B is 


equal to the projection on the X-axis of the same forces and oF p equals 
ox 
the negative projection on the Y-axis. 


It is often convenient to determine two functions, Fp and Fy, the first 


representing the boundary forces and the second the body forces; then, 
Fr = Fp + Fy. If no mass forces exist the boundary conditions can be given 
a simple form on the portions of the boundary where no forces are acting, 
namely: 
OW Ae — OR By nent meter oni e Smeg (8a) 

and, the derivative normal to the boundary, 
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aot FDS 5. Sas Oo ccm Ate aes ae em 
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When the stress function, F, is known, the displacements, wu, parallel to the 


X-axis, and v, parallel to the Y-axis, can be determined in the case of plane 
stress by the expressions: 


La 


7” oF 
Bu = =O Wp Get ory + 20) +htAy + B..(9a) 


a ee 


ee Te 


4 Sat a 4 r xg Smit +. 
VARS 5 y} 
rs hes 


we SE oust — yt + 2920) th Ae 


ch, E is Young’s modulus; p is Poisson’s ratio; f, = real 


ae 


4 
fe : a ; RS a 
f(z) dz; fe = imaginary part of she f(z) dz; and, f(z) is the analytic 


tion of the complex variable: 


2=at+ ty = PR OGES ha cour eee 


_ such that v* F = the real part of f(z). The constants A, B, and 7 in, Foun ; 
tions (9) determine the reference axis for u and v and only effect a solid- 
_ body rotation and translation. ¥; e 
In polar co-ordinate’ the radial and tangential displacements, ur and vp», aa 


,; are given by: 
pal ‘ 


Buy, = —(1'+ p) or Star os (6 — B) + f1+ Boos + Csin @.(11) 


a and, : 
Eu =—(1+ pz) a - = gr?sin (6: By + fat Air 
3 $C cos 0 — Bisin OO OR 0A) CA 


in which f = real part of e” te f(z) dz; f’2 = imaginary part of e? 4 f(z) dz < 


and, f(z) is defined as in connection with Equations (9). In Equations 


* (11) and (12), the constants, A, B, and C, have the same meaning as in Equa- — be 
a tions (9). If no mass forces occur, g = 0. Displacements in the case of 
7 plane strain are found by Equations (9), (11), and (12), by replacing H with me f 
ana p with —Y—. ze 
1 1—up pape 4 


ad Appiication I.—SrTrESsES IN THE Upper Part or TRIANGULAR Gravity Dams 


In order to indicate the general procedure for determining the stress func- 
tion in the case of two-dimensional structures for a given loading the case 
of the infinite wedge with hydrostatic pressure on one side, and no body 
forces, will be developed in some detail. The solution was first offered by 
_ Mz. Levy.* 

4 Assume the hydrostatic pressure on the up-stream face, y = 0, to be p at 
a unit distance from the top measured along the face; then, by Equation (11) 
with y = 0, the boundary conditions are (remembering that no mass forces 


occur.so that gz = gy = 0): 


ae EE eS ee 
_* Comptes Rendus, Vol. 127, 1908, pp. 10-15. 
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Fy = Ad’ + By’ + Ory + Dey. oie 


. 5 pattntien Equations (3), and, therefore, F, is an Airy stress function. Tf the 
2 arbitrary constants A, B, C, and D, can be determined such that the boundary — 
- conditions are satisfied, Equation (16) is the solution. 
} By substitution of Fp into Equations (13) and (14) and with y = 0: 
27) 6 Az = = — px. Consequently, ay ik —2 Cx = 0; and, i= ae th Bs suis 
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1 stitution of Fy into Equations (15) and with y = Kz; 6 Bak? + A Dek Kee 0; 
_= +2D2K=0;B=——? ;andD= tapes 
3 K* 


ie these values of A, B, C, and D, chokes into Equation (16) the 
stress function for the aban forces becomes: 


e a ko = = aa (fK° + 2y° — 8ay’K).. see 2D) 

_ By the definitions expressed as Equations (1), with gz = gy = 0, the 

nf __ stresses are: ; 
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The displacements corresponding to Equation (17) are derived in some _ 


detail in order to illustrate ae use of the general equations (9); thus: 
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E quation (19); and, f. = azy — 3 (a? — y*) = imaginary part ofthat 7 
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In Equations (20) #, and p, are Young’s modulus and Poisson’s ratio, 
respectively, for the dam in case of plane stress (cz = 0). In case of plane 
strain, HF, and py are found, as explained in connection with Equations (11) 
and (12). ‘aes 
~ -- Mass Forces.—Let the components of the mass forces (weight and earth- 
quake forces) be ga and gy along the X-axis and Y-axis, and assume that no 


4 forces are acting on the faces of the dam. It may be verified that the 
_ function: 
“ Py = —& [ge K*x' + (ge K + gy K* — 2 gv) + 39K ay’].. (21) 


2 6K? 
will satisfy the boundary conditions expressed by Equations (5) or Equa- 
’ tions (7). By Equations (1) the stresses are: 
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Sele Equations (9). 
‘It will be seen by Equations (18) and (22) that all stresses are linear ies, 
any straight line in the triangular dam of infinite height. This agrees with 
the assumptions of engineering practice. It will be shown subsequently that 
this is not true near the base in a dam of finite height. 
The boundary stresses computed by Equations (17) and (21) for Moe 


Dam, in California, are plotted in Fig. 2. 
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Additional Solutions in Polar Co-Ordinates® for the Triangular Dam of 
Infinite Height——Simple expressions are obtained when the bisector of the top 
angle, y = 2a, is chosen as the X-axis, and the Y-axis is positive down 
stream; as, for example: 

(1) When the full hydrostatic pressure, pr, is exerted on the up-stream 
face, and with no forces on the down+stream face: 


+3 
r= val (3 sin? a — sin? @) sin 0 — (8 cos® « — cos? @) cos 8 
24 sin’ « cos® « 


5For greater detail refer to Carothers’ “Plain Strains in a Wedge, with Applications 
for Masonry Dams’, Proceedings, Royal Soc., Bdinburgh, Vol. 58 (1913), pp. 292-303. 
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oy nae pr? (sin 20 — 26 cos 2a) 
= _ . 2 (sin 22 — 2a cos 22) 


r 


moe (4) When a uniform pressure is exerted on the up- Patron face, at no s 
Pierce on the down-stream face: Wage tncae 


Z pad a Sin 26 — sin 2x — (20 — 20) cos 20 
4 sin 2a — 2a cos 2a 


6) With a pressure, pr, on the up-stream face and tension, Pie on oo x 
_ down-stream face: > Sehig 

_ Pr cosa sin 30 — 3 cos 3a sin 0 (27) 

4 6 / cos a sin 82 — 3 cos 3a Sin « ee | 

a (6) With a pressure, pr, on both faces use Equation (27), ery 

“gin” for “cos” and “cos” for “sin” throughout. : 

- (7) With a moment, M, per unit of the Z-axis (positive in direction +0) 
acting at the vertex: 


“ tg 


M sin 20 — 28 cos 20 Se eo 
2 sin 2a — 2a cos 2¢ 
(8) With a force, P, per unit length of the Z-axis acting at the vertex 


a : 4 F=— 


- and directed along the positive Y-axis: 

: pep De at ee Ae aes sp 7 
a 2a — sin 2a 

3 (6) With a pressure, pr, on’ both faces use Equation (27) pubseinayae 


directed along the positive X-axis, use Equation (29), substituting “sin 6” — 
for “cos 6”; and 


e (10) With body forces, g, per unit volume, directed at an angle, 8, with 1M 
_ the X-axis: 
; 
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The stresses are obtained by Equations (2) with g = 0 if no mass forces exist. 


Appiication Il.—Srresses in Founpations.” 
Let a tangential line load, Q, act per unit thickness of the infinite half 
plane (foundation). With the co-ordinate system in Fig. 3 it is easily verified 
PLO NS Se a ae SS aS Oe 


Pets Ge) Se eee 
‘ 6 Wor a more detailed solution, see “Applications Potentiels 4 l’Htude de Hquilibre et du 
Mouvements des Solides Blastiques’’, by J. Bouissinesq, Paris, 1885. 
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62 -%. Therefore, by Equations (7) no tangential force distribution occurs 
_ when r > 0, or when r < 0, but a force, +@Q, occurs at r = 0 because here 


makes a jump equal to +Q. 


The stresses are found by Equations (2) with g = 0; thus: 
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The rectangular stress components in the foundation may ie found by the 
formulas: 
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Fo=is pipet USo se eet 
Tv tT— 


in ahich the angles, a, and a2, are as shown in Fig. 3. Equations (84) are 


§ tions (34) becomes invalid; ¢ = or is then obtained by direct integration of | 
_ Equations (31); thus: APs 


ey = = ee [+ cos 2a, — cos Bait 


te, = oa ne (2a, — sin 20,) — (2a, — sin 2a,)]..... 


easily extended to a number of uniform distributions, q, q, qs, etc. 
At Point a of the surface, cy = + = 0 (see Fig. 2(b)) and og by nia 


maw ath atl 
ane f cdr == 2h poe 2% pg 2 ape 
a Wt o r T a van ae 
in which minus is to be used if a is on the positive side of the distribu- 
tion, q. eth: 
' By Equations (9) (remembering that g = 0) the displacements due to 


Equation (33) are: 


a + pg) 2 2)" _ + log [e — 2)? + y' — Ay + ch. 38 


(z—27)?+y¥ 
and, 
p= a+u[ be Petar L] 
aA "@—aety © — 2, 
: — 2 tant —Y +42 +B} POM ph th Satie (36b) 
r— ty, 


Next, let a normal compression, P, act per unit thickness. The stress 
: F i P z= P ae | : F 
function for this load, F = + —r@cos#? = + — a tan” [=], will satisfy 
‘ Tv Tv x 
the boundary conditions. The stresses by Equations (2), with g = 0, are: 
2 P sin 0 


ees gy = 0; and + =.0. 
ese 
The rectangular components of stress due to P are: | 
ra: 2P sin @co? 6. _ 2P sin’ ae __ 2P sin? 6 cos 0 
ee Sy SS GRE; Ogee 
T Tr ieee T T 
The influence function, 
F = de (x ers cy) tan ( y ) Ree cr akeicieJohaienploreeum (87) 
{ar L— ay 
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Toy = Epes [+ cos 2a, — cos 2a]. We) 
| ae Ls dole aie 
We Faaciite ee can He extended to ice sin ited distributions Ds Py Po 
oe ‘The displacements due to Equation (87) are: — 


a | 1 eS ee = -1 ] oe . 
= 1 eS = hy es R 
‘a uw E, {« sc ls (@—2r+y¥ j eerie 
ere +2 tant —¥ + ay +t 

; t— Ty ‘ f : 
and ae lend: 


~ 
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An interesting case occurs when po is the uniform reservoir pressure on 
the foundation above a dam. If the X-axis is chosen positive down stream, the 
stresses are: 


= a [7 650] "sin 2 aa] So eee (39a) 
2 
oy = ee Pe + Do, + Bim 9a]... ~ cs ee ee ROE 


and, 


tay = —22 [1 + cos 2 Gi}. va suns see seen Oe 
20 


The corresponding stress function is F = + we (sin 2 6 — 2 6), 
with the origin at the heel. The corresponding displacements are: 


u= 


- {a + ps) x tant +22 tant” +— y log (x? + y*) + Ay+ at a (40a) 


2 
and, 


= Po {0 + pe) y tant¥ — 9 y tant?Y 4+ x log (2? + y?) — Ar + c . . (40b) 
ae oF £ 2 ‘ 
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Finally, the stresses in the foundation due to the weight, g, per unit volume 
are: oe = oy = — gy; and t,y = 0. 


at i a ic rt i ii eis fe 


ey, 


\ 


ASS *. 
* ba hs NY 


Tay As 


tions” for a wedge or a corner with an a 
tion is defined such as to give no forces on the straight bour 
: or the corner. For the present purpose it is convenient to. 
co-ordinate system shown in Fig. 4. These functions appear in conne 


tion with the exact solutions by the writer’ for stresses in two-dimensional 


corners. H. M. Westergaard, M. Am. Soc. C. E., in a communication to 


Professor yon KArman, has called attention to the existence of such fune- 


tions, which he used for the investigation of stresses in the region of the toe 


and heel of dams, cracks in concrete, ete. 
The definition of this function is expressed by Equations (8); hence, when 


5 6-0 and 0 = y: ° 


5 ta aN aS i Man toa mele Ss) 


‘ on r 06 
Consider the real part of the following function which, by Equation (4), 
is known to be a solution of Equations (3): 


OS ee Tas ok rhe Se 


F=[(B—iA) cr +(C—-iD)y]) A cane eaten hd 


eatwhich x =) cos 62 y = rsin 6; 2.= re? = r (cos 0 + 4 sin 6); and, 
= re"? =r (cosnd+isinn 6). By substitution of these expressions 


into Equation (42) and taking the real part only: 
F =r [Bcosn 6 cos 6 + Asin n 6 cos 8 + CO cos n 6 sin 0 
+ D sin n 6 sin i crise pb, Sea oP) 


T Presented to the California Institute of Technology in 1932, in partial fulfillment 
of the requirement for the degree of Doctor of Philosophy; also contained in Technical 
Memorandum No. 420, U. S. Bureau of Reclamation, Denver, Colo. ; excerpt published in 
Applied Mechanics, April-June, 1933, Vol. 1, No. 2; and in Physics, Vol. 4, No. 2, 


February, 1933. 


eae Ate [n cos n 9 sin 8 oh msinn@sin@ — — sin ng cos 0]... .(44) 
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. ahs The’ Siledeae produced by Equation (44) are found by Equations (3) with — 
A g = 0. The displacements corresponding to Equation (44) are: 


te =F hatn) [n cos #-sin (n—1) 6 + sin n 6 — n’ sin @ cos (n—1) 0 


4} 
— mnsin 6 sin (n— 1) 6] — 2sinn @ —2nsinn@ 


Sh omnia Bhar glo A ce 


a ja +H [n cos 6 cos (n — 1) 0 —necosn# — msinn #6 
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+ n> sin @ sin (n — 1) 6 — mnsin @ cos (n — 1) 6] 


Bgl st +2cosn0+2necosn0+2msinn@—Creos@+ HH 


for wand v. The values of n that satisfy Equation (45) will be referred to as “cor- 
ner-values.” It will be found that when t < y < 2 x, there are two significant real 


i in which C, D, and F are arbitrary constants used to fix the reference system | 
i 
roots, No and 7, the others are complex; when y = ¢, all roots are real: 1,2,3,ete.; ; 


when y = 2 7, all roots are real: =, 1 > 2, ete.; and, when 0 < y < 2, all roots 


are Pemnplex: The real roots, mo and m, for various values of y are given in 
Fig. 5, based on the formulas: 


BIN No Y Se No BIR. t,o. ok. vt ine Rae ETE 
and, 3 
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The complex roots may be found as follows: Let n = a + ib, which, substi- 
tuted in Equation (45) gives: 


+ (a + ab) sin y = sina y cosh b y + icosaysinh by 


= 
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j - waa} 
; r ex > eee 
‘ — ? hn . » 4 3 wh 5 
eae yo < * Fis Joe 
t y ‘ 
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Ayia 


cosa y sinh by = + Bsin yee 


Values of ¥ 
Fie. 5. } 


When y is known, approximate values for a and b are found by: 


shia (A ath): Sete RL DOr EAS 


aa = 
27 
and, 
b= tl log {CED sin es eae 
Y y | i 
in which, 


ih 
a (2k —1) 


Pine (5 (2k — 1) sin x) 


&, 


and, k = 2, 3, 4, 5, ete. 

The values, a, and bx, obtained by Equations (49a) and (49b) will not 
satisfy Equations (48) exactly; therefore, they must be further adjusted.° 

When nz is known the corresponding value of mz is found by Equa- 
tion (46). Let m= d+ te which, substituted into Equation (46), separating 
real and imaginary parts and making use of Equations (48), gives e 
=btanay; and d = cot y — a cota y. 

The various quantities in Equations (44) have now been expressed in . 
terms of y so that, when the angle, y, is known, the “eorner-function” is fully 
determined except for the arbitrary constant, A’. Each complex value of n 


will furnish a real and imaginary part of Equation (44). Both parts will 
satisfy the required boundary conditions; therefore, two functions are obtained 


8 Adjustment demonstrated in record manus 
Societies Library. New York, N. Y. 


cript, filed for reference in Engineering 


iA 


0, 1, 2, 3, ete. 
r Shieh Fy, is given by Equation (44). 
The function, F y, contains a double infinity of arbitrary constants, Ax and | 


forces acting on the straight boundaries of the corner, y, irrespective of the 


_ values of Ax and B;. At points inside the corner the stresses will not be 
_ zero, but the stress distribution will be such that the resultant of the stresses 
acting on any line terminating at points of the boundaries will be zero. In 
other words the function, F',, will only affect the distribution of the stresses. 
It is important that the full significance of F is realized. It is extremely 


useful for the investigation of stresses at sharp corners and for the correc- 
tion of stresses in structures if the boundary lines differ from those assumed 
in the original computations. 


Application of F', for the Determination of Stresses Near the Base of a 
Straight Triangular Gravity Dam on an Infinite Elastic Foundation 


Two cases will serve to demonstrate the foregoing theory. -In Case 1, the 
elastic properties of the dam and foundation are assumed to be the same 
throughout; and, in Case 2, they are different in the dam and foundation, 
but constant throughout in each of the two. 

Case 1—Consider an infinite elastic homogeneous plate of unit thickness. 
On this plate imagine the boundary lines of the cross-section of the dam 
and its foundation. These lines will divide the plane into two separate regions 
one of which is a slice of the dam and its foundation. The problem is to 
determine a function, F, that will deliver stresses between the two regions 
equal to the loads on the dam and its foundation. This function then will 
deliver the correct stresses throughout. Incidentally, it is interesting to note 
that F’ will give the stresses in both regions of the infinite plate. 

Select the height of the dam (see Fig. 2), measured along the up-stream 
face, as the unit of length. Assume a full hydrostatic load on the up- 
stream face and constant pressure, p, on the up-stream foundation. Let the 
uniform mass forces (weight and inertia) be g per unit volume acting at an 
angle, 8, with the X-axis. It is assumed that the forces, if any, acting on the 
rear face and the down-stream foundation are known. 

Then the boundary conditions which must be satisfied by the stress func- 


tion, F, are: (1) When 0 = 0; y = 0; (up-stream face): ¢, = — pA—r); 
and, t = 0; (2) when @ = y; y = 2 tan y = wK; (up-stream foundation) : 
o9 = — p; and +t = 0; and, (8) the rear face and down-stream foundation 


will be considered later. Let, 


P= — Py ROS Fy eRe ae Ae ey 


F;, = real part of 5. ap eee oe 


By, and is constituted so that it corresponds to the case in which there are no — 


EEE 
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eaten +e 4 
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h Fp, Fy, Fo, and F., are given in Equations (17), (21), (24), and (E 
7 is the re-entrant angle at the heel (y > =); and, K = tan y. Ao oat ee a 
2 By adding the stresses when y = 0 and y = Kz, respectively, due to the | 
_ various functions in Equation (51), it will be found that Conditions (1) and. 
3 (2) are satisfied exactly. (It is to be noted that F, was determined ao Mate ; 
_ give no contribution on these lines.) — a ‘ts 
| There remains the matter of satisfying the conditions along the line, — 
_ TBDA, in Fig. 2(b). For this purpose the form of the boundary conditions 
_ derived in Equations (7) is utilized. : - cla 
f In Equation (51), F is determinate with the exception of the arbitrary : 
constants, A; and Bx, in the function, /,. This double infinity of constants 
is used to satisfy Equations (7a) and (7b), at all points of the line, T.BD A. 
3 o The right-hand sides of Equations (7a) and (7b) are completely known and 
ean be computed at all points when the forces on the dam are known. A> 
a double infinity of simultaneous equations is obtained with a double infinity — 
; of unknowns, A, and By. This problem is solved, mathematically, by expand- — 
ing the right-hand sides of Equations (7a) and (7b) in terms of the “corner- 
functions,” Fy. In the practical solution of the problem it is only necessary 
to satisfy Equations (7) for a finite number of points. In this case it is 
“- better to include all three of Equations (7). It will be found that if they 
are satisfied at three points, B, D, and A (Fig. 2(b)), Equations (7) will be 
closely satisfied at intermediate points. This can be checked by plotting 


le mee rad oF 
ony 


, along the line, 7 BDA, and comparing the results with 


the values computed by the right-hand sides of Equations (7). 
Three equations of conditions are thus obtained at each point, B, D, and A, 
_a total of nine equations. This means that nine constants, Ax and Bx, must 
be included in F’y. ; 
When only a finite number of points, B, D; A, etc., are used, it is neces 
sary to compute the stress functions for both the heel and the toe; that is, 
for the two different angles, y. As soon as the constants, A, and By are 
2 determined, F', is known by Equation (50). Thus, the total stress function 
valid in the region of the corner is known by Equation (51) and the stresses 
can be determined by Equations (2). 
The stresses computed in one of the preliminary studies of Grand Coulee 
Dam’ are given in Figs. 6 and 7. The stresses in the foundation were found 
by Equations (34), (85), (38), and (39), after the normal stress and shear 
distributions along the base had been computed. The stresses in the upper 
two-thirds of the dam are not shown. They may be computed by Equa- 
tions (18) and (22), or by the usual engineering methods. 
Tt will be seen that all stresses are compressive and tend to become infinite 
Of course, this is not possible because the material will be- 
magnitudes are reached and redistributions will 
ch must be considered as singular points where 


at the heel and toe. 
come plastic before such 
occur near the corners, whi 


See See EES 
° The details are contained in Technical Memorandum No. 408 by Chief Designing Fingi- 
neer J. L. Savage, U. S. Bureau of Reclamation, Denver, Colo, 
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are v very small. The effect of a redistribution, therefore, can only be slight. 


the theoretical stresses within the elastic limit at all times. This was done 
in other studies of this dam. 
ase 2.—Referring to Fig. 2, consider the dam as a triangle with elastic 
24 tus F,, and Poisson’s ratio, p., resting on an elastic foundation which 
is considered as an infinite half plane with modulus, E., and Poisson’s ratio, ps. 
A slice, one unit wide, across the dam and into the foundation is considered. 
It will be assumed that no separation or slippage will occur between the base 
and the foundation. This implies that, when all effects are considered, the 
, resulting normal stresses are compression, and that the shear stresses are 
i within the permissible relation to the normal stresses. 
it All forces acting on the dam and on the foundations are known, except 
the reactions between the two. The problem is to determine these reactions 
; such that the elastically deformed base and foundation make a perfect fit 
’ when equilibrium exists. This is expressed by demanding that the relative 
displacements of the dam and the foundation must be zero along the base. 
o -. With the co-ordinate system as in Fig. 4, let the stress function for the 
| dam be F = Fp + Fy + F,, corresponding to full hydrostatic pressure on the 
' up-stream face, and the mass forces, g, per unit volume. Functions Fy and Fy 
: are given by Equations (17) and (21); and F’, is the “corner-function” derived 
in Equation (50) with y equal to the top angle of the dam (y < 2). In this 
case F’, may be considered as a correction to the basic function, Fp + Fy, 
for the infinite wedge., It has already been indicated that Fy does not alter 
the foree system which acts on the faces of the dam. Its value is known 
except for the arbitrary constants, Ax and Bx, which must be determined by 
making the deformed base of the dam congruent with the deformed foundation 
line. By this method it is possible to determine: (1) The effect of a varia- 


tion of'the ratio, —, on the stresses, in which FZ, and EZ, are Young’s modulus 
‘ 2 


for the dam and foundation, respectively; (2) the effect of uniform shrinkage 
of the dam relative to the foundation; and (3) it includes Case 1 when 
Ey = Hs and pa = po 


Appuication I[V.—Approximate Bounpary Stresses at SHarp AND Rounpep 
Corners or FILurrs 


Formulas are developed herein for the boundary stresses at sharp and 
rounded corners of elastic structures, using only the real “corner- values,” no 
and m, given in Fig. 5. The method is applicable to any structure with any 
loading. Numerical examples are computed i in the case of a triangular dam. 


distance from the corner. It will be noticed that the theoretical stresses are 
extremely high over only a minute distance so that the actual forces involved _ 


In order to be able to predict the elastic stresses at the heel and toe it” 
~ would be necessary to fillet the corner with a large enough curvature to keep — 


ee ee 


0 of the corner in E 


represented by a function, F'., and that the stresses, ga and gp, at the points, 


a and r = b, of the boundaries are known, or have been comput 


- ordinary engineering methods. It is required to determine a function, F 


_ represent the state of stress at the corner such that, when 6 = 0: 
oe | . : 


i: 


2 Then by Equations (2), 
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3 vat Ee + 420 sg c0s9—8) | tiga dosk ee 


ee ty OF (9=0 
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; q On the boundary, 6 = y, Equations (52) and (53) apply likewise except, of 


Ry 


Sl a 


course, that here 6 = y, r = b, and ga is replaced by ov. Let, 
BT A! oeteacosis tse hs 42 os oie ee ee 


Function F, is the “corner-function” for y given by Equation (50), with 
k = 0 and 1; that is, only the real values, no and m, found by Fig. 5 are 
involved. Hence, Fy = Fo + F,, in which Fo and F, are given by Equa- 
tion (44) for n = no and n = m. 

It will be seen that F defined by Equation (54) satisfies Equation (52). It 
only remains to determine the constants, A, and A;, such that Equation (53) 
will be satisfied. 

By substituting F into this formula and re-arranging: 


Ee 1 2Fy] oe -[2* yer Oe 
roe i OME 7==0 : 770 F T On 
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— Tg cos @— 6) | BAM daicis tore (55) 

; G =O 

. PS! 
and, the analogous equation for 9 = v, r = }, in which every term is known 
except Ao and A; The brackets on the right side of Equation (55) are the 
stresses or, computed from the function F, by Equations (2). The derivation 
of F, will be given subsequently herein. Denote these terms by Qo,a and 


Q,,»- In general, let 
Die eer econ b= (8) fgloina hoagie, (56) 
IO. r OT 


and let ‘subscripts denote the co-ordinates of the point at which the quantities 


are to be computed. 


Cae oN, 


and, 


: tat hare, wa tl pert Fiat de 


ory, 1 Lary Sail’ pe et in senate ach ea rwhte oon hd 


\s 22 


tuting back iho, Seton “(55), the following mis ent are abt 
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2 Ao no Mo b™ — 2 As th ma bo = oh Webs: oe (688) 
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eo ee “)In most cases it is possible to choose a = }, and, if so, further simplification } 


is Sela Placing a = b in Equations (58) and solving: 


A,= 2 Fa — Qa + oo — Qyro tt (59a) : 
er eas Sh ays 


Ss fa — Qou — os + Qro ! (590) 


The bias function, F = F, + Fy = Fo + F, + Py is now completely known 


/ and the stresses can be found by the definitions in Equations (2); thus: 


(o;)oco = 28 = Qoa + 0 — Or (z)" 
2 \a 


Hs oa — Qosa 35 = = Qy.v (2)" + 2 Ree ONS StS wo ee 5 ee (60a) 


and, 


(0;)ocy = 24 Qoa + o> — Qvo ( zn 
a 


2 
7 fa — Gove — oF Ore (£)™ Sh Dg esa) Se (60b) 


The function, F., will now be given for a dam. Let Fs = Fp + Fy, in 
which Fp is the part of F, due to the external forces acting on the boundaries, 
@ = 0 and 6 = y (Fig. 4); and, Fy is the part of F'; entirely due to the body 
forces, g, per unit volume. 

Case 1.—In this case there are no external forces acting on the boundaries, 


6 = 0 and 6 = y, on the stretch, a to b. Hence, Fp = 0, and by Equa- 
tion (56), Qr,9 = 0. ia 
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7S — [8cos@sin2y + 6sindsin*y — 4sin 6] — 
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conditions of this case. Then by Equation (56): nite 


imbrrad ob 


7 
a o- 


sin? 


ie Case 3.—Only mass forces are considered, and the function, My, must be 
_ such that oe and t,y,g are zero along the straight boundaries, 6 = 0, an 
@ = y. Substituting polar co-ordinates into Equation (21) and placing 

ge =geosPandg=gsin£g: - 4) Se, 


<< 


4 : ; eer 

q Py = | SEE sins 9 4 SERB cso cit 9 — ETE aint 0 

a 6 tan y tan y tan? y 

a + sin B sin’ # + cos 8 cos’ 6 | prea s GAbs ‘ah a4 agen 

By Equation (56): | 
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- sin Y 

2 and, 

od Qyr = py i See hae Sere ty AoE) 

in sin ° 

A It will be seen in Fig. 5 that no is always less than unity for angles between 


180 and 360 degrees. This means that the stresses found by Equations (60) 
- _ will approach infinity as the radius vector, 7, approaches zero. In other words, 
a mathematically sharp re-entrant corner cannot occur in a structure without 
- having the stresses exceed the elastic limit, contrary to the assumptions of 
} the mathematical theory of elasticity. The material becomes plastic and the, 
x stresses in the immediate neighborhood must be redistributed. 
Theory Pertaining to Fillets—In order to predict the elastic stresses it is 
necessary to abandon the sharp corner and assume a finite curvature. In this 
case the preceding formulas no longer are valid in the immediate neighbor- 
hood of the corner. Photo-elastic experimentation has shown, as might be 
expected, that a slight curvature only will effect the stresses in the immediate 
neighborhood (that is, at distances of the order of magnitude of the radius 
of curvature). The stresses, or, at the tangent points, d and e, Fig. 4, may 
‘still be computed (on the side of safety) as if the corner was sharp: 
With the foregoing assumption an approximate method is developed herein 

for the determination of the stresses on the boundary of the fillet. The origin 


xerted on as tirved: Sie piss R. If aes is no pressure, p = 0. 
‘ow, consider the stress function: 2 ' 


= (B, cos ¢ + B; sin 4) (2 + 2eRee) Fp 
p Bh. oe 2 
+ = oft p sin (@ — B) ..--. ot. ee (65) 


ols. With Equations (2) applied to Equation (65) it is easily verified a 
when p = R, op = — p; and 7,,¢ = 0 for any value of ¢. 


The arbitrary constants, B, and B., are now to be determined such that 
when p = R, and ¢ a : 


o¢ = Tig sin nie hate ek oe oe ee COR . 
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Fs ¥ ; By applying BauRodas (2) to Equation (65): . 
B. o¢ = (B, cos ¢ + B, sin 6) (24 =) — p — pg cos (¢ — 8). .(67) | 
oa Hence, by the conditions expressed in Equations (66): . 
7 . By = [ce + p + (ca + p) cos y + gR cos B sin y] R .-(68a) 9 
A 4sin y : 
: and, 

. oe 

; ; B, = [ca + p + gR sin 8] Fe ee FORIOR (68b) 


Substituting back into Equation (67) with p = R: . 
(o¢)exr = [oe + P + (ca + P) cos yj S2k + (og + p) sind — p ..(69) ~ 
ey 
The method of procedure is best shown by a few typical examples. 


EXxaMpLes CompuTEep By APPROXIMATE THEORY 


The underlying assumptions are based, to a large extent, on results of 
photo-elastic experiments conducted in connection with +he Morris Dam (see 
Part IT). It is natural, therefore, to base the example on computations for 
this dam, because verification is possible. A cross-section of the Morris Dam 
is shown in Fig. 2. Theory and experiments have shown that the stresses in 
the upper two-thirds of the dam can be computed accurately by Equations (17 ) 


er with talk fects: ‘presiure/ pl — ae , on Pine up-stream ‘teed qd. 
- constant pressure, p, on the foundation above the dam. The body forces will 
ig considered subsequently; that is, g = 0. 

For this loading the value of F, = Fp is given in Equation (61). rome 
; Ba ond b (Figs. 2 and 4) are chosen as r = a = b = 0.3; the stress, ga, can 
be computed with reasonable accuracy by superposing Hqhetions (35) and 
(89a); q: in Equation (35) is the total horizontal water pressure on the dam ; 
uniformly distributed over the base 4; and po in Equation (39a) is ae 
pressure, p, on the up-stream foundation. 


P 


Z oe In the present application, a is an up-stream point; q@ = oh : and, — at 
eo pai ey 0) p 1.16 ea 
1, = 0.86. The total stress, then, is og = + log=—— p= 0: 500 p. 
0.867 0.3 


- Point } is sufficiently far from the base so that op may be computed as oz by 
4 Equation (18a). In the present application, K = tan 41°52’ = 0.8963; 
ec — 0.7; y = 0; and os = + 0.871 p. ; 
The up-stream corner angle, y, = 267° 08’. By Fig. 5, no = 0.55; and 


. nm, = 0.93; and, by Equation (46), mo = 0.895; and m, = — 0.818. The quan- 
e tities, Qo,a and Qy,», are now found by placing a = b = 0.8 in Equa- 
, tions (62); thus: Qo,r = — (1 + 0.057) p; Qur = — (1 + 0.00257) p; 
Qo a = — 1.015 p; and Q,,» = —1.001p. By substitution in Equations (60) ; 
4 : 1\"* 1\o" om 
‘ (7) e=0 = [0.0045 (2) — 0.6078 (2) — 0.05 r — i|p ...(70a) fs 
Tr r 
sand, e 
: 1 0.45 il 0.07 ie 
(or)exy = [0.6045 () + 0.6078 (4) — 0.0025 r — |p . .(70b) 
1 r | 


In Equations (70), it is to be noted that: r is measured in units of height 
of dam and p is the actual hydrostatic pressure at the base. By Equa- 
tions (70), or becomes infinite as r > 0. To eliminate this condition assume 
a finite radius of curvature, R = 0.025, at the corner with tangent points at 
d= e — R cot 0.57 = 0.0238 (see Figs. 2 and 4). 

Now, compute ca by Equation (70a) and oe by Equation (70b) for 
r=e=d = 0.0238; thus, og = + 1.943 p; and oe = + 3.523 p. Substitute 


in Ep ahatien (69) with g = 0: 
(og)oan = [— 4.8818 cos @ + 2.943 sin § — 1] p.....-.-(71) 


The maximum occurs when ¢ = 146° 07’; and og (max.) = + 4. 278 p. The 

stresses computed by Equations (70) and (71) are plotted in Fig. 2(a). 
Up-Stream Oorner, Mass Forces. —Now, consider mass forces only; that 

is, let p = 0. Only vertical forces are assumed to exist; hence, B = 90°; 


G2 = 0; and g = gy = weight per unit volume. The stress, oa, computed by 


‘compute the quantities, Q, ete., by (6: 

37° 08’; and a = b = 0.3; thus: ue = s 0.050 ngs Qrie= th 
ooa = + 0.015 g; and ee efi: BOOB gh Woy OF) mh senriieda . 
i ‘By substitution in Equations (60) with p = ys and no and NM as before: gi 3 


(oes = ie 0. aor (3 Ly" ers so5(2)" 4 +0. ae rd ; 1 (P20) 


i 
- 


: ee by 1\°8 4 \0% Rae fe 

ets aimee ~0.0is(2)" 4 + 1.001 r a0 . . .(72b) 

in which r is measured in units of H. Again, the stresses become infinite 

for the sharp corner when r > 0. The radius of curvature, R = 0.025, is 
assumed as before, and d =e = 0.0238. By Equations (72) with r= a= | 
e = 0.0238: og = — 0.9596 9g; and oe = — 2.0657 g. Substitute in Equa- © 
tion (69) with p = 0, thus: 


(o;)oxr = [+ 2.0202 cos $ — 0.9695 sin $] g...---.---.- (78) 


_ The maximum occurs when ¢ = 154° 36’; and maximum g¢ = — 2.287 g. 
3 The stresses computed by Equations (72) and (78) are plotted in Fig. 2(b). 
is u For the section under consideration let p = 108.5 lb per sq in., and g/= 2.5 .p 
ae: = 271 lb per sq in., corresponding to 156-lb concrete. 

aes _ In computing ie stresses for the down-stream corner the procedure is 
rps exactly the same as for the up-stream corner. The origin is chosen at the 
ae toe with the Y-axis along the down-stream foundation and the Y-axis posi- 

; tive into the foundation. If no pressures act on the down-stream face and — 
foundation, Fp = 0; and Fy, is again given by Equation (63), with y equal to the 
eS re-entrant corner angle of the toe, and 8 = 90°, if the foundation is horizontal © 
ee and if inertia forces do not exist: Point b (Fig. 2 (6)) is chosen on the 
ae! down-stream face and Point a on the down-stream foundation, such that 


i= = — als . The radial stress, op = oz + oy, is computed by Equa- 


\ sin y 

tions (18) in the case of hydrostatic loading and by Equations (22) in the 
case of mass forces. The stress, oa, is computed by Equation (35) and becomes 
negative. The results are plotted in Fig. 2 for the Morris Dam, assuming a ~ 
fillet radius of 0.025 H as at the up-stream corner. 

The Morris Dam is not designed with fillets, but it was necessary to 
manufacture the models used for the photo-elastic experiments (scale, 1: 600) 
with a 4-in. radius (0.025 x height) at the corners to prevent failure when 
hydrostatic load was applied in the absence of body forces. For this reason 
the same radius was used in the computations for stresses in the fillets. 

The approximate method of analysis is recommended only for preliminary 
work, in order to fix suitable dimensions in the base region of the section. In 
the final analysis one of the methods explained in Application III should 
be used. 


ee ee ee re ee 


d- at, for stable gravity a i 04H to 0 


itable at the hee 0 on the heig 
the dam and the quality of the masonry. The up-stream fillet is deter- 
ned for the case of an empty reservoir, the down-stream one for the case 
fa full reservoir. genet’ aie Mite ee 


“PART IL—PHOTO-ELASTIC EXPERIMENTS IN CONNECTION — 
‘WITH THE MORRIS DAM, IN CALIFORNIA 


Brier OuTLINE oF THE PuHoto-Exastic PHENOMENON seen 4 
; Although the principles involved in photo-elastic stress analysis may be 
_ found in treatises on photo-elasticity® a short description is given herewith, 
_ for the sake of completeness. ee 
Tz About 125 years ago the English physicist, Sir David Brewster, discovered 
that isotropic media become doubly refracting when placed under stress. 
_ They then behave optically as crystalline media. The planes of principal 
stress correspond to the crystal planes. If a plane-polarized monochromatic 
ray (see Fig. 8), is sent normally through a plate stressed in its own plane, 
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Through the Model 


é ae RIG. 8: 


Nicol’s Prism, 
Polarizer; Crossed 


fae ks 


the ray will be resolved into two component rays each vibrating in a plane of 
principal stress. The two ray components will usually have different veloci- : 
ties in the plates; therefore, while passing through the plate, the wave front i 
of one ray will get ahead of the wave front of the other ray a small dis- 

tance, 8, which is called the relative optical displacement of the two ray 
components. If, after leaving the plate, the two components are combined so 

as to vibrate again in the same plane, interference will take place due to the 

relative displacement, 8. Complete interference will occur when 6 = 0, 

1, .... n, wave lengths. This furnishes a means of measuring the displace- 

ments. It has been shown theoretically by Clerk Maxwell", and verified 
experimentally that the difference in magnitudes of the two principal stresses, 

op and go; at the line of passage is in direct proportion to the relative dis- 
placement. The latter can be measured in terms of wave lengths of the 

light by counting the order, n, of interference, hence: 


also Say 


Pp IT oma 1 YP 


; 10 “photo-Hlasticity”, by Coker and Filon. 
“1“On the Equilibrium of Blastic Solids”. by Clerk Maxwell, Fransactions, Royal Soc. 
Edinburgh, Vol. XX, Pt. I; also, “Photo-Hlasticity’, by Coker and Filon, p. 198. 
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in which C is a constant which depends on the type of material, the thick- 
ness of the plate, and the wave length of the light. It will be shown later 
how C is determined. (It is to be noted in passing that the number of fringes 
is directly proportional to the thickness of the plate.) 


The principal stress difference, cg» — oe, is equal to twice the maximum 
shear at the line of passage; hence, this quantity is known when the fringe 
order, n, is observed. Furthermore, op — go will change continuously from 


point to point in the plate; therefore, the optical displacements, 6, will do 
likewise, and equal integral values of n will occur along continuous interfer- 
ence bands ealled “isochromatics” which then may be defined as lines of 
equal maximum shear. 
The plane polariscope is shown schematically in Fig. 8, with explanatory 
‘notes. The two Nicol’s prisms, the polarizer and the analyzer, are “crossed”; 
that is, their polarizing planes are set at right angles so that no light will 
reach the screen, when the plate is not in the field or when it is unstressed. 
As already explained alternate dark and light bands, “isochromatics”, will 
appear on the screen when the stressed plate is placed in the path of the light. 
Along the edges of the model one of the stresses, gp Or oo, is known, and 
the fringe order, n, can be observed so that the other stress (og or op) may be 
computed by Equation (74) when C is known, 


CALIBRATION 


The simplest way to determine the value of C is to subject-a beam of 
rectangular cross-section to uniform bending. The beam must be cut from 
the model material. The fringe pattern obtained in a bakelite test beam is 
shown in Fig. 9. The dark band at mid-height is the isochromatic of zero 


Fig. 9. 


ISOCHROMATICS IN CALIBRATION BEAM. 


order (n = 0), or the neutral axis. The constant bending moment and the 
dimensions of the beam section can be measured so that the horizontal stress 
at the extreme fibers can be computed. In the present experiment, gp = 1055 
Ib per sq in. No load occurs at mid-span; that is, cg = 0. The fringe order, 
n, is counted at the top or bottom (in this case, 6), the neutral axis being 
of zero order. Hence, by Equation (74), 


oO — 
C= pe a's 5-0 3'n artic Re int aS 
nN 


1 055 fs 3 
= 176 lb per sq in. The even spacing of the parallel fringes shows 


or, 


eh ie eee 2 
the engineering assumption ar stress distribution hol 
for a beam of uniform thicknes 
; thickness of the beam was 0.346 in. and white light with a blue Wratten filte 
- (about 4100 Angstrom) was used. This filter and models, 0.346 in. thick, 
__ were used in both the experiments to be described. If the thickness, gee Ate 
: _ the model were different from that of the test beam, 7’, merely multiply C in ag 


- Equation (75) by ©. | ; 
e : 3 | 


-- It may happen that Ray V, (Fig. 8) is vibrating in the plane of the prin- ee 
3 cipal stress, op, at the line of passage. In this case Ray Va (= Ray Vpeowill = 
pass unmolested through the plate and the orthogonal component, Ve, will be a: 
_ zero. Therefore, nothing will happen on the screen and darkness will prevail — “g 
ee. in the image of the point of passage, the Nicol’s prisms being “crossed”. ta 
When the polarizing plane of the polarizer is set at a known angle the rays a 
automatically will seek the points in the model at which one of the principal % 
a stresses is parallel to the polarizing plane. The stress directions in a model 
fe will generally vary continuously (except at singular points) so that the dark a : 
o points just described will also form a continuous band in the image known 
e 
4 


\ 


as an “isoclinic” or locus of points having their principal stresses parallel re | 
to the known plane of polarization. If the polarizer and analyzer are nowt 
rotated together through a certain angle, another isoclinic will appear, etc. 
In this manner the direction of the principal stresses can be found through- EY 
out the model. : 
Although the location of the isoclinic changes with each prism setting, the : : 
system of isochromatics will remain stationary and, therefore, are easily a 
distinguished. 
By using circularly polarized light, instead of plane polarized light the 
isoclinics can be eliminated. In this experiment the light vector, V,, is made 
to rotate with a high, uniform, angular velocity while passing through the 
model and thus it can have no directional preferences; therefore, isoclinics 
‘cannot oecur. Circularly polarized light is produced by inserting what is 
known as a “quarter wave plate” between each Nicol prism and the model. 
These crystalline plates are usually obtained by splitting mica to a thickness 
corresponding to a relative retardation of the two component rays, equal to 
one-fourth the wave length used in the experiment. The theory of this 
phenomenon may be found in treatises on passage of light through crystalline 
media.” 
By using materials of low optical sensitivity (that is, a very large C in 
Equation (74)) such as plate glass, the isoclinic lines which are only a 
function of the stress direction can be made to appear sharply long before 
the model is stressed sufficiently. for the jsochromatics to show except as a 
slight illumination of the glass. The latter effect makes the isoclinic stand 
out clearly as a black line on a light background (see Fig. 8). Glass without 
initial stresses must be used, and is distinguished by holding the unstressed 
glass plate in the field of the polariscope before cutting it into the desired 
model. ‘A 
2 “photo-Blasticity”, by Coker and Filon, p. 73. 
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Notes: 
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¥ lines of principal stress can then be integrated by use of the isochromatic lines, 

i the known boundary forces, and Equation (74). Methods for this work are 
given in treatises on photo-elasticity.® 

In many types of problems, even when body forces are to be dealt with, the 

following method is convenient for determining the magnitudes, op and go, 

along given lines of principal stress. The method is essentially structural 

and was developed by considering that the structure is composed of two 

systems of orthogonal arches bounded by the lines of principal stress. These 


18 “Photo-Hlasticity,” by Coker and Filon, pp. 2, 47. 


certain reactions. 5 4 
Consider the element, A,, B,, As, Bs, bounded by the principal lines 


#nr 


_ stress shown in Fig. 11(a). The only forces acting on the element 
Puy P,, 2, Qu1, and Qs, 1, normal to Surfaces A, B, ete., and perhaps the mass 
_ force, M fet TE Pant, O12 and M,,, are assumed to be known in magnitude ‘ 
and direction, then P:,, and Qs can be found by constructing the force 

polygon. ; poe 
Proceed to the element, A. B. As Bs, etec., until all the forces, P:,n and 

5 Q2, n, are found. It is necessary to start the construction at the boundary, — a4 
_ Ax, An, As, ete., where the forces, @1,n and P,, are known or can be found by | 
j the isochromatics and boundary forces. When all the forces, Q2,n, are found, — 
consider B,, B,, B;, etc., as a new boundary, apply the now known forces, Qs, n, 

reversed, and the known forces, Ps,., and find the forces, Qan and Pan, by 
q construction, etc. The average stresses over the elementary sides are found _ q 
_ by dividing the forces by the corresponding sides. Finally, the principal = 
stresses at B., say, will be: . 


7 == 
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Tt will be seen that Equations (74) furnish a check on op — og at all points, 
which should be made after the completion of the stresses along each line 
of principal stress. 3 
The construction is easily arranged in a force diagram so that duplications 
of lines are avoided. Note that if the lines curve rapidly small elements must 
be used. 
, The integration of stresses along principa 
. often not very accurate in complicated stress patterns. 
have been made to obtain other methods. 
A two-dimensional state of stress ig determined completely when the prin- 
i cipal stresses, op and go, and their directions are known at all points; op — a 
. is known throughout by the isochromatics and by Equation (74). The direc- 
tions are known by the isoclinies; hence, if op + oe can be found, the state 
of stress is determinate; og + op can be obtained by measuring the lateral dila- 
tions of the model with an extensometer,® or by the membrane analogy.” 


1 lines of stress is tedious and ark. 
Consequently, efforts 
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-__ photo-elastic experimentation. 
“i proceedings, Am. Soc. C. E.,.May, 1935, p. 597. 
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The weight of the model itself is too slight to have any effect unless it is 
magnified many times. The easiest method would be to rotate the model 
and observe it whenever it passes through the optical field. This method has 
not yet been attempted, and it will require a steady rotor mechanism, free of 
vibration, to obtain a good- photographic exposure of the photo-elastic lines. 
The polariscope used by the writer would be admirably suited for such an 
experiment. The model could be attached to one end of a rigid rotor and the 
photographic plate to the other end, each passing through the optical field 
simultaneously with the same relative angular velocity. A new method has been 
developed by the writer in connection with experimentation on the Grand 
Coulee Dam, by which only the corrections to the straight-line stress distribu- 
tion are involved. An interesting equivalent boundary loading has recently 
been pointed out by M. A. Biot.” ; 

Stresses have been determined experimentally for the section shown in 
Fig. 2, under all critical loading conditions, including earthquake accelera- 
tions and uniform shrinkage relative to the foundation. 

Hydrostatic Load on the Up-Stream Face and Foundation: Isochromatics. 
—A bakelite model, with blue filter and quarter-wave plates, was used to 
produce isochromatics. The model scale was 1: 600; the photographic scale, 
1: 900; and the loading factor, 4.6, which means that the hydrostatic pressure 
at any point in the model was 4.6 times as great as at the corresponding point 
of the actual dam. The up-stream face was loaded linearly by three levers, 
each applying the correct force through knife-edges to a steel shoe. The uni- 
form distribution on the up-stream foundation was applied by a single lever 


Transactions, A.S.M.E., Vol. 2, No. 2, June, 1935. 
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tion (74) ;. thus: tmx. = ———* = — = 
aR Sy ree 2. ieee : 
_ by Equation (75). The shear stresses at corresponding poiuts of the actual 
dam will then be, in pounds per square inch, BV) < 


a : 
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‘@ 


88 n 88 n | " 
a aes LN teria «sa eK 
Load factor 4.6 m ‘pene ce) 


in which n is the fringe order given in Fig. 12(b). The actual boundary ~ 
_ stresses are also determined by Equation (74): . 


: 2 SE ee ec LE ee aera a ef ee 78 

_ Load factor ; ibe : 

ig On the rear face and foundation, P = 0; hence: 

4 Q = — BBM. 0 re eee ere tee ees (79) 

- On the up-stream face and foundation of the dam, the hydrostatic pressure is, 

Phe. abs ve 
+ in pounds per square inch, Q = = 62.5 h’ = — 0.434’, in which h’ is the . " 


_ 144 : = a 


28 Ras 


_ depth, in feet, to the point in the protetype, = hh xX 900 = 75h, and h is 
12 


the depth, in inches, to the corresponding point in the photograph. Hence, 
by Equation (78), ’ 


Fah S= 38a — BO DN soe Sete sine Foe oC 80) 


in which n is the fringe order at the point in question. | 
The experimental results obtained by substituting the fringe order, n, x 

along the boundaries of Fig. 12(b) into Equations (79) and (80) are plotted 

in Fig. 2(a) with the theoretical values obtained in Part I of this paper. 
Tt will be noticed that in the upper two-thirds of the dam the stresses 

obtained by Equations (18) agree with experimentation. In the base region 

the boundary stresses are in close agreement with the approximate methods 

discussed under the heading, “Application IV”, in Part T. It can be shown 

that the isochromatics for the infinite wedge are concentric ellipses with their 

centers at the vertex, agreeing with the upper isochromatics in Fig. 12(b). 


A slight disagreement between the computed and experimental values may be 

expected due to the finite crown width, which is neglected in Equation (17). 
Hydrostatic Loading: Isoclinics—A plate-glass model and white light 

without the quarter-wave plates were used to produce isoclinics. The photo- 
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ao graphic scale again was 1:900. It is to be noted that the isoclinics are inde- 
yf pendent of the loading factor as defined previously, except that these lines 
: become sharper with increased factor. In order to cover the entire region it 
ty is only necessary to obtain the isoclinics for various prism settings, a, through 
: a range of 90°, as shown in Fig. 10(b). A separate photograph is necessary 
| of course for each prism setting, a sample of which is shown in Fig. 13 cor- 


Fie. 138—IsocLinic LINnS FOR PRISM SETTING, 
a, IN Fie. 10(b). 


responding to the prism setting, as, in Fig. 10(b). The isoclinics are plotted 
in Fig. 10(a), and the lines of principal stress are obtained in the following 
manner: Through a, of Fig. 10(a)draw a, a, perpendicular to a, of Fig. 10(b); 
then draw a,a; perpendicular to a, ete, The points, aza3, etc., are taken 
approximately midway between the corresponding isoclinies, 1, 2, ete. 

It can be shown that the isoclinies for the wedge are straight lines radiat- 
ing from the vertex, agreeing approximately with the upper two-thirds of 
Fig. 10(a). The isoclinics perhaps show more clearly than any other argu- 
ment the effect of the foundation on the stresses in the base region.. The 
corners in the models were rounded; with sharp corners, all isoclinics would 
be represented at each corner. 


eg ee - the base and in the foundation were deter- 
: The forces representing the weight of the dam above the’ base we 
ores through a single lever distributing the weight to four points ‘in su 
a manner that the resultant was correct at the base. The isochromatics a 
D otted in Fig. 12(a) near the base and in the foundation. The photographic — 
scale was 1:900, and the loading factor was 1.75, based ‘on 156-Ib conerete. 
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: The experimental results near the base obtained by substituting the fringe 
order, n, at the boundaries of Fig. 12(a) into Equation (82) are plotted in- 
Fig. 2(b) together with the theoretical values. It will be seen that the 
stresses obtained in Part I agree fairly well with the experimental results. 
It is of interest to note that the compression at the down-stream corner is of 
considerable magnitude. 

The results in Fig. 2 can be applied to a dam of similar proportions, but of 
height, H’, merely by changing the stress scale in the ratio, 250: H’. If the 
concrete in the new dam weighs w lb per cu ft, the stress scale in Fig. 2(b) : 

is changed further in the ratio, 156: w. The stress concentration in the base z 
region at the down-stream face, due to weight, extends only a short distance a 
into the dam. 

The foundation was clamped between rubber gaskets to the steel frame ~ 
along three edges. This arrangement allowed elastic deformation to take 
place so as to produce, as closely as possible, the effect of an infinite half- 
plane. © The photo-elastic results in the foundation, of course, are only valid 
in the region near the base. ‘ 

In the interpretation of experimental results from model to prototype it is 
important to note that in two-dimensional problems: (1) The stress distribu- 
tion is independent of the elastic properties; and (2) the stresses in the pro- 
totype are directly proportional to the stresses at corresponding points in the 


model. 

The discrepancy between t 
Fig. 2(a) may be due partly t 
fillet of the model and partly to a s 


bs eal ed eas. Ration) Reh 


uu. 
OA 


ere 


he experimental and theoretical curves of 
o incorrectness in the pressure, Q, at the 
light shop error in the radius of the fillet. 
Other experiments with larger fillet radii showed much better agreement 
between experimental and theoretical values of fillet stresses. Very small 
radii were used in these experiments in order to produce as closely as possible 


the effect of sharp corners. 


eyond es scope me this paper to eae on ee subject. A ae pla 
stem of electrical strain-gauges was installed in a section of Morris Dam 
‘ich, in time, will throw further light on this point, on the question of | 
rinkage effects, and on lateral stresses in 1 long gravity dams. ie. 


5 LOOD AND EROSION CONTROL PROBLEMS AND 


THEIR SOLUTION 


By E. COURTLAND EATON, M. Am. Soc. C. E. 


Synopsis 


The flood and erosion control problem treated in this paper occur in highly 
developed areas where there is a relatively low seasonal run-off; these areas, 
however, are subject to brief, although violent, torrential storms, that result. 


in floods of exceptionally high intensity. When fires denude the sparse 


vegetation from steep mountain slopes the flood intensities increase and © 


added débris hazards occur, due to erosion. 

Basic precipitation and run-off records are given in this paper, as are 
engineering methods of constructing hydrographs of expected Acéas and 
the regulation needed for control. Similarly, measurements of erosion quan- 
tities are presented with suggested solutions of control and a method of 
avoiding unnecessary capital expenditures in advance of requirements. The 
incidental conservation of flood waters for domestic use and for irrigation 


is discussed. 


INTRODUCTION 


Los Angeles County in California, with a population of about 2 250 000 
(of whom more than 85% are residents of cities), covers an area of 4115 
sq miles. Its population increase has been nearly 180% since the last major 
flood (1914), and assessed property values rose, in the same period (1914 to 
1931), from $850 000000 to more than $4000 000000. The County’s petro- 
leum, agricultural, motion picture, and manufacturing industries have an 


aggregate annual value of more than $1000000000. It is estimated that 


only 8% of the present population have experienced, or have a realization of, 


the 1914 flood. 


Property valued at 
is subject to possible inundation due to floods.’ 


more than $300 000000 and containing 3880 000 per- 
Flazards may be classi- 


er, 1935, Proceedings. 


sons, 


Norr.—Discussion on this paper will be closed in Decemb 


1 Cons. Engr., Los Angeles, Calif. 
2“Comprehensive Plan for Flood Control and Conservation,” by B. C. Eaton, 1931. 


ndation by ei er the overflow of present 
_ of the rivers from existing to new locations; aaa in the foot-hill ar 4 
vila) ‘are from sudden flash floods, and débris flows, particularly following 
- -water-shed denudation by fires. Valley inundation may be anticipated, by at : 
east hours from distant mountain rainfall records. Foot-hill débris flows | 
a may result if burned catchment areas, after becoming saturated with water, 
‘i are subjected to high-intensity storms. Such flow occurs without warning — 
and may follow storms of high intensity or cloudbursts, almost immediately. 
ae The factors that combine to menace life and property in this region are: 
eA (a) Exceptionally precipitous gradients; (6) short mountain streams; (c) the 
unstable character of deep weathered mountain cover held from erosion by 
Sit sparse vegetation; and (d) the characteristic violent torrential storms. Pro- 
tection from floods may operate to conserve the average annual 120000 
acre-ft of flood water that is now wasted. 
The paper describes physical conditions and methods of solving hota . 
especially those dealing with the relatively new engineering field of erosion ~ 
control. ‘Necessarily, full “details cannot be included, but the foot-note 
references serve as a bibliography of public reports. 


ey PuysicaL ConpITions 
< a Of the 4115 sq miles that constitute Los Angeles County, 70%, or 2758 
# ; sq miles, is included in the south slopes which drain to the Pacific Ocean f 


(see Fig. 1). The area comprises 1589 sq miles of mountain water-shed 
f (809 sq miles of which is in Federal forests); 285 sq miles of small water- 
eh sheds; and 884 sq miles of valley plains sloping toward the ocean. About — 

3 210 000 acres of valley land is under irrigation, and 235000 acres is devoted 
_to industrial and domestic purposes. 

The high mountain areas form a northern boundary of approximately 45 
miles, crest length, from which the two major collecting rivers—the San 
Gabriel and Los Angeles—originate and, after traversing the valleys, they 
converge and discharge into the ocean at points only six miles apart. Moun- 
tain peaks rise to Elevation 10 080 (United States Geological Survey datum). 
The main rivers are forced to converge by the foot-hills that project on to 
the plains—the Santa Monica hills from the west and the Puente and San 

) José hills from the east. 

The San Gabriel River has a tributary drainage of about 800 sq miles 
and the Los Angeles River one of about 900 sq miles. The largest catchment 
of the San Gabriel River is 213 sq miles of the San Gabriel water-shed; 
that of the Los Angeles River is the Big and Little Tujunga water- 
sheds with 137 sq miles. The remaining drainages of both rivers are made 
up of numerous small streams originating in more than 100 mountain and 
foot-hill water-sheds ranging in size from 0.5 sq mile to 30 Sq miles. 

A third drainage system is Ballona Creek fed from 134 Sq miles, its 
principal supply being the run-off from the streets of Los Angeles, Holly- 
wood, and Beverly Hills, transported to Ballona Creek through storm drains. 
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is near San Pedro, Calif. te: a a 
Jong the coast from Nenies Bie ste: Ls ‘he ee nee Yount 
Pes bout twenty individual streams flow directly into the Pacific Ocean. 


rgest are Malibu Creek, with a drainage area of 67 sq miles; Topanga | 


(20 sq miles), and Mandeville, Rustic, and Sullivan Canyons (12 sq miles), 
the latter group passing through Santa Monica, Calif. 


Mountain and foot-hill drainages have well-defined canyon channels “on 


dence of active channels, and, during the dry cycle (1914 to 1931), home sites _ 


; were sold readily in foot-hill territories by real estate agents who were. them- 


selves possibly ‘ignorant of hazards. Normally, mountain and foot-hill water- 


sheds are covered with a chaparral type of vegetation, efficiently functioning 


as retarders of flood flows and controlling erosion. 


Historic Fioops 


The most recent major floods occurred in 1914 and 1916. In 1914 there 
were two storm periods: One on January 25, with a maximum intensity 
of 2.6 in. in 24 hr; and the other on February 18, with an intensity of 4.26 in. 
Exclusive of harbor damages the estimated loss from these two floods was 
$10 000 000. Conservatively estimated representative flood peaks were® as 
shown in Table 1. Many lives were lost, thousands of people were made 


TABLE 1.—ReEprESENTATIVE PEAK FLows 


Peak run-off, in 


Stream Drainage area, cubic feet per 
in square miles second per square 
mile 
SamdGabriell River.)s.cc.cs a oemiae aaa sk «tk Maen ete eck 229 117 
Big Tujunga 118 115 
Arroyo Seco 39 366 
EA WOLOL tbh tela ain lee, Co oh eickate steieRre toe oteeietaa nei serene Seek bi 550 


homeless, thirty-five bridges were destroyed, and, for six days, there was little 
communication with the outside world. 

Engineers concede that this was by no means a record flood. During pre- 
vious floods within a 70-yr period, main rivers have materially changed their 
courses, the Los Angeles River changing from westerly to southwesterly, 
and the San Gabriel River cutting its new channel from 3 to 6 miles south- 
easterly. A new river, the Rio Hondo, cut its way and connected the San 
Gabriel and the Los Angeles River. 

Early settlers testify to the occurrence of severe floods in 1815, and again 


in 1825, and, successively, five major flood periods are recorded in 1861-62, 
1884, 1886, 1889, and 1914.? 


* Rept. of Board of Engrs. on Flood Control to the Los Angeles County B 
Bee eet January 25, 1917, by H. Hawgood, Ch A aninecce se 
Ptumiien Monten? A dee £4 Charles T. Leeds, J. B. Lippincott, and 


as ‘i precipitous grades of 20 to 30%, sharply flattening from 8 to 10% where they e 
-debouch on to unstable delta areas. In many instances cones show little evi- 


SS 


a 


yment of 9 

5 55 sq miles has increased the rates and quantities of run-o i 
_ decreased the natural absorption and consequent replenishment of subter 
_ ranean storages. ; eee 
3 _ A curve prepared in 1931 (see Fig. 2) shows past trends of seasonal rain 
- fall. After this curve was plotted a prediction was made of an uptrend in 

a wet cycle with increasing wet periods, comparable to the period, 1881-90, 
> in which interval one or more major floods would occur.’ | bit 
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Fig. 2—-GuNERAL TREND, SHASONAL RAINFALL. 


RAINFALL—RUN-OFF 


A co-ordinated system of 237 rainfall and 197 run-off stations was installed 
by the Los Angeles County Flood Control District commencing in 1927. The. 
more important stations were equipped with recorders. The year, 1931-32, 
not a year of general flood, had valley rainfall totaling 114% of the normal, 
that in the mountains totaling 113% of its normal. The seasonal run-off 
from the San Gabriel River exceeded the combined run-off of the three pre- 
vious years. Summarized records from three drainages follow: (a) San 
Gabriel-West Fork—area, 49 sq miles; (b) Little Tujunga Creek—area, 21 
sq miles; and, (c) Ballona Creek—area, 112 sq miles. , 

The San Gabriel and Tujunga water-sheds are mountain catchments with 
cover in good condition. The San Gabriel water-shed was burned in 1924, 
but its brush growth is now (1935) restored. The Tujunga water-shed has 
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ps January 3 | February 7 


aH 


Days of << itation, 
precipitation a fechas . 


(a) San GABRIEL WATER-SHED 


October 1 | December 24 18 
December 25 | January 2 


7 
Total for first storm period... . 25 
7 
8 
5 


February 8 | February 28 
March 1 | May 1 


20 


Practically none 


October 1 


December oa 


_ December 25 | January 0.12 tor 
Total for first storm period. . 0.12 0.6 
Febrasy 8| February 28 31.3 
Total for second storm period. 15.0 
ona Alar cages ag 3.40 200 0.88 28.0 
Total for first storm period... 15 cat tas i tee | paar ak ieee 
Poway (81 Heoniaee ae a Pe 1:95 1:33 27.0 
Total for second storm period. reat GAAS Re Cae aeereray ao ae > paed.c ue 


tion, maximum 24-hr intensities, and the resulting run-off. The distribution 
of the San Gabriel run-off was affected by snow. There was 6 to 8 in. of snow 
above Elevation 4000 during December, 1931; and 4 to 8 in., during January, 
1932; the snow disappeared by February 26, 1932. 

Both mountain catchments (particularly the Little Tujunga, the flow 
of which was not affected by snowfall), show the relatively large proportion of 
annual precipitation required to produce sufficient saturation to permit run- 
off. Mountain and foot-hill catchments are deeply covered with weathered 
material, ranging from the porous disintegrated granites to the less porous 
soils with higher clay content; consequently, separate values for any water- 
shed under consideration must be determined. The flood flows will depend 
upon: (1) Slopes, character of soils and cover (organie growth, litter, ete); 
(2) the quantity and rate of rainfall required to produce saturation; (8) the. 


; e gs: li mn 3 on Ohi eases MEAD OTE Bae: 
- On the Little Tujunga water-shed from February 6 to 10, 1932, the 
: tal 5.13 in. with 24-hr intensities, on three consecutive days of 2.10, 
~ and (1.15 in. The run-off in the succeeding five days totaled 1.05 in., | 
ratio of run-off to rainfall, of —20 per cent. ae 


Be quick response of run-off to rainfall on Ballona Creek. A total rain of 0.17 in. ; 
was followed (November 15, 1931) by a 24-hr rain of 0.96 in., of which 17% — 

appeared in a few hours as run-off. Table 3 shows the rainfall-run-off rela- . 
tion for maximum month, day, and hour, for the three catchments. Gs 


eo TABLE 3.—Comparison Between Ratratt AND Run-Orr, 1931-32 


ai, - Maximum Monts Maximum Day Maximum Hour 
Water-shed | Rainfall, | Run-off, | Per- | Rainfall,)Run-off,| Per- | Rainfall,|Run-off,| Per- ae 

A _ ip eit centage | | in yan centage in in centage 

on inches inches |ofrun-off| inches | inches |of run-off] inches | inches |of run-off © 
nS nr 
«San Gabriel...| 19.74 5.8 29 3.80 1.58 41 0.69 | 0.08 16 

¢ Little Tujunga 7.71 1.47 19 2.10 0.48 23 | Bera hee bis 

Ballona Creek 5.95 1.36 23 2.09 0.60 28 -0.27 0.07 26 
. i 


A representative number of scattered records of 24-hr and corresponding 
1-hr intensities is given in Table 4. ; 


TABLE 4.—REeEpRESENTATIVE Twenty Four-Hovur Anp CoRRESPONDING ~ 


a One-Hour RaAtNraLu INTENSITIES - 
ee 


Elevation, in Twenty-Four Hour 


: Years Fost CU. RAINFALL One-hour 

Station of e sa nae xeintel 

in inches 

ages Survey) Date In inches 

Opids Camp.........--+-++- 17, 4 480 4— 5-26 12.30 2.00 
IML GOOVVLIBGE ice cities =< ares 30 5 850 12-19-21 11.26 1.10 
Alder Creek.........----+-- OP Oe ae ee 12-19-21 7.40 0.80 
San Gabriel Intake......... Ee 1 250 4- 5-26 6.76 1.16 
Baines fel. see cae vice os 17 2 500 4— 5-26 5.95 1.10 
Valley Forge....-..-- thre .; 3 400 2-15-27 5.58 0.87 
Sister sPldters. acts cee some to's <4: aes hehe tait ® 4— 5-26 4.77> 0.62 
iy 3 300 4— 5-26 6.44 0.76 


Coldbrook Camp 
SS SSS SS nn ee 


FLoop HyprocraPpH DETERMINATION 


Studies of the Big Tujunga water-shed are given herein, showing a method 
the results, and the regulation required. 


of computing possible flood flows, 
t individual tribu- 


This catchment, as yet only partly controlled, is the larges 
tary of the Los Angeles River drainage. Its control is a necessity because 
of the flood hazard to the thousands living adjacent to the Tujunga Wash and 
to those living along the Lower Los Angeles River. That most of these people © 


have not, personally, observed conditions causing or accompanying a major 


flood makes the hazard none the less real. 
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Note: 
Run-off Season Ends Sept. 31, 
of Year Designated 
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Indices of Rainfall and Seasonal Wetness 
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15 20 25 30 
~ Seasonal Run-off in Inches of Depth 


Fie. 3.—Big Tosunca Watmr- Surep—INDICATED SBASONAL Run-OFF, 
In INcHES OF DEPTH, 


*Rept. on Big Tujunga Creek Flood Control and Conservation Posstnilivien by 
Franklin Thomas, M. Am. Soc, C. E., September 22, 1931. 


infall records « are availaple! at thirteen stations in, or apple os the 


A correlation was prepared from the longer rainfall-run-off records on the 
San Gabriel water-shed. The results, together with actual Tujunga records, 
were plotted against the master indices of mean seasonal rainfall and are 


shown in Fig. 3 in which the seasonal run-off is expressed as depth, in inches, 
on the water-shed. 
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2 - Seasonal Run-Off, in Inches of Depth: 
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; Seasonal Rainfall, in Inches of Depth (Composite 
4 : for Tujunga Water-Shed) : 
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Fifty-eight year aeetace ee ts he a Die 
Maximum .... E aim eektte oa wees 65.8 
6.8 


_) Minimums ¢.. 2 ssf. 5% Bee Felt Biot BN afte Pas ME 
Ratio of Run-Off to Rainfall, in Percentage: 
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Rainfall indices for years of largest flows are shown on Fig. 4. In arriving 
at the expectancy of a major flood it was decided to use the flood of 1888-84 
ie as a criterion for protection, because all indications pointed to that year as 
the one of major flood in the 58-yr period. 


a. 250 


Rainfall Indices Expressed as Percentages of the Mean 
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Year 2 
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Fig. 4..-RAINFALL IDICHS FOR YRARS oF LARGE FLOODS. 


The mountain soil cover is disintegrated to great depth; practically no 
rain falls between May and October, both inclusive; and hydrographic experi- 
ence has shown that on this area a total of about 10 in. of rainfall causes 
saturation. If this is closely followed by high-intensity storms, floods are 

4 inevitable, the magnitude rising rapidly with increasing intensities. 


on Index 
Run-off | of | ° 
depth, in | season- 
| . inches. able 


vale des 26.5 10 2.8 BOSS. 
42.3 36 1§.2 > 5263 >} CROO S20 Swede spa nriar 
38.8 28 10 Bi 1ST |= PAOD AT Ee of eae a 
47.8 47 22.3° | 886 | © 148 OO08 cde ees 
aati ate) 64 6.8 2 0.2 cae Sara 
30.3 €.25 10 73 ON SE BQO Ste erericy oh pee wos 
15.3 0.5) pees FOOT RAE CPF ae as atas 
32.0 4.9. fe. BB Sk CRORE nD ie aes 
18.8 5 20° TT ES B60 i eae 
18.3. 5 L.Os sf) AY <[- o= 6 CHO ee Bh Vics oe 
19-1 5 BO e180} 2 6 BaO essa) ae eee 
65.8 72 BPE 4 fh BZ adi: B16: TOD bi cate,2 5 elem 
13.9 3 5 aw ie ee ee ee Re ee Oe ork - 
39.1 3 11.2 SeKld- lee ROSS oe 
23.5 8 5S Pe | ae il ne OE ete aon Rll Ae 
29.8 13 359u (PSST 25.800 iy aes Ma ve 
35.8 21 1.8 E130" 40 800) cc: oo, Eee 


————— | cq |e qc“e $s) mm |_| —_— 


0 64 SEO sSh> G7 je Bor B60. Fcc fo a ee 
9 9 2.2 $4 BOD 0 SS e.g eee 
Tag. ct 2 6 1.2 7 970 
Se re 5 34 14.1 93 620 a 
ee 9 2 0.3 2 000 ee 
SPai 0 15 4.9 ele e's 
ermal nae. 8 3 0.4 teeee - 
aka cae B: 13 3.9 bl son | 
Ce ak ‘8 3 0.3 a ee 
eee 0 2 0.2 if Le ; 
eet. t 8 3 0.4 | 
ky en: 0 12 3.6 | 
Aan Fe 2 4 0.7 : 
a a 9 10 3. : 
me eee toe Tong fee | 3 4 0. : 
Beak ere 1905: ik Gert ae 15 5. 
ee ity | BEGG. NE | 14 24 se 
yer Ne ae 1807 ee: a 38 16. 
) SE SS Oe ee oe 0 8 2. 
A ed eet G00 wk 14 6. : 
AM VOLO. (is tot | BY 4. 
Py 1911 9 te 
,* PS RRR Gee .0 2. 
rh ee ara 6 it 
141s... 3 i > 
190508. 8. a 4. 
7 1916..$.5.., 2 9. 
19075. a 2. 
TOUS. & | 0 ay 
TF ae atetes 9 vir 
1920. 5 2. 
15) a a 6 a 
THOIA Noes 7 8 
NOOR ee ck, 0.8 2.6 
DORA IE ts at 0. : 
TOOK ea 2 5.7 0.6 ‘ 
TOG ry, 8 Se 
1927, 2 Le 
MODS too wets 4 0. ; 
T0500 <i be sa 0. 
1920-1930......,. ‘ 0 0. : y 
Summary: A> y 
AUOUARNicersemlne sib i. NEN Sy RT rs ep ay FD : 490 505 ~| 236-315 \ 
MEGAN wave stele ats z P 35 036 16 880 
Maximum 1883-84 91 650 | 103 000 ; 
Minimum 1876-77 J ‘ E 3 428 2 441 7 


ber lis 

istribution is such ‘that years having indices less” 
120 would produce large floods; those with indices as high as 130 t 

rf would almost invariably ‘produce large floods; and indices above 200 rep @ 
_ so much rainfall that major floods would be inevitable. oS 
i In computing probable flood flow quantities it is necessary, therefore, 
_ examine not only daily precipitation records but hourly intensities. Results 
Vy of studies of probable maximum 24-hr precipitation of a frequency of ‘one 3 
, zi in 50 yr are shown as shaded zones in Fig. 5. Like studies of maximm 
sustained intensities, in inches per hour, are designated by lines and quan-- 
tities also in Fig. 5. Wales 


ss A value for the possible 24-hr rainfall averaging 10.4 in. over the catch- 
_~ ment, was adopted as well as hourly intensities ranging from 1.8 to 2.7 in. in 
various parts of the same water-shed. Isolated instances exist in small zones : 

. of 24-hr precipitation in excess of 12 in., but it was not expected that this 
2 : quantity would be general over the entire water-shed. . erp 
Further experience showed that a critical 4-day storm may be expected, — 


i and that the run-off during that storm will not exceed 50% of the total — 


& seasonal run-off. Percentage ratios of maximum 4-day flows to seasonal run- 
a off for the same years, are: 

- : Year Percentage of ratio Year Percentage of ratio 

s AOOGe aoa 23 TE A ea 

4 1914..... 12 L926 es tes 35 

m UNG pees: 25 199 ey DE 


The flows from 1906 to 1916, inclusive, are derived by correlation with 
San Gabriel records, andthe remainder are from actual records. — 

In a determination of desirable flood-regulating capacity it is necessary to E 
know the effect of the lesser peak flows that will immediately precede the major =a 
flood peak. Daily discharges were plotted for critical periods during seven 
successive days on a number of water-sheds where rainfall-run-off records 
were available, which indicated that a 4-day hydrograph would involve as 
extensive a period as would materially affect flood storage; from such records 
the relation of successive peak flows was derived. The estimated distribution 
of run-off during a 4-day period on the Big Tujunga water-shed is given in 
Table 6. The successive peaks are shown on Fig. 6. The quantities refer to 
estimated flows at a point six miles down stream from the mouth of the can- 


yon at Dam Site No. 5. 


TABLE 6.—EstimaTep DisTRIBUTION OF Tora, Run-Orr Durine 
Four-Day Prriop 


=e 


Run-off, Rainfall, Run-off, Ratio: Run-off 
Flood. day in pose test in inches in incbes to rainfall 

a ues a ae 
Maximum day.....+-+++:++% 30 026 10.4 4.5 0.43 
One day preceding...-.----- 23 236 8.7 3.5 oe 
Two days preceding....-:--+ 13.316 6.7 are e-30 
Three days preceding...----- 5 884 4.5 : : 

72 462 30.3 10.9 0.36 


j mile; and the ratio of 4-day to seasonal run-off as > 2 OF 
, AL the values appear reasonable and quite readily basadaele d 
excess precipitation. oS nat Sere eeness Sean ends 
| Five reservoir sites are available in the canyon and below its mouth. — 
( peration studies were made of all combinations resulting in a decision that 
lating storage totaling 27 360 acre-ft would reduce the flood peak of 48 300 _ 
cu ft per sec to 14 500 cu ft per sec, at a capital cost of about $152 per acre-ft, 
uivalent to $125 per cu ft per sec, of regulation affected. This regu- 
lating storage can be created by dams at three locations, developing capacities 
3240, 7170, and'13 950 acre-ft, respectively. ene 


Erosion AND Erosion Controu 
oe HS Erosion, in this instance, is defined as the washing of loose material from 
_ weathered mountain and foot-hill areas, and its deposition upon the lower 
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Fic. 5.—Bie Tousunea WATER-SHED: ZonnS 
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Ii tt Ae 
steep) depend p 
I y on the condition of natural mountain vegetative cover. Destructioy 
_ of this protective cover by fires, or otherwise, will permit enormous déb: 
- flows to result from rainfall intensities which would create, normally, | 
small débris movement. This abnormal condition diminishes with re-grov 
_ of Nature’s protective covering, which is often a slow. process, since erosion 
= has partly removed the top-soil capable of supporting vegetation. ios 


Sad 


. Hazarps rrom Désris Movements s; 
Excessive débris waves first affect the lives and property of foot-hill resi- 
dents. Run-off from smaller tributaries responds more quickly to rainfall. 
than run-off in main rivers and, consequently, a débris-laden stream, upon © 
entering a larger main channel, deposits temporary débris cones, forming — 
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OF PRECIPITATION AND LINES OF INTENSITIES. 


mt channel under conditions of amnual stream How varying from 8 to 860 
moans, is most difficult; Seas problem is complicated by | Beard 
construction costs increase enormously and e erosion sae at ‘or 1 
origin of the débris | becomes essential. — as ise. EN ae: ; 


Fic. 6—CONSTRUCTION OF FouR-DAY HYDROGRAPH. 


- Recorps or Warer-Suep Fires 
Financial considerations govern in the matter of fire prevention and 
although Federal, State, and County agencies are entirely efficient, at present 
(1935), with the funds allocated, far from 100% effectiveness in fire prevention 
has been attained, as shown by forest fire records during the fourteen years 
since 1919. Major fire records are summarized® in Table 7. 


TABLE 7.—Burnep Water-SHEep ArEAs—1 000 Acres or MorE—1919 To 1933 


No. Year Area burned, No. Year Area burned, 

in acres in acres 

Je Sink 6 (OOM e ater 1919 60 000 HS eases cerectohs 1925 5 650 
erecta ctlaye Wns 1919 60 000 PLS Rote ete < 1925 4 000 
tS PaCS Co eno CREAR 1921 4 000 Dk. Se coke es wink 1927 4 630 
a yee See a Tee 1921 3 500 Lio eek ae 6 we 1927 9 000 
TR SS re 1921 3 300 tt eS er, 1927 2 960 
Oeste etiersis. ties 1921 1 500 US a ee es 1928 40 000 
Uo? 0 BS sae Se 1922 20 000 Se ge 1930 23 000 
aretha nce SRT GW. «s\'eis, ar 1924 52 000 AYs.. oUhwtoewoss oe 1933 4 860 
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What expenditures are warranted for fire prevention varies widely in 
different areas, depending upon population concentrations and values menaced. 
The foot-hill water-sheds of Los Angeles County are made up of many hun- 
dred individually small catchments, aggregating 285 sq miles. Of this, areas 
totaling nearly 60 sq miles, as yet not denuded, are located above incorporated 
territory having a combined population of more than 100000 persons and 
assessed values of more than $200 000 000. 


5 Rept. on New Year’s Foothill Débris Flood, by BR. C. Eaton, March 19, 1934. 
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» erosion is r latively slight. Quantitative records of erosion from 
A = les County water-sheds that have a normal cover are not as plentiful 
a8 accurate as those from water-sheds from which the débris flows are m 
readily measured; therefore, only estimates are available. is 
One estimate (cited as ample) of the extent to which the storage on the 
San Gabriel water-shed would be depleted, was 3200 cu yd per sq mile 
annually.” This value probably includes assumptions of occasional burnings 
since the estimate was made in 1927, within the 3-yr critical period after 
more than 30% of the area was burned in 1924. nea ee 
- The measured débris deposited in Sweetwater Reservoir, San Diego County, 
is reported as averaging 1400 cu yd annually per sq mile in 39 yr,’ and 
Gibraltar Reservoir, in Santa Barbara County, averaged 1380 cu yd per 
sq mile in its 1l-yr period.” A major fire in 1923, burning over most of the 
water-shed, makes this value high. Twelve scattered records collected by a 
the writer average 190 yd per storm per sq mile, ranging from 30 to 500 yd. 
No single rate is applicable to all the catchments of Los Angeles County 
and, although admittedly based upon insufficient data, it is the writer’s opinion 
that, with a water-shed undisturbed by fires for 10 yr or more, 1 500 yd per 
sq mile is a maximum quantity. : 
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a Désris QuaNTITIES FROM BuRNED AREAS 


¢ One quantitative example of débris movement is that given by the case - 
of Sunset and Brand Canyons® which debouch directly on the streets of Bur- — "4 
bank, a city of 16 600 population and more than $28 000 000 assessed valua- 
tion in 1931. Brand Canyon discharges on a part of Glendale which has a 
total population of 60 000 persons and a total assessed property valuation of 
more than $80000000. The catchment area of Sunset Canyon is 1.10 : 
sq miles; its length is about 1.5 miles, and its width % mile; and in elevation, ie 
it ranges from 886 to 3 126 (U. S. Geological Survey). Brand Canyon has a 
similar catchment of 1.03 sq miles. A fire that burned a total of 4600 acres 
occurred from December 3 to 5, 1927, and included all of Sunset and Brand 
Canyons as well as adjacent canyons. One year later (October, 1928), the 
first rains occurred, recording a total depth of 0.40 in. During November, 
9.10 in. of rain fell, of which more than 65% occurred in one day; a 0.59-in. 
rain of November 13 was followed on November 14 by a downpour with a 
sustained 1-hr intensity of 0.43 in., between 8:30 and 9:30 A.M., making 
a total for that day of 1.41 in. Intensities of 10-min duration were recorded 


as high as 0.35 in. — 
Canyon observers reported the resulting débris flows to have followed 


closely after the high-intensity period, with two distinct flash flows. The 


So Na ae Ee et 
® Rept. on Flood Control and Conservation, San Gabriel River, by F, H. Fowler, 
M. Am. Soc. C. B., C. D. Marx. Past-President and Hon. M. Am. Soc. C. B., and C, H. Paul, 


M. Am. Soc. GC. E., March, 1927. 
T“The Silt Problem’, by J. C. 


Boh i er 8 
_ te Rept. on Sunset and Brand Park Canyons Flood of November 14, 1928, by E. ic 
Haton-——-7- seme : : PaaS — ; : 


Stevens, Proceedings, Am. Soc. C. E., October, 1934, 


ay per. sq mile of water- ce a eel that wetted closely with as : 
truck removals. from. streets and lawns. Débris removal costs by public 
encies averaged $1 per cu yd, exclusive of Me; cost, to residents for removals. 
m cellars and lawns. 
Another example was afforded by the Arroyo Sequis, on the vente ce 
~ County boundary, which has a water-shed of 11.4 sq miles, discharging directly 
into the ocean. It is roughly triangular in shape and ranges in elevation 
from sea level to 3060 (U. S. Geological Survey). An intense fire from 
October 29 to November 6, 1930, burned 90% of this area, leaving its decom- 
posed granite slopes, of 300 to. 500 ft per thousand, exposed to rain. 

~The nearest recording rain-gauge was 12 miles inland, at Elevation 600. 
Since storms approach from the west they strike the coast water-shed an hour 
earlier than the gauge and it is probable that the actual rainfall on the burned 
‘- area was greater in quantity and intensity than that recorded. Rainfall 
es, records are listed in Table 8.. Beginning at 1:00 P. M., on January 7, 1931, 


TABLE 8.—RatnraLt Recorps, Arroyo Srequts. WaTER-SHED 


; 2 Cumulative 

Day Month and year ages rainfall, 

in inches 
ho Ockaber, LHSO. O80 ec ss. ccoma ls sie sitesi ens sie cee 0.16 0.16 
ty November, 193C 2.09 2.25 
~*~ | December, 1930 0.00 2.25 
= January, 1931-....... ce GL 3.36 
; January, 1931... 2... 0.47 3.83 

January, LOOT... cna 1.30 4.13 4 


and continuing to 4:00 P.M., the intensity for 10-min periods, expressed 
i in equivalent inches per hour, are listed in Table 9. The total rain falling in 
the 170 min. was 1.24 in., or an average rate for the period of 0.44 in. per hr. 


TABLE 9.—Arroyo Srquis Water-Suep; Ten-Minute Rarraut Intensities 


Time, Rate of rainfall Time? Rate of rainfall Tim 
Janery 7, | for 10-min periods, January 7, | for 10-min periods, J Sete fi ogres NolallA 
in inches per hour in inches per hour 1931 in inches per hour 


P.M 0.30 
123 0.33 
Sa EN 0.24 
P.M 0.36 
P.M 0.30 
Py Mi. 0.45 
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i) Quantitative measurements could not be made since much débris passed 
to the Pacific Ocean, but estimates based on mud marks on highway bridge 
piers (which were as much as 15 in. higher than adjacent bank marks), gave 
velocities of between 5 and 9 ft per sec. At least 50% by volume was solid 
matter and the combined débris and water flow was estimated as 700 cu ft 
per sec per sq mile of water-shed. 

A third example—Delta Canyon, a 0.63-sq mile tributary, of Big Tujunga 
Wash—has a particularly large cone of débris at its mouth (see Fig. 7). In 


Fic. 7—D£eris Conn at MoutsH or DELTA CANYON. 


elevation it ranges from 1 650 to 4500 ft (U. S. Geological Survey) ; its length 
is 1.9 miles from crest to mouth; its average width, 0.3 mile; and the grade 
increases from the minimum slope of 625 ft per mile at the lower end. The 
canyon forks 3 mile from the mouth, and at the end of the small fork, run- 
ning west, 1s a steep'mountain peak which is considered as having contributed 
the major portion of the débris. A large slip along the entire side of this 
mountain was first noted in 1913, but some sliding probably occurred previous 
to that datex The area is badly faulted. and the slopes are in an unstable 
condition. 

Fire destroyed.80% of the canyon cover on September 13, 1913, and, dur- 
ing the 1914 flood, débris was deposited at its mouth, shifting the former 
stream bed to the west. Although slips had been noted prior to 1914 there 
is no record of débris at its mouth prior to the rains. The heavy rains of 
1926 brought down additional deposits. No differentiation can be made 
between the relative quantities transported by the two débris flows, and a 


owever, ile eens transported in d 
, conservatively, at 123000 cu yd per Be nile, or 
ua yd per sq mile per storm. ‘ 


Sete most accurate, as well as the largest unit débris flow of authentic 


length of 6 miles and a total area of 7.5 sq miles. The catchment areas in 
this region are steep and rugged, rising 2 500 ft in 1.75 miles. The catch- 
, “ments range, in size, from Pickens, the largest, with an area of 1.6 sq miles, 
oe ‘to several as small as 0.2 sq mile. They debouch, through sixteen separate 
ig channels, on to developed foot-hill, city, and urban areas below the canyon 
mouths, aggregating about 7 sq miles. The area, which is mainly residential, 
extends from Tujunga City, the largest and at the western boundary, easterly 
‘through the unincorporated cities and towns of Highway Highland, La Cres- 
-eenta, and Verdugo City to La Canada, the easterly boundary. 
Kine aforementioned water-sheds were completely denuded by fire from 
‘November 21 to 24, 1938. They had not been burned previously since 1878. 
The early storms that began December 14, at 5:00 P. M., and ended on Decem- 
ber 15, yielded 4 in. of rain which packed the ash residue into the surface 
pores of the heavy soil cover, built up during decades of undisturbed brush. 
Recording rainfall records are given in Fig. 8(a) and Fig. 8(b) from Flint- 
ridge Fire Station at Elevation 1325 (U. S. Geological Survey). Following 
this preliminary saturation, came succeeding storms culminating with Decem- 
ber 31, on which date successive 1-hr intensities increased from 0.5 in. per 
‘hr at noon to 0.78 in. at 1:00 P. M.; 1.8 in. at 2:00 P. M.; and 1.14 in. at 3:00 
P.M. After 3:00 P.M. the rainfall subsided to 0.3 in. per hr until midnight 
when it increased suddenly to 1.28 in. per hr. A flash, 5-min intensity at a 
rate of 2.16 in. per hr began at 11:47 P.M. Three separate recording gauges, 
distributed over the area, registered sustained 1-hr intensities at midnight 
of 1.28, 0.85, and 0.87 in. per hr, respectively; and it is probable that a 1-in. 
intensity was representative of the entire area. No particular 1-hr intensity 
was a maximum, all-time record, but the successive cycles of sustained high 
intensity were without known parallel. 
Beginning almost exactly at midnight on December 31, 1933, and lasting 
for an hour, the resulting succession of débris flows caused 30 deaths; 483 
homes were either completely swept away or rendered uninhabitable; and the 
total property damage aggregated $5 000000. All the loss of life and 80% 
of the damage occurred in the eastern part below Pickens and Hall-Beckley 
Canyons. The largest city, Tujunga, was completely protected by a débris 
basin at the mouth of Haines Canyon (see Figs. 9 and 10). Another basin 
was under construction at the head of Verdugo Wash, the main collector 
at which the smaller tributaries eventually terminate. 
Falling on a completely saturated soil, the rain had caused sheet erosion 
over practically the entire mountain area concurrently with a slippage of 
numerous masses from steep canyon sides of weathered rock, into the narrow 


~The La Crescenta-Montrose débris flow, on 1 anuary 1, 1934, ree: the : 


—— record.* In this case seventeen contiguous water-sheds have a combined crest 
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Fic. 8.— RAINFALL RECORDS, STORM OF DrecymsBpR 30, 1933, To January 1, 1934. 


earth barriers thus created, to saturate them, they moved down stream until 
they either drained sufficiently to retard progress, cr were stopped by other 
slides farther down stream. As these movements were repeated, débris was 


collected in increasing volumes and the masses moved with increasing 
amomentum. Simultaneously, the side tributaries attained peak flows and 


the combined effect was to release débris and water from canyons in waves ~ 


15 ft or more in height. Reaching open territory these waves flattened to 


ie aa 
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5 or 6 ft in height and spread to widths of 100 to 200 ft scouring new chan- 
nels and thus picking up new débris loads. In canyons velocities were suffi- 


' Fic. 9.—Vinw or Haines CANYON D&BRIS BASIN, BEFORE STORM OF JANUARY 1, 1934. 


cient to carry along many 10-ton boulders. One 594-ton boulder was rolled 
on to a paved street at the mouth of Dunsmuir Canyon. 


Fig. 10.—Virw or Haines Canyon Diipris Conn, AFTER STORM or JANUARY 1, 1934 


. ‘ 
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-. In Haines Débris Basin. ....1...5...... 
In Verdugo Débris Basin..............+--0+5 92 000 
On property and streets.....................- 5385 000 


1A Oe RR ea eee eer Abies ... 659 000 


ra : 
_ The foot-hill and urban area below the canyon mouths covers about: 7 sq. 
_ miles; of this area 1040 acres, or 23%, were affected by débris (see Fig.l). 
a os _ Practically all the débris was deposited within the 1-hr period after mid- 
night, occurring in a succession of about fifteen sharp débris flows, moving i 
‘id e: ‘velocities estimated at from 5 to 10 ft per sec. Canyon velocities were 
considerably greater. Each sharp peak of major débris flow would be suc- 
ceeded by a rapidly moving stream which contained less débris, cutting a 
channel to one side or through the deposited débris masses. In some x 
~. instances these masses would drain out and, successively, become saturated _ 
by stream flow to a point at which they would again begin to move; and, 
encountering a rapidly moving stream, they would be carried along in a suc- 
cession of surges. It is estimated that the flows containing maximum per- 
centages of débris occupied about 30 min. ie 
* The east half of Haines Canyon (the westerly limit of the burned area) 
A had been denuded. Practically the entire débris movement from this source 
: - originated from 0.47 sq miles, indicating a débris rate of 67000 cu yd per 
sq mile scoured out by a single storm. Débris on property and streets and 
in the Verdugo Débris Basin came from 7.08 sq miles, equivalent to nearly 
- 90 000 cu yd per sq mile. ae 
An analysis of smaller sizes of débris collected from streets, is presented a 
in Table 10. 


TABLE 10.—Anatyses or DEpris IN Montrose SECTION 
(Average weight per cubie foot dry = 103 Ib.) ; 
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of water by volume, flowing at an average velocity of 5 ft 
tive quantities of débris and water discharged in the 1-hr period were 


pr. 


eam flows of higher velocities, which ground the stream bed far below its 
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gained from those witnessing the flows at the time of their activity. 
The writer’s conclusions were that a sustained 1-hr run-off occurred, which 


was equivalent to an average depth of 0.9 in., or 580 cu ft per sec per sq mile, 
and that the maximum débris flow contained 70% of saturated débris to 30% 


estimated as 48 acre-ft of water and 55 acre-ft of débris. 


/ 
i 


Désris ContTROL 


For major débris movements, positive protection requires direct storage. 
The best protection thus far advised is the débris basin, which is admittedly 


expensive and not sightly. It is positive, however, and its first cost is often 


less than that of one débris removal from streets. 

A case in point is the Haines Canyon Débris Basin (see Figs. 9 and 10). 
Its low capital cost of about $0.25 per cu yd was due mainly: (1) To its 
location in a Federal Forest area, consequently involving no cost for right 
of way; and (2) to the basin excavation created by an operating rock plant. 
Immediately after it was filled, the débris was removed with an average 
500-ft haul, for a cost of $0.30 per cu yd. Débris basins at more than one 
hundred locations will be needed ultimately. Their capital costs per cubic 
yard of débris capacity will range from $0.30 to.$1.00 per yd, depending 
mainly on excayation costs. Their physical locations are closely fixed, between 
the canyon mouths and the deltas, at points sufficiently toward the canyon 
mouths to prevent channel-cutting around them and as far down on the 
deltas as feasible toward the flatter slopes where excavation costs will be least. 

Their storage capacity is obtained, in part, by cut-and-fill excavation in 
which considerable waste is necessary due to the safety requirement that 
the deepest part toward the outlet or spillway should be in solid cut. The 
sides of the basin may be leveed from part of the excavated material. Where 
feasible, side-channel spillways are advisable, thus forcing the débris flows 
to enter from the upper end and follow circular paths, facilitating maximum 
deposition of débris. Sufficient adjacent areas should be provided for 
future deposition of material. 

Where water-sheds are unburned, the operating costs of débris basins are 
negligible, because little or no débrig removal is necessary. After a fire and 
succeeding high-intensity storms, the basins must be emptied promptly. 

The necessity of locating the basin largely in cut, results jin excavation 
costs ranging from 50 to 15% of the total, the other costs being for spillways’ 


st channel sections as left resulted from cutting by the later and receding © 


riginal grade and below that at the time of the making of the highest mud — 
marks. Computations of peak flows based upon these highest marks and the | 
final channel cross-sections must be modified by judgment, of course, as Pa a 
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Fie. 12.—Pickuns CANYON: PLAN AND PROFILE OF CANYON, JANUARY 1 1934 
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gene ally occupies limited space and requires consi ti 
construct. — On the other hand, the basin excavation covers a large 
acre to 2 acres, or more), and thus permits the use of large equipm: 
d rapid construction. Two large power shovels, supplemented by le 
ks, can excavate, on a conservative estimate, 20000 cu yd per three-sh 
_ day. Under emergency conditions, a basin for protection from an area : 
4 of 2 sq miles may be excavated in from 15 to 30 days. Although it is 
urgently needed during the 8 to 5-yr period following a fire, a débris basin 
_ will be inactive over the period of years when no fires have occurred. | At 
5 best, it cannot be made attractive in appearance and must be located not far 
_ from residence property. 

: ; These factors suggest the possibility of saving in capital investment by 
setting aside definite areas, preparing detailed plans, and constructing spill- 
_ ways, channels, and auxiliary structures, but deferring the excavation of pits 
_ until fires have occurred. The large excavation equipment can be used later _ 
_ for prompt removal of débris deposits. After mountain cover has regrown, | 
_ basins may be permitted to fill until fires again indicate that they should be 
"a re-excavated. 


Water Conservation IncipentaL To Fioop ReGuLation 


- Conservation of flood wastes will result from flood regulation for pro- 
“4 - tective purposes. a . 
Increased drafts on local water supplies have been accentuated by decreases 
4 in original replenishment due to paving and building activities which have 
changed former pervious surfaces to impervious surfaces. River channels 
have been improved reducing thus the original wetted areas. The over-draft 
from wells tapping the basins between 1914 (the last major flood year) and 
1931, has created a present underground storage space of more than 2250000. 38 « 
acre-ft. The effects of pumping in excess of replenishment is especially : 
marked in the basins of the coastal plain where serious intrusions of salt 
water have resulted. Here, water from wells on an area of more than 57 sq 
miles shows a salt content of more than 15 grains per gal, which is considered 
to be the limit for human consumption. 


Surrace Hoxip-Over ConservATION ExpENSIVE 


To provide conservation by surface hold-over storage under Los Angeles 
County water-shed conditions would run into enormous costs, particularly 
since hold-over of flocds for five years or more would have to be provided. 
High surface storage costs are here due to the steep slopes and large number 
of individually small drainages that are characteristic of the mountain 


< 


territory. 
As indicating the costs involved, the twelve reservoirs constructed or 


under construction by the Los Angeles County Flood Control District will 
regulate the floods from 400 sq miles, having a combined capacity of 109 000 
acre-ft, at an average capital cost of $200 per acre-ft. Under such condi- 


: ust 

1 e Bae: into” en asegsoned basins R ‘ 
0. unately, there are available unparalleled ae peel Se of high 
pacity. Created by the deposition of sand, gravel, and boulders through 
: the | ‘ages, these basins are twenty-six in number. Dikes and geologic ae 
‘i ave created underground barriers, thus forming a series of reservoirs which 
13) retard underflow, the water passing over and through these dikes into suc- 
es - ceeding basins down stream. Void spaces in the gravels range from 20 to | 
. 10%, the larger values applying to the main valley basins, the smaller to the 
MN coastal areas. py 
- The logical solution is to conserve the present wastes by vegaluctes and 
passage into the underground basins. Measured 1931-32 flood-water waste 
- to the Pacific Ocean from 1841 sq miles of major drainages, was 101 578 
- acre-ft, indicating conservation possibilities when flood regulation has been 

provided sufficient for percolation into existing underground basins. The 

distribution of this flood-water waste is shown in Table 11. 


TABLE 11.—Measurep Fioop-Water Wastes To THE Ocean, 1931-32 


Waste, in 
Waste, in acre-feet Waste. Spy 
* of = 


acre-feet per square 
mile 


a Los Angeles River 50 958* 47.9 0.94 
a : San Gabriel River 6 555 * IS.% 0.26 . 
oy ae Ballonal Oreelow rey. wan as hots hea Seopa re 21 785 194.5 3.80 4 
ue Malia ‘Greekyh: Sais tile boca tee ee ae 14 665 142.2. 2.78 | 
= IeweM SlGUEHS ... Seion sas SoMa Ad Ute Se a 4 029 60.7 1.18 4q 
iy Topanga Greek te. ck ee poe Ut 3 586 199.0 3.88 
4 PEOUAIC tte chee es ate mer eee cae ee Sn me ee 101 578 54.9 1.07 
: 
be * Includes part of San Gabriel run-off transported to Los Angeles River via Rio Hondo. 
“~~ 
4 SPREADING TO INCREASE PERCOLATION . 


Regulated flood flows may be encouraged to percolate in increasing 
quantities by extending the wetted areas on overlying gravel cones by “spread- 
ing.” A representative set of percolation measurements made in 1930-31 
and 1932 at fifty-three different locations is given in Table 192. 

Rates of percolation through surface gravels are mainly proportional to ~ 
areas wetted and only to a minor degree upon depth of water. Flood flows, 
even when the main débris has been removed, contain sufficient fine suspended 
matter, with collodial content, so that a depth of deposit of from 75 to 4 in. 
is sufficient to reduce the percolation rates 50%, or more, or practically to 
seal the surface. In stream-bed areas the “scarification” ef a few inches 
is sufficient to increase percolation rates markedly. 

Searification by plowing and harrowing of test sections on the San Gabriel 
River and Rio Hondo gave results indicated in Table 13. 

Increased percolation may be effected by the construction of “spreading | 
works.” In their design, ponding alone, or the creation of large shallow 
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receive a scrubbing action, and high percolation rates can be maintained. ‘e 
The velocities. necessary range between 2 and 4 ft per sec and must be supple- 
mented by flushing immediately before, and immediately after,arun, The quan- 
tity necessary to waste by flushing should not exceed 5% of the total 
quantity spread. Spreading works include head-gates from which main and 
lateral ditches take off through gates in effect similar to an irrigation 
system. The object of spreading is to provide a wetted area of water in 
motion covering as large a total area as possible. 


TABLE 13.—PerconaTion RatTEs IN THE San Gaprie, River anp Rio Honpo 
Seas 


PurcovatTion Ratss, In Cusic FEET PER 
Length of section, in feet Srconp ppR WeTTHp ACRE 


Before After Percentage increase 
Geeb feet genes sess Vee oh pear aa are AS aT SS 
Dei Medes 6 ae BOSS eA OO OTRO AROS oe A 5 ae 1.75 1.90 10 
OUGLOR Pee rca ere Oak ere. Mere n 2ee ety ose i" Af 0.52 0.88 69 
BAD 7 Oc ribts cath dialed iene te sccte eo sate ep eres oat 4.77 13 
DU Sree alee eo cieacPege exe aici Spite = nits . 
aC ees inet ei tsaceirrateue Doisitay tt fy Ee eye 1.99 28 
BrOBO ecie atteotea © acs usta tt eee ae” 2.72 5.50 102 


ging from $200 to $300 per eat acre of s 


ay range from $0.15 to $0. 50. per acre-ft of water percolated, 4 tk 


Ye ie oe : Concusioxs 


Bader physical and weather conditions as described the conclusions are, 


as follows: 


Flood Flows—On the larger mountains with vegetation undisturbed, a - 


the saturation period, will produce floods, their magnitude rising rapidly 
with increasing intensities. Twenty-four hour intensities may be expected, 
varying from 38 in. on the coastal plains to 13 in. on the mountains, with 
sustained 1-hr intensities of 1 in. to 3 in. Under the higher intensities, 
following saturation, flood flows of 400 acre-ft per day per sq mile are possible. 
With the characteristic sharp peaks of these areas this quantity will repre- 
sent a 1-hr peak of 400 cu ft per sec per sq mile. The regulatory storage 
needed will depend upon the economic balance between channel costs and 


total of 10 in. of rainfall will produce sufficient saturation to permit marked ; | 
run-off. Storms of intensities of 1 in. per hr and more, if closely following | 


regulation costs. In Los Angeles County its range is from 150 to 300 -— 


acre-ft per sq mile. 

On the lower foot-hills and valley areas with high percentages of impervi- 
ous surfaces, little or no saturation is needed to produce an immediate 
response of run-off to rainfall. It is necessary to anticipate future develop- 
ment both domestic and industrial in designing drainage channels in such 
areas and to set aside adequate areas for the ultimately required drainages. 

Débris Flows.—After a normal water-shed cover has been destroyed, erosion 
rates will increase from 50 to 100 times that with undisturbed vegetation. 
The best protection from débris is the normal vegetative cover. After 
denudation, any protective program is a secondary, far less efficient, and 
much more costly, defense.~ Erosion of a water-shed removes portions of the 
soil capable of supporting vegetation. Sufficient regrowth to afford protec- 


. tion will take from 5 to 10 yr, depending upon weather conditions and the 


extent to which the soil has been removed. With vegetation destroyed pro- 
tection requires direct débris storage; this can be afforded by débris basins 
at or near the mouths of the canyons. Following denudation, erosion may 
produce from 50000 to 100000 cu yd of débris per sq mile, the quantity 
depending upon slopes and intensities. The larger quantity may result from 
intensities as low as 1 in per hr. Little preliminary saturation is necessary 
and may be as low as 5 in. if occurring in a short period. 


ton and ihe Sriniee of spillways, ete., should be built, so that in curate 
 gencies they may be constructed promptly. od 

-Conservation.—Provision for flood regulation for protective purposes will he 
give capacities in excess of that required for conservation of the flood waters 
gee will result from increased percolation, which is principally in propor-_ 
‘tion to the stream-bed areas wetted. Where increased percolation above that 
3 of the natural channel is needed, spreading systems may be built at costs of © 

~ $200 to $300 per gross acre, with percolation capacities of som Ito 
~ aere-ft per day per acre. 
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_ AN ASYMMETRIC PROBABILITY FUNCTION _ 


Y 


Discussion 


a . By J. J. SLADE, JR., Esa. 


4 J. J. Suave, Jr., Ese., (by letter).“*—In developing the function presented 
in his paper the writer has necessarily had to give considerable weight to the 
% partly bounded function because it was by means of his analysis of it that 
_ he was able to construct the elements of the totally bounded curve. In his 
“a work, however, he has not had occasion to make use of it. This former func- 
‘ tion is merely a limiting case of the latter, in which one of the bounds be- 
_ comes infinite, and in the problems that he has considered the bounds, 
* obviously, have been finite. Perhaps one of the reasons that much of the 
work of statisticians has been unavailable to the engineer for hydrological 
investigations is that statisticians are interested usually in graduating the 
middle range of a frequency distribution, the extreme outer 1%, say, of | 
- the distribution being of relatively slight importance to them. The hydro- 2 
 logist, on the other hand, is frequently interested mostly in this outer 1% . : 
of the range. It is for this reason that one must insist on the determina- 
tion of the limits of the range. It is absurd for a probability function to . 
give a finite probability to a Noah’s flood. In modern times, there simply is a 
not enough energy resident in the atmosphere to produce anything like it. 
Equations (45), as given in the writer’s paper, are based on an easy gen- 
eralization of Equations (44). It is apparent, at once, that the factor, p, 
belongs therein and that it is not merely “a kind of artificial tour de force,” 
as Mr. Fisher expresses it. In the paragraph following these equations, 
the writer states that their form will probably need modification. A more 
extended analysis has shown that this is indeed the case. Their correct form, 


when b < g, is: 
:=peloea (+2) Aix. ga uielinndealecte AGeou 
g—2« 


goth 


Notu.—The paper by J. J. Slade, Jr., Hsq., was published in October, 1934, Proceedings. 
Pieenssion on this, eines has appeared in Proceedings as follows: January, 1935, by 
Messrs. Gordon R, Williams and H. Alden Foster; February, 1935, by Messrs. R. D. 
Goodrich, and F. T. Mavis; March, 1935, by L. Standish Hall, Assoc. M. Am. Soc. C. E.: 
April, 1935, by Arne Fisher, Hsq.; and May, 1935, by Arthur WwW. Kempert. 
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Sia Received by the Secretary July 22, 1935. 


/ 


via: 


“ 


In his discussion Mr. iia states that “any purely mathematical ‘teal 
ment of stream-flow data presumes at the outset that all the data are equally 
reliable,” and, from this, he infers that there are limitations in the analytical 


method which are partly absent in the graphical treatment. The analysis of 


a physical problem begins when all the data relevant to the problem are 
known, and these data necessarily include a knowledge of the reliability of 
the observations. The purely mathematical treatment begins, first of all, 
with the proper weighting of these observations, for which various methods 
are well known and widely used. The writer can not agree, therefore, with 
the idea that the analytical method is inferior to the graphical because the 


free-hand curve “can be shifted to give more or less weight to certain points.” 


Mr. Williams’ objection to listing the complete set of peaks in a stream- 
flow record is well taken. In the first place, it is difficult to define a peak, 
particularly in the range of small flows, as he observes; second, it is difficult 
to estimate above what datum the peak should be taken (that is, statistically, 
a peak is an accidental irregularity above a smoothed hydrograph, and it is 
difficult to smooth the hydrograph properly) ; and, third, peaks are not uncor- 
related quantities and, consequently, a simple analysis of their distribution 
is not possible. In presenting this paper the writer limited himself to the 
analysis of the generalized probability function, and he did not attempt to 
offer a method for the analysis of stream-flow characteristics. The illustra- 
tion of the peaks of the Tennessee River record was given merely to show 
that an extremely skewed and platykurtic distribution could be closely fitted 
by means of this curve. Using the corrected Equations (88) and an esti- 
mate of the extreme flood (415000 cu ft per sec) based on the following 
method, a much closer fit is obtained. 

The use of the partial series of peaks above the “one-year flood” is open, 
at least statistically, to a fourth objection, namely, that an error, which may 
be serious, is introduced in assigning a probability to the “one-year flood.” 
Because of the indeterminacy of this error the writer has not been able to 
make an ‘estimate of the reliability of the probabilities obtained from a curve 
fitted to such series. 

The writer concurs with Mr. Foster’s opinion that Pearson’s curves satisfy 
all the desiderata presented in his paper for a generalized frequency curve, 
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100 items) are so unreliable that they are useless. This statement nee 


pared with the wide fluctuations in the latter. £7 


Table 22 meets the objections to the labor involved in graduating dis-— af 


p tribsttions by means of this function. If greater detail is needed than is 
_ given in this table, then, by perfectly straightforward computation and the 


if fill in this detail. Tables of the probability integral have been published fre- 
; quently and distributed widely; Mr. Foster’s statement that these tables are 


“not always available to the practicing engineer” and that “even if avail-_ 


r able, might easily lead the average engineer into difficulties” is, therefore, 
_ quite disturbing. The answer to this objection is probably that the “average 
3 engineer” should not attempt, with his small equipment, to analyze a highly 
, specialized and difficult problem. Undoubtedly, it is true that some kind 
_ of an answer may always be found easily to any statistical problem, but it is 
in the interest of “practical purposes” that this answer should be questioned 
and its reliability tested. It is unfortunate that statistical analysis does not 
offer many short cuts, particularly since a wide variety of important problems 
yield to no other kind of analysis. 
Both Professors Goodrich and Mavis disagree with the writer’s statements 
- regarding the undesirability of graphical methods of analysis. These state- 
ments the writer intends, of course, to apply to the graphical analysis of sta- 
_ tistical problems specifically, because he is well aware of the value of 
- graphical analysis in general. Without entering into elaborate arguments, 
his objections to the graphical method may be stated briefly, as follows: ike 
an n-parameter curve is to be fitted to a set of data, then exactly n inde- 
pendent statements must be used to determine these parameters. If a statis- 
tical series is capable of furnishing four independent, significant character- 
istics (for instance, the mean, standard deviation, skewness, and kurtosis, 
or, as the writer prefers for small samples, the mean, standard deviation, and 


the two end points of the range of variation), then these statements are 


sufficient to determine reliably a four-parameter curve; more parameters may 
be determined by making use of moments higher than the fourth, or by 
the method of least squares; but if the significant characteristics of the dis- 
tribution are only four, the added work is meaningless; that is, there is a 
definite point beyond which the flexibility of a curve is undesirable. Now, 
a free-hand curve is a many-parameter curve (theoretically the number of 


no mathematical proof (although such ‘could be readily given); the experi- « 
. ment may be tried of drawing small samples from a fairly large population, — 
and, from these samples, the mean, standard deviation, and the skewness may i 
be computed. It will then be found that, although the first two character- f 
 isties will vary, they will show a remarkable uniformity, nevertheless, com- bs 


4 use of tables of the probability integral, Equations (88) may be utilized to — 3 


ry i 7 a SS LL ne ae, ce 

: 68 \1— go'I—| 891] 29°I—| 99°I—| oo'I—| eo'I—| 19°] S8hI—| TI] 71: 
S 9 I—| 9S I—| Le I—| 2a] 21] 22'I—| 22-1] Let] 82'I—| o8'I—| vet] 17rt7! 
£0'I—| ¥0'1—| 0°I—| 990'I—| 90°I—| 90'I—| 240°I—| L0°I—| 80°I—| TIT] STI) 967T-] 0TTT! 
8'0—| 98°0—| 98:0—| 28°0—| 18°0—| 88'0—| 68'0—| 06°0—| 16°0—| $46°0—| IO'T-) ATI] 77177) 
¢c'0—| s78'0—| 99°0—| 99°0—| 29°0—| 89':0—| 98c'0—| 09°0—| 29'0—| 99°0—| 84°0—-] 26°0—-| 71" 
-62:0—| G8z'0—| 6z'0—| 0g'0—| I8°0-—| ze0—| ze'0—| Fe'0—| 9e°0—| 98°0—| S9F°0—] ZL'0—-] 777 "7! 
B0°0—} 20'0—) 80:0-—| #00} 20'0-) 90:0-| 40°0-— 80'0—| O1'0—| €1:0—| s81'0—| r'0—| 77"! 
wae areca occ | 4270 | 00 | 02°0.| 0G) Gt'O | 41'0' | 9170" | et:0 | wO'O | 17777! 
oso | ogo | 09:0 | oso | 6F'0 | 6F'0 | 6r'0 | 8h'0 | 8F'0 | 8h'0 | BF'O | BHO | TTT! 
‘0 | s8'0 | sso | 7:0 | #8°0 | #8°0 | 98°0 | 98:0 | 98°0 | 880 | 260 | OT | 77°"7: 
ma joie Tt) |ocet | cet | eer) | ye | get | Zot | oe't | wet | .eo:n | 99L°1 ee 

@ oe | srs | to: | a's | occ | soc | coo | 90° | OLS | 94'S | BA's | 77°! 
| wee | vee | see | tee | cers | oye | cure | sre | ere | Tee | os | ooie 
% | wy | 920% | OL% | ScI'e | Gzt'y | €1'e | T1'F | GO'% | SoBe | BSE | 26S | To”: 
+ | 99°> |ce0% | oL% | 1L°6 | OL: | 99°F | GOGH | GOh'H | OF | GLB | 66'S | oI”! 

Td se a a 86°) | LL'% | Teh | ¥8' | Ge6' | °° 

g | ox. ZL | 2L°9 | g9°9 | 9699 | g¥'o | O89 | S009 | 999'y | 6872 | OO7e | 77°": 
Pore cr ees) #18 | 909 | e6°9 | 82°9 | Go°9 | BI’? | 29°% | #6°2 | 0078 | °°°°°! 
gain 09 Onl ee Gui veo | seo | 9¢°9 | 9o'9 | oz'e | ceo | e4'% | gee | one 1’ 


- Z=6 (9) 

2 98" T=—|-<-99° T Zot) 89° | 69 I) et) en 2I—| se°t—| 08 I—| 988'T—-] ct ie o6 
io Ce eee eee ce te Set OSe =| BRT) Fe I) SOT) TL Tp en | eT 06 

90 oj te Z00I—=) S01), 60 I) Bhi sl I—| 91'T—|) Oc I—| O8°T—| O9°I—| te Pc ctr cs 
8 /93°0—| 98°0—| 28°0—| 988°0—| 06°0—| 76°0—| 76°0—| 66°0—} GO'I—| ABI | In 08 
eg" ee°0—| 72°0—| F9°0-—| so'0—| 99°0—| L29°0—| 69°0—| £9°0—| 0L'0—| S96°0—-] | i OL 
Se Oz O02 0—| 96 0—| Jz'0—|- 46°0—| 82:0—| 6c°0—| 08'0—| S¥8°0—|, Br'O—-| | tt: 09 
oe a 0 0 0 0 0 0 0 0 Oo fetset te | eee [eee ecee see nes 0g 

‘ me OG Or |" 960 96'°0 L120 160 82'°0 66 °0 0g'O | SFE'O Sy cOrateee ce ie G en  e aeaae OF 
eee oO #0 rg 0 gs°0 9¢°0 19°00 690 £9°0 OLC0 MIR eeG Ot ie ee e A e M e 0£ 
| ¢8 88°0 98'0 248°0 | 988°0 06°0 660 $60 66°0 60° T ORM ah te ee oh. ot. ee eae 02 
ve oe 'T Te ‘1 GET 8° T cert 96° 1 68° 1 tT eo°L PARI gik ee ed oe Ri as eer ea Ot 

“ 18% 18'S Tes 08'S 62% 82'S 92'S GSS SI'S (EME ce |S = aS i mee i a Maa” oT 

| 80° 10'€ 86'S 66'S 68 °3 8s 81'S 69°% co's FES OO tier Bena WER > “iter AN ga ee ol TO 

sg" go's og sé CPs 0g s Go & OL € 6S OLS (Gana NE Ch He SF SPs ek” AR ee ee 10°0 

| +0" | 86°78 06's SL's 198 8h 'E GE'S Ole 68% 97'S CAC Gee me Rigi aie ae eed 100°0 

CP cs 7 FOP 10°F HE 19° Lv'E | GOS'E 83'S 8h'S OD izmir ot EY ge Bat ss 1000 °0 

92 is 99°F Gchye| 18h 0 °F 68°& 8G°€ 86° 66S G8h'S WONE ie at eee ae eee 10000 “0 

al- 90° i t6'F 9L'F yh 9L'F F6'8 99°€ ce € $63 6h 'S 032) ie ee ee oe Sa eee 100000 °0O 

oe" 8g 96> | 89% 82 F 10'% €1°& Le°€ 96'S 09°% (Oi a ton lel (a hie a= 1000000 0 

T=0 (2) , 
eee ee Obed Naga ye ol eel. 8 ahige 1.6 ko ogi | tT 
a —s = . oul} 
t= Y ‘GO01VA XEN] WHS MOU ALIIAvaOUg -7o-% 
oni e “NY GNV @ JO SHNIVA SOHOLUV A HOX 7 AO SHNIVA—G6E WIAVL 4 
a eres 
7 ¥ . ars 
a = te woes | ee 


—| 18:1) 09:I—| 8rii—| srI—| eer 
€S'I—| SsI—| 2o'i—| 121I—| 6tI— 
£0'I—| $0°I—| seo'1—| ¢0°I—| zo I— 
98'0—|} 28°0—| 28°0—| 28°0—| 28°0— 
1¢°0—| 6¢3'0—| 6¢°0—| 09°0—| 19°0— 
ze'0—| e8'0—| re'0—-| se'0—| 2e°0— 
— 20°0—| 60°0—| OT'0—| IT'0—| &1'0— 
“0 61'0 91'0 91'0 STO €T'0 
a) L¥°0 SPO 770 +r'0 oO 
0 €8°0 T8°0 | ¢08°0 180 08°0 
‘T FET ee rT GET | $9e°T 
“4 99°% 69'S 71'S 64'S | 998°S 
e GLa SL's 98°¢ 76'S | ScO'F 
oh $9'F IZ’? | SI8'> €6°F 10° 
2 979 | ggg | Ggo'¢ oss | see's 
a4 Gzs'9 ze'9 cP9 1g°9 | S29°9 
9 26°9 | StO°L | SOT'ZL | SZa'Z 082 
9 09°2 89°2 082 83° 2 G82 
06'8 828 |-o2e°s | Son's Ges 
b=0 (P) 
eo'I—| e¢'I—] of I—| 6rI—| oF I— 
¥oI—| Fo I—| So I—| ss I—| so I— 
¥O'I—| SO'I—| 7O'I—| F0'TI—| FO'T— 
98°0—| 28°0—| 28°0—| 88'0—| 68°0— 
1S'0—| 8¢'0—| ¢8¢'0—| 09°0—|} ST9‘0— 
Teo—| e80—| e8'0—| se'0—| soe‘0o— 
90'0—| 80°0—| 80°0—| OL°0—] Zt°0— 
61'0 810 L1'0 91'0 +10 
8F'0 L¥'0 L¥'0 SP'0 $0 0 
0 £8°0 | 928°0 | sz8'0 | ¢z8'0 | sz8'o | Szs"0 28 °0 0 
ise T20 -1 eVEt | 98es pe rT Ge'T Le°T oh 
IG '% Ho's 89'S 29'S 99'S 69'S | SeL's 61° 88'S 86'S ore ae : Le GP 
ave Bhs 78 £9'¢ 19'S TL’ 84° 98°8 36'8 96° ee Soe ape asa 
61h 08'F 98°F LY'Y oS 'F 8o°7 | SZ9'F 89°F OL | 98e'F PLP ; ee 10'0 
| 98°F 00'S 80'S 06'S 93° 08's ge'g | ogee | g2z'G | S66°F 08'F #: eee TOOK 
8b'¢ cog gL'g 98'9 | 26'S G6'S 86'S 06°9 |-S%Z°¢ LZ'S (LET ee) | aor Dy eT OOO LO) 
| $20°9 Gz'9 Ge'9 oF'9 o¢'9 | Sze'9 | SLPr'9 989 10°9 9P'S eh Pp me Mie Ne TORONTO 
¢09'9 8.9 069 00'2 | 820'Z 00°L | 926°9 GL°9 | 9ze'9 09°¢ cory |: “= **""T99000'0 
204 eZ, | 9er-d og'L 03°24 Gp, 1 ¢¢e'L | 9S0'2 | See'9 OL’ tao las Te ap |ae *T000000°0 
oe of 2 e=0 (0) 
ee see eee epee vy ore Hoek gh) Nee ge et | ot 
- - ourty 
i= YX ‘HO1IVA XGON] AHL WOU ALITIAVEOUT -jo-% a 


*(penuyu0g)—ese ATAVL 


SE AUR? Ol 


= 
8 
Soococcn 


Zoi) oF I—| 8FI—| oF I—| ehI—] Ge'I—| Pe'I—| 99c'I| STI] HO. 
¥o1—| Ze'I—| g1z'I—| oz i—| 61°I—| 9tI—| FI'I—| OL't—| e0'1-| 68°0-] 
‘t-| g0'I-| 80'I—| Zo‘;I—| 20°1—| I0°I—| 00°I—| 66:0—| 86°0-) ¥6°0—| S8°0-|) 
‘o—| 98:0—| 998°0—-| 98°0—| 98°0—| 98°0—| 98°0—| 998°0—| 98'°0—| ¥8'0—| 180-7) 
‘O—| 19°0—| 29°0—| ¢89°0-—| 69°0—| 09°0—| 09°0—| ST9'0—|_¥9'0—| 999°0—| ILO 
‘o—|  Te:0—| aze'0—| ve'0—| ge'0—| 98°0—| 2e°0—| 68°0—| &h'0—| 80] 62:0] 7 
-o—| 10°0—| 80°0—| 60°0—| 90T'0—| ZI°O—| SEI°O—| 9T°0—} 1z'0—-| 82:0—-| FO 
' oo 1 sto. | ZNO | 910 | 91:0 | 81:0, |..ct:0 | 60°0 | O'0 | WO] FO | 
0 e170 | 2¥'0 | 997'0 | S¥'0 | 0 | G0 | GoV'O | 88:0 | ee'0 | 9:0 | WO TT, 
0 0 | aiso0 | 180 | 1870 | 080 | 62:0 | 982°0 | 84:0 | 24:0 | £2°0 | 29°0 | 6F'O J 
T Bh | Mar get [over | ee't | vel | Vel | 98'T | get | get | ob fire: 
°& iz ‘2 | cou | o:% | 4:2 | 82's | 98's | 86's | Go's | ST's | eee | B9'E [ori 
a ¢ |stoe | soe |ocse | ose | 66's | SOt'y | 20°F | OFF | BH | 966% | GGFR Pe 
F + | 09% | 99'% | 918% | 986% | 920'9 | 990'9 | Bro | 84:2 | 920'9 | 81°90 | S099 | 
‘9 ¢ | te¢ | ogg | e249 | 06'¢ | 80'9 | Ze9 | 969'9 | 06°9 | OTL | 869 | S4'o Pe, 
$ ‘¢ | so'9 | 1e'9 | 00:9 | 08°9 | co's | O82 | 9892 | BBL | BOL | Ses | oLB eg |. 
“¢ G@ | eo | Lox | Ser, | 9902 | coz | Oc | O98 | C18 | Zs | 8612 | 9869 |. 
‘2h Z| ween | sh | 128 | OFS | 9r4'8 | 90° | G26 | M776 | 826 | Ses°8 | 969 PT 
L 1 \oats |ones | o68 |ope6 | 296 | 986 | oror! ocorl 46 lens | 269 1 
I—| 9¢'1—] 62'I—| 21:1] 96°:0— 
t—| 9t'I-| Zt'1—| 90'1—| £6:0— 
0 T—| 10't—| 10°1—| 66'0—| 96°0—| 06°0— 
| 98° 98°0—| 98°0—| 98'0—| 98°0—| 28°0—| 98°0—| 28°0—| 18°0—| 18°0—| 98°0— 
o¢0—| 9¢°0—| 29°0—| 29°0—| 82:0—| 69°0—| 09'0—| 19°:0—| £9'°0—|, $9°0—| 69°07) 410 
ezo-| of'o0—| 1s:0—| ze‘o—| ee'0—| ¥8'0—| 9e°0—| L8°0-—| OO] &h'0—-) 0S'0—| 99'0— 
0° 0— 10°0—| $10°0—| 80°0—| OT'0—| ZI'0—| eT:0—| 9T'0—| 02'0—| 6z:0—| I2'0— 
02'0 60 | 20:0 | 210 | Sto | €1°0 | 2t:0 | of'0 | 30:0 | wO'0-| 08 0— 
60 | 170 | 97:0 | Se:0 | 0 | ZO | STv'0 | 68:0 | ge:0 | Lz:0 | 0 
€8°0 | 9180 | 80 | 18°0 | 18°0 | 08:0 | 64°0 | 64:0 | 44°0 | 94°0 | 69'0 | 6F'0 
een | eet |} eer | eer | Set | per | wet | 98:E | Set | Sze'r | Set | oe8't 
oo'z | soc | 29:2 | 89's | oea's | 24°S | 9e8's | 06S | Zoe | Ste | gee | 2978 
See ieeee | jeoe | mre | ose | roe | toy | oly | 2e7 | us'% | Fey | S6FF 
lore | Gee | 098 | GA:y | Gee | oe | stg | Tes | e's | OL'e | 29'S | 08% 
ce | oe | io | sree | O29 | ceo | 60:9 | 18:9 | Ze:9 | O99 | 2:9 | IGF 
9g'¢ | 10:9 | 61:9 | ¥°9 | 09°9 | 08:9 | 926°9 | Oc:2 | S182 | S2e'2 | 49°9 | 86% 
#0 | co'9 | 26:9 | ooz | IH2 | cod | 642 | 966: | 20'S | BL | O89 | 66% 
que | eek | 69:2 -| $6 ers | 9e:8 | 99'S | 04°8 | O24°8 | SoS | 90:2 | 266°% 
gra | 60g | ec8 | sos | O88 | coe | tee | ee6 | oe'6 | sus | sors | OOS 


omeper ered fet foee bey it ae) (et eel 


© j= X ‘GNIVA XUAN] HHL WOe ALITIAVaOUT 


+ 


— (ponuyuon)—ee ATAVI 


it 
gs 


a 4 
SRESASSSSS 


Rrewwetionnosoe 


eh 


* 


SANSNBBRAT 


. 


betel atata eae bik 


S 


al 


bas 
iw 


SSEHeease 


BREA 


ice) 
1) 


te) 
ito} 
x 


pt] 
a 


Tit 


ON GSESHTHAnOSS 


gees. 


ie} 


SRSRIBSSRBIASaS 


* 


TTT 


ON GOSBIAIAROOS 


MANOON 
THOR OAD 


rea | 


BegStee: 


AMHKSO 
ZSS528 


1D 
WOROATWHHPOWoO 


TE TE 
Tt oo 
‘QE -o ie os 
pee -G— fei AS ire a a 
Oo apa Soe) woo tO we ee 
98:0—| 18°0—| 0F-0— Q— eso ae 
o> Sei0—|"- FE 0—| - 25 0— ee (0 RO alae eee 
0 eL'0 IL'0 80°0 0 ‘o— Qo 
0 £70 1¥'0 or 0 Le°0 ce'0 G3'0 LOROR NS ate Boe, | 0g 
0 64°0 84°0 L1L°0 g2°0 1L°0 ¢g9°0 Seer eras 02 
T ee'T od rer Pe Li | SESiL CoplMecoeuL, (tee: atest OT 
4 18'S 88°S | $163 80°€ FOE LYS FE ae ee RL lad OT 
€ 80°F oo F Iv? co°F L6°¥ ges ORGS alii dacs as eee TO 
*g | See's cr ¢ G1'9 40°9 | gt9°9 Se ea ec) AS) a a a 100 
‘9 GE '9 09°9 | $246'°9 CFL 06°L Che eel eG, Ae ek ks 100°0 
L oe°h TLL ST'8 $98 OT 6 90°6 GOR esas tke: T0000 
‘8 | gg8°8 81°83 66°6 | $42°6 GTO | $L2°6 PA Ne seks: **T0000'6 
‘8 08 6 816 ASLOTAMGLLAOL | OOrETo|) OS20T Io 68. |. 01s e 100000 0 
6 92°OL | 92°OT | O€ TT PTT lee Clr LE jes Ou sOb SES 2 (oti 1000000 ° 0 
8=o0 ty) 
Gr i—4. Zr 1—) “seI—) - e6:1—| sse-I—| ef I—| S16°0—-\"* oe 6 
61 i) Skt) OLI—| «ett 60° I—}| St0'T ATS cel ae ees a eres 06 
Z0°I—|. 10°I—| 0o°I—| 66'0—| 26°0—| 26°0—| 88'0-] 8 
98'0—| 98°0—| 98°0—| 998°0—| $8'0- €8°0 SROa2: cu bebe ke cee 08 
63 '0—| S6¢°0—| 209 0—) ¢19°0 £9°0- ¢9°0- BO Ooo bree ae OL 
ce0o—| eoe'0—| 2e'0—| 6E°0-| crO—| 2F0—-) O90) lt tle: 09 
Tost 0 | 0 — | 9b 0— |. 1s O19 AS10— |" TPO fe 0¢ 
0 F1'0 €1°0 Ilo 60°0 | $80'0 BORO ra seh: Onalame e w e OF 
‘O | s8%'0O | SIF'O | S0v'0 8€°0 €e°0 &z'0 SOPO: Syme ep a Se 0€ 
‘0 08°0 SLO |-SLL°0. -| $94°0 6L'0 990 Fie Se | esp ae em 0G 
\ eer ee. cet | SPe'T ve 'T ee T OGG) Are dace. ORR OT 
6 08°% 18°3 $63 90°S G'S vrs SE EN Seo NTN: OT 
$ G0°% LV} | SPE 'P L19°% 98°F 06'S ATOR Paling yet ois 10 
oe] LI'g 98°9 69'S G6'9 0g 9 OBO Cy GON Oe na ce co BS 100 
9 | g&'9 LP9 08'9 | 812 Hod Oke SE828) te T00 0 
9 ToL | Gee'h | 9682 T8'8 09°38 ges O20) Wet. y ***1000°0 
4 | SLT°8 Is'8 £6°8 ce '6 8h 6 168 SBTO VN a T0000 °0 
8 80°6 S¥'6 | $186 Pe Olle 0c Ob | 08:0") GNG.8 Wien: ** 100000 °0 
6 66 SOT | GL°OL | €0'TT 8°Or | 096 ¥6°9 | °° ~*°T000000°0 
L=0 (9) 
eae | eel NEO e's lt see eT 
; — oulry 
‘TO1VA XMANJ THL WOU ALITIAV GOUT -Jo-% 


*(penuyu0g)—ee WIAVL 


ag 


8883 


BER 
19219 
gugSdegqeeaase23s 


«(| SaS8R888 


; ¥°0-| 

| 0g 62'0— 
520) SPO '0— 
ige-o" | 0S'O 

: @P'0 | GEO: | 
zs'o | €8;0 

| get | zevr 
ges | 69's 

| eg'e | og'e 
core | Lhe | 
gst'¢ | og's 
zo'¢ | 20°9 
-o9'9 | 289 
2EL | O9'L 
162 | 61's 


g 
aa 


€o°I—| 09't 
SS. T=) 1S T— 
CO}, 4.60 | 
¢$8'0—| 98'0— 
OS Op LEO 
Te 0—| ceo 
90°0—| 40°0— 
61°0 8T'0 
L¥'0 9F'0 
S18'0 18'0 
set | 98" t 
9% 19°S 
49°8: | 94'S 
19° | WSL°h 
LYS | 69°9 
soe '9 Lo’ 9 
ar ehh 
S8°L G28 
09'8 90°6 
fj; 
aa], OF T— 
Fo V—|, 46. T— 
COL} 0. 
98°0—| 9F8'0— 
oc" 0 1980— 
Re Gea Ge, O— 
20° 0} -20°0— 
8T'0 810 
Lv'0 v0 
18'0 e8'0 
Te'T Ge'T 
19'S 89'S 
99° LL 
S6o'F 91°F 
FPS cog 
ga'°9 gg°9 
40h 8e°L 
SLL L GS '8 
see 8 86°8 


$19'0—| 


8 


te) 
PMO QDIANNOIWOS 


HDRBRAHHAMMDGA 


mopenibinetis " 


Sj 


BANANA 


osocsocHn 


‘HO1VA XUMGN] GHL AOd ALITIAVAOUT 


oer'I—| IF I—| 93°'I—| 18'I—| es I—- J (E(k Sesaioee 

Brats a St 20 986 0 880 oe 
To°’I—| 0o'I—| 66°0—| 86'0—| 96°0—| 968°0—| 984°0—|" pee 
98°0—| 9c8'0—| ss'0o—| #8'0—| £8'0— : ‘TART Ve ae a a 
62'0—| 62°0—| 09°0—| 19°0—|} 229°0— OO ays ae pe ie 
9¢'0—| 98°0—| L18°0—| 6&°0—| ZF 0— (7 | a tea ai 
Il'0—| Sat‘0—| SFL'0O—| L21'0—-| 12'0— Pee Oanle cer rie 4 
#10 £10 O1'o 10'0 £00 GET Oe tee 
£b'0 10 680 92°0 Te'0 “| exact, Spee 
8L°0 8L'0 92'0 %L'0 0L'0 2 nae lip nore 

cet eet | $ze't £e'T oe i! Eis olga 
€8°S 16'S 86°% IT's 92'8 i ee EXO: 
Ob | G12'F OF $L'F 60'S 9 ory ST 
9ze'9 | gua‘ 98'S 62'9 | $z8'9 Re eee 10° 
obo 08'°9 | GZT’L CLL | S38 Oe a Rr ae 100° 
gg’ 66°L 09'8 8I'6 68°6 EO Wer cece yok 1000° 
99°38 | SSI'6 9L'6 19°01 | OZ IT ae AER eS T0000 
89°6 62'OL | F6'OL | SLIT | 98'SI ‘6 *****190000 
OL'OLr | 98'It | L20°SE 6°21 | 98°ST Biraleence s 1000000 

oT =a (f) 
Hr i—| SIh'I—| LeI—| ce i—| Fz oe ti CN ep as 
6L'I—| GZE'T—| St't—| st't—| 20° ‘o-— Pips Pe inecisionic sys © 5s 
“o'T-| 10°T—| 66°0—| 86°0—| So 7 oe eee 
es'0—| 98'o—| 98'0—| s8-0—| | &8° ‘o— sp |eceeeeeresenes 
ges-o—| 09°0—| 09°0—| 19°0—|” 20: a jl eee 
92'0— 92° 0— L2°0— or‘ 0o— 4a ‘o— sae (Or eae Bt 
To—| et-o—|-oFt'o—| zto—| iz la > on eres 
er'0 ZO OL'0 10°0 20° pa > el eee 
$F'0 OF'0 620 92°0 ze" 6 os Peer eee 
6L°0 LL 0 9L°0 cLo Te ‘0 ‘0 eeeseecereoeres 
ze'1 2e°1 ee'T ze" t 28° + Pees eenvaamnne- 
Gaud 68°S 86'S OT's oT ‘g | eR ea Sa 0 
60°F oF ehh OL ZO" g a. age et a Saal 
82'S 6F'S 8L°S 819 OL" “2 MecR scirtee eee 10 
eye'9 | 969°9 Irs | 969° LE: °8 ie Weak. os 2 T00 
SPL e8°2 Ges 06°8 oo" 6 or (aa T000 
G6F'8 96°8 $56 IT'O1 | 69° “OT Peal Oy T0000 
0S ‘6 ZO°0% {a20:00 | O8s bE] TZ: AT. ce Pa T00000 
Gr OL ZOTL | SkIL | Pe:eT| “09: “LE A Se a eal 1000000 
4 6=8 - 
pas |i rege -|" 38 FZ IS 


TR ahs -(panurquoy) )— ee ATavVil 


. the flexibility that the statis 


many cases than are obtained by the function the writer proposes (the partly veh 
_ bounded function is obviously not flexible enough to represent the earvationsl 
= in the samples given as illustrations), but this added flexibility is scarcely 
a desirable feature when dealing with inadequate samples. : Fee 
% +f Tn functional analysis (as distinct from statistical analysis) the probable Fy 
4g errors of the constants of a curve fitted to a set of data may be computed by 
- the methods developed from the theory of least squares because of the reason-— ret ; 

4 able assumption (in most cases, at any rate) that the variations in the obser- 
-s vations from the true functional relation follow a symmetrical law closely 
_~ represented by the Gaussian normal. In fitting a curve to a statistical 
q series, however, the errors arise in an entirely different manner; they are 
. 


(co 
~~ 


errors in the characteristics of the distribution function itself. These errors 

are seldom small quantities of the first order, and so one is scarcely justified 

in computing the probable errors of the constants of the curve from them, 
even though by means of them the difficulties of correlation emphasized by 

Messrs. Fisher and Kemper are obviated (since statistical characteristics 

are mutually independent). In 1902, Pearson gave™ the probable errors of 
statistical characteristics. ‘The writer uses these values, reducing the order 

of the moments to any desired degree by a suitable recursion formula (Equa- 

tions (22), for instance), and computes two curves, one with the errors 

added to the characteristics as computed from the sample and the other with ae 
the errors subtracted. This gives a range within which the curve character- . 
istie of the population has something better than an even chance of lying. © oe 
The important and difficult subject of errors of sampling has not yet 
received adequate attention from engineers. 

Professor Mavis raises an important question with regard to the repre- q 
sentation of finite and discrete variates by means of a smooth frequency 
curve. Few phenomena in applied mathematics have been so misunderstood 
as the nature of a frequency function. A smooth curve does not imply that 
the variate is continuous or that the variations are infinite in number, 
although in so-called derivations, still abundant in treatises, such motions 
appear implicitly and explicitly. The fact is that the ordinates of a frequency 
curve do not give frequencies. The frequencies are always given by the 
area included between the curve, the class interval, and the two class limits. 
These frequencies, therefore, are always discrete. It is true that the mid- 
ordinate of a class (when the class interval js taken as 1) is usually numer- 
ically close to the frequency of the class, so that this ordinate is often taken 
as the frequency; but the frequency is always given by an area. When the 
variate is continuous and the number of variations infinite, one may only 
estimate the frequency of the values of the variate within the class limits; 
jn the continuous case, the frequency of a single value is always zero. 


? 


nO eRe oe tke 


58 “Qn the Mathematical Theory of the Brrors of Judgment,” Philosophical Transac- 
tions, Royal Soc, (London), A, 198, 1902, pp. 235-299. 
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Tn calling the totally bounded function the most general homograde fune- 
tion, the writer is stating a definition and, consequently, cannot be falling 
into a logical error. He is, in effect, formulating a physical law. — This 
’ formulation is empirical, and it will take experience and not mathematics 
to prove it right or wrong. It is well known that the small errors of precise 
measurements are distributed, to a high order of approximation, according 
Ny to the normal law. The mean and the standard deviation are the two char- 
acteristics necessary to determine. such distributions and, for them, there 
are no other independent characteristics. When the causes producing varia- 
tions, are not infinite in number and are not mutually independent, then 
these two characteristics do not suffice and two others must be introduced. 
The writer then considers the four-parameter generalization of the Gaussian 
normal the most general homograde distribution, more complicated distribu- 
tions being merely superpositions of this type. 

Thiele’s observations, as mentioned by Mr. Fisher, regarding the “pseudo- 
normal” variate, t = h(a), are not discoveries in the theory of statistics, but 
merely special cases of theorems from the general theory of curves, Any 
two curves, provided they satisfy certain conditions of uniformity, may 
always be transformed into one another by suitable transformations. A 
theory of “pseudo-normal” variates, therefore, is quite as general as the 
integral equation or the Gram-Charlier formulations. The great use of such 
an observation (and the writer has availed himself of it) is that tables of the 
probability integral may be used with the more general function. 


The writer was in error in attributing Equation (10) to Mr. Fisher, and 
he did not notice, until it was too late to acknowledge it, that the partly 
bounded function is exactly equivalent to his generic function; he was con- 
fused, on reading Mr. Fisher’s treatise, by the discussion of the “mathe- 
matical zero.” The writer objects to the use of it as a generic function, at 
least when dealing with small samples, for reasons already stated. 

The work of the Danish statisticians is known to the writer only through 
Mr. Fisher’s work and he was not aware of Thiele’s generalization, Equa- 
tion (68). Thiele. however, does not show how to compute the constants of 
this curve in terms of the statistical characteristics; in this respect, at any 
rate, the writer has shown some originality. The awkwardness in the writer’s 
presentation, referred to by Mr. Fisher, was to a large extent, unavoidable. 
In the first place, the writer has definite limitations of expression and knowl- 
edge; in the second place, the “obvious” detail in the analysis and the choice 
of mathematical procedure was deemed desirable in presenting a paper of 
this sort to engineers instead of mathematicians. In defense of the analysis 
of Article 6, Section II, which Mr. Fisher considers. superfluous, the writer 
should like to state that the approach of the logarithmic function to the 


limit, —, as the skewness vanishes is by no means self-evident (although it 
o 


may be well known to mathematicians), If the transformation chosen had 
been the simple “pseudo-normal” variate, ¢ = h(x) = aa*, then for no value 
of the parameter, a, would the resulting curve have reduced to the normal, 


m hecaéenr Svalent Both ropvesstit, any function, hate oh 7 
or not, which is continuous and vanishes, together with its first deviratiges ite 
. infinity in both directions. fod 
ih Mr. Fisher is inaccurate in inferring hint ‘the writer’s stahern abies cree 5 
img the generality of these two equations implies that a Fourier series is 
, _ inferior to a Taylor series. What they do imply is that both the Fourier 
and the Taylor series are mathematical devices which sometimes offer a — 
comparatively easy means of solving a difficult problem — more often than * 
not the solution thus obtained will be purely formal and of little practical 2 
value. aa ce? 
Elsewhere,” the ‘writer has proposed a method for determining ies maxi- — 
mum and minimum floods from a record of yearly floods. In this method — 
the constants have been determined only for use with yearly floods and as 
it stands, it is not applicable to the example furnished by Mr. Hall. The 
form of the equations has been suggested by theory—physical and statistical— 
but their use in their present state can be justified only by the consistency - 
_ of the results obtained by means of them in the problems attempted by the 
writer. The following quantities are used: g = maximum upper deviation 
from the mean, where the mean is taken = 1; ) = maximum lower devia- 
tion from mean; go = maximum observed upper deviation; 6). = maximum 
___ observed lower deviation; A = drainage area, in square miles; M = mean K 
flood, in feet per second; and, N = number of years in record. From 
0.4 go 


0 


- 5 sea ; 
0 \ as flood ‘index. : 


Ba a wae eee es at 


these the following two quantities are computed: = storm index; and, 


aes a 


M 
Using logarithms to the base 10, compute: 


Seb Gee (os ee Ne 1h top “4a: ) me Va he (89) 


from which the maximum upper deviation is computed by the formula: 


Paty 1) (0.74 +038) 10. bollelaps) Beg) 


In cubic feet per second, the maximum flood will be given by M(1 + g). 
For the maximum lower deviation the same two equations are used, except 


Dibra i Jo 
h appears, ———— 18 substituted, and where bo appears —~—— 
Ce EBD @ +m) 


wil Un i ete 
59 Report to the Mississippi Valley Committee, July, 1935. 
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¥ ~The | Stsiy Allen Hae; M. om Soe. C: By Stitt an sieentipil dint | 
- shows the effect of the large flood of 1913 on the yearly flood record of the 
Hudson River, at Mechanicsville, N. Y. If the 23-yr record™ ending in 1912 
is taken, then NV =' 23; A = 4500; M = 41700; g = 0.43; and, b = 0.37. 
Substituting these values in Equations (89) and (90), it is found that 
g = 2.225, from which 3.225 X 41700 = 134600 cu ft per sec for the maxi- 
mum flood. If the 34-yr record” ending in 1923 is taken, then the quantities 
are: N = 34; A = 4500; M = 44000; g = 1.56; and b = 0.41; and thus 
i - g = 2.285 and the maximum flood is 144 600 eu ft per sec. Although the two 
‘graphs show quite different characteristics, the values computed from them ~ 
for the maximum flood are seen ‘to be fairly close. Likewise, from two 
records of the Arkansas River,” the maximum flood is computed to be close 
to 50000 cu ft per sec.” These two rivers have equal drainage areas, but 
quite different flood and storm indices. (The storm index as herein given 
bre is a very variable quantity and characteristic of the particular record only.) 
; To use Table 22: (a) Find the mean, M, and the standard deviation, o; 


es (b) estimate b and g; and, (c) compute A = i and @ = <. Table 22 gives 


am : o 


\ 

4 

values of ¢ corresponding to percentages of time for various values of and @. 1 

Bi ‘ First-difference. interpolations will usually be necessary. Then, compute — 

at: (d) X = (ct + M), which is plotted against the corresponding %-of-time. t 

Table 2 of the paper simply gives the particular case in which g is infinite. © 

3 It may be stated explicitly that the tables are computed for ¢ = 1, but that, | 

os because of the homogeneity of the function, any standard deviation is taken | 

into account in Step (d) of the foregoing instructions. When } is greater 
a than g, the function is left skewed. Because of the symmetry in b and g of 

B. \ the logarithmic linear fractional transformation, the same tables may be — 

used in that case by merely interchanging b and g. 

For instance, taking the 34-yr record for the Hudson River, = — 9,985; 
and, b = 0.525. Hazen gave the value, 0.37, for the standard deviation, Cv, 
in this case, and if an attempt is made to fit the curve with this value, it is 
found to be too far off the main body of the statistics—the 1913 flood pulls 
it up considerably. Even this poor fit gives to this flood a probability of 
being equalled or exceeded but once in about 1000 yr, so that since this flood 
contributes vastly more than its share to the standard deviation, one may 
conclude that a better estimate of this value may be made by taking only the 
forty-two values that exclude it. In this way the standard deviation is found 


™ “Flood Flows,’ by Allen Hazen, Wiley & Sons, New York, 1930. 
& Loc. cit., Fig. 88, p. 89. 

Doc, ctt., ‘Fig. 15, p. 78. 

% Loc, cit., Fig. 29, p. 89, and Fig. 34, p. 84. 
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acl : * 1 the standard deviation, fir t 


Equation (93) is obtained from Pearson’s relations for the error in the | 


a “moments in connection with the recursion formulas, Equations (22), and 


F 4 E i 
TABLE 23.—Compuration oF Prosasitity Curve. yg 


- ¢9) (2) (3) (4) (5) 
“4 a 
: 6.41 7.66 6.86 3.06 0.000001 
A 6.22 7.36 6.63 2195 0.00001 
ae 5.99 6.96 6.34 2.90 00001 
e 5.63 6.42 5.55 2:77 0.001 
5.11 5.66 5.31 2:60 0.01 
ee 4.34 4.66 4.46 2.34 0.1 
- 3.17 3.25 3.20 1.96 1 
4 1.41 1.38 1.40 1.42 10 
6 0.75 0.72 0.74 1.22 20 
as 0.32 0231. 0.32 1.10 30 
“ 0.01 0.00 0.01 1.00 40 
on —0.25 —0'24 —0.25 0.92 50 
—0.48 —0.47 —0.48 0.86 60 
—0.69 —0_67 —0.69 0.79 70 
—0.90 Sire —0.89 0.73 80 
—1.10 —0 197 —1.05 0.69 85 
Sak: —1.09 Sgt 0.67 90 
S31 96 = 2s —1.25 0.63 95 


modified so as to be applicable to the general function. These expressions are 
za not exact, but they give values for the error comparable to those obtained 
by the cumbersome exact formulas. Using the foregoing values for the 
ie Hudson River, the error is found to be 0.04. The curve for ¢ = 0.26 + 0.04 
= 0.3 will be computed. With this value, \ = 1.75 and @ = 4.36. From 
Table 22(d) by interpolating between the columns for , = 1.5 and X= 2.0; 


Upper Limit = 3.285} 
+ =a 


Ratio of Flood to Mean 


ann OOF 


0.01 0.001 0.0001 
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Fig. 9.—YxEaRLY FLOOD Record, 1881-1923, Hupson RIvER aT MECHANICSVILLE, N. Y. 
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ANALYSIS OF CONTINUOUS STRUCTURES BY 


© ‘TRAVERSING THE ELASTIC CURVES 


a 5 “i = oY grikhia Discussion 
By RALPH W. STEWART, M. AM. Soc. C. E. 


Ratpu W. Srewart,* M. Am. Soc. C. E. (by letter).“*—The discussion has | 
os brought up criticisms and questions relating to specific items not satisfac- 
__torily covered in the paper. Tt has also raised the broader question as to 
1 whether the traverse method is a special method of limited scope, or whether 
_ it is a general method which may be used for the analysis of structures for — 
which slope deflection and end-moment distribution are not suitable. To 
answer the latter query it will be necessary first to show its application 
to some cases of single-span beams. 

Figs. 12(a), 12(b), and 12(c) show the simple moment areas and the posi- 


tions of their centers of gravity for all conditions of normal loading; a 
3 Figs. 12(d), 12(e), 12(f), 12(g), and 12(h) show traverse diagrams for beams ; 
_. of constant cross-section having various conditions of end restraint. The b 
__ - double hatching indicates a fixed end; the single hatching a restrained (but Bee 
“not fixed) end; and, absence of hatching, a hinged end. In the following , 
equations, A is the area of the simple diagram and A, the area of a u 


triangular Me -diagram (not shown), with a base extending the full length 


of the beam and an altitude which is the end moment due to restraint, divided 
by # and J. If the beam has a constant cross-section and is of the same 
material throughout, and movements of supports are not involved, both 


_ Y M 
E and I may be omitted, and the moment areas used instead of the ae 


Notn.—The paper by Ralph W. Stewart, M. Am. Soc. C. E.. was published in October, 

1934, Proceedings. Discussion. on the paper has appeared in Proceedings, as follows : 

i December, 1934, by Messrs. Garrett B. Drummond, Austin H. Reeves, HB. G. Paulet, 
A Adolphus Mitchell, and David M. Wilson; March. 1935, hy Messrs. W._H. Kirkbride. 
R. B. Ketchum. A. Floris, and Ivan M. Nelidov: and May. 1935, by Fang-Yin Tsai. 


Assoc. M. Am. Soe. C. E. 
*%Fngr. of Bridge and Structural Design, City of Los Angeles, 


24a Received by the Secretary June 25, 1935. 
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the appurtenant end moment. 
_ For cases in which a traverse 
, fetes a single triangle, it is 
often more convenient to use the - 
angle relationships of the triangle 
for obtaining values of the un- 
knowns than to write a traverse 
equation. In Fig. 12(d), the isos- 
eeles triangle formed by the 
_ traverse gives A = 3 A, and from 
ie Ml Pe wh 


for uniform load- 


ing, or M (the end moment) 


Therefore, M, = 24 (39 Bayh 


In Fig. 12(f), a(2 l \= Agt; 
andieid’ soit eed Ebenalee eee 


= 


==) et Tithe Jeadisves an 
Fig. 12(b), then g = 4 and 

2 ‘ 
A= =~ giving M, = + Pi. : 


In Fig. 12(¢) (assuming counter- 


Fie. 12. . i ~ 
‘ : clockwise rotations to be nega- } 
tive, as in slope deflection), for symmetrical loading from left to right: , 
+ § Ay =P A BAe SOs oo cn cen nie OSD 
and right to left, 
— 05 — 4 Ag = gAy rt Ay Seine. Vantico 6m) 


If the loading is unsymmetrical, use A g instead of 4 A in Equation (66) 
and A(i — g) in Equation (67). Solving these equations for A, and As, and 
substituting the M-values of the A-angles, standard slope deflection equations 


: ge botations 
coeffi ; of J are 


x 


(67) only the 


a8 13 illustrates the traverse of a beam of variable section and ; 
; dissymmetry having fixed ends. All dimensions and properties ‘of this 
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Y Los Ay Vas 
SJ 4, 
4 = 89.15 
19.62 cf 
ot nee | piee F 
ie: Fie. 13. ie 
2 beam have been presented™ by George BE. Large, Assoc. M. Am. Soc. OC. E., $f. 
__ jneluding the computations for the A and A locations and magnitudes. Nae 
g _ From the triangle in the traverse: Ai: = ae % 89.15 = 25.56; M, a 
ee : . it 
— = ase <x 10 = 28.00; A. = ssid x 89.15 = 63.59; and M, = 63.50 <x 10 Rar 
9.12 5.65 PEATE - 


— 51.83. This solution gives the fixed-end moments by a simpler procedure, 
Z involving less computation, than the previously published solution.” Prob- 
eZ lems of the fixed-end: moments in restrained beams cannot be solved by slope 
deflection or by end moment distribution. Their solution as herein demon- 
strated shows that the traverse method is applicable to fundamentals. 

Frequently, in a continuous beam, the direction of rotation at a support 
cannot be forecast by inspection, and the question arises as to how to apply 
This condition is illustrated by Fig. 14, in which the slopes 


the traverse. 


f° %* “The Analysis of Continuo 

s Hardy Cross, M. Am. Soc. C. E., Transactions, Am. Soc. Cc. E., 
also, Bulletin No. 66, Ohio State Univ. 

26 Transactions, Am. Soc. C. B., Vol. 96. (1932), p. 105, Equations (86) and (87), 


us Frames by Distributing Fixed-Hnd Moments,” by 
Vol. 96 (1932), p. 102; 


veh tft for the fll jength of yore ale om 
and the coefficients of the respective A-angles; 


3 an. Tee 


fi OAs Be 2160; and, ct usr pa 26 Ao Ss sy 1620. _ (The si : 
quations follow the rules’ for slope deflection.) _ ‘Solving, Maz 


= 17,61. Continuing the solution it is found that 6, = — 8.80 En 


: Oe = + 19.57 E I. The solution was not affected by the error in the 
diagram. Attention is invited to the fact that a mechanical solution of this 


problem could be obtained by drawing the A-angles to a suitable scale on 


separate slips of tracing paper and manipulating them with their legs inter-— 


secting on the 7-lines until proper closure is obtained. The result would be 
rather rough, but probably more accurate than that obtained by using the 
arbitrary coefficients of various building ordinances. 


P=1000. 
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pa discussion of Mr. Nelidov shows comprehensive ‘understanding of | e 
traverse method and also shows that it can be applied to frames with slo 
_ roof members. In order to do the method justice, however, it will be neces 
Bone to show a shorter and simpler solution of the problem that he pre- 
& sented. Fig. 15 shows (in its deflected position) the same frame and loading 
as Fig. 9, the right half indicating the elastic curvature and the left half, 
a the traverse diagram. The value of A; results from the relative stiffness 
Be factors of the column and the rafter. The value of A, is obtained by summa- 
_ tion of angles from A, upward. ee 
__- Beginning at A, the following traverse equation for the horizontal move- 
ments of A, may now be written: Vo 


ap Rae ot noe a mao 68) 


For the left half of the frame two moment-shear equilibrium equations— 
one for the column and one for the rafter—may be written as follows: 


Wings Sf eS ep ee, SR ae 


‘ and, 
Z My Maj 600 X 10) — 5 Beets ee oe ty IO) 


a in which AH equals the horizontal thrust caused by the load. After trans-— 
forming: M-values to A-values, Equations (68), (69), and (70) are easily 
solved, giving H = 220.1, and moments which agree with those found by 
Mr. Nelidov (Moment-shear equilibrium formulas of the type of Equations 
(69) and (70) are used in the analysis of wind stresses in buildings. The 
principle is that the shear in a column is equal to the sum of its top and 
y bottom moments divided by its height.) 

t No trigonometric functions are necessary in this solution, the number 
and size of equations are less impressive, and the labor is reduced. 

' Several discussers thought that the original paper was confined too much 
to symmetrical structures and loading. The following illustrations will show 
that the traverse method readily solves unsymmetrical cases. Fig. 16 is 
introduced because its solution may be checked conveniently by other methods 
or by handbooks. It will serve to illustrate the method of approach to the 
general case which succeeds it. 

Tf the columns in Fig. 16(a) were on rollers the beam would deflect as if ; 
it were simply supported. The joint rotations would be equal to the end 
slopes of the simple beam, and the spread, d, = di + d, at the bottoms of 

the columns would be equal to the sum of the end slopes times the height 
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ce the frame, which is equal to the area of the = aineram for the beam 


ght of the frame. For equilib- 


(considered as simply supported) times the hei 
H, must be exerted with 


rium, since the column bases do not spread, a force, 


‘ Gite ithe Shade sacar Fie Bao! 
d 
2 x a) + oab= AEM we rktted eae 
oan Eek tl Le Ltettta Bods lade. « soe 
uating x to ds letting k = = the ratio of the stiffness of rei ee 


“Share fos ade 
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= tah Fos and solving for H 1 © 
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si" lech is a andbook formula.” With H known, the moments are sadeenell 
determinate. ye solution takes into account the effect of side-sway. 


Fic. 17. 

es, Fig. 17 shows unsymmetrical loading for a hinged-base, gable roof frame 

of variable moment of inertia. FoNowing the method of approach indicated 
% for Fig. 16, the equations by which this frame may be solved are, as follows: 
a In Fig. 17(a), 
5 dt dy Sd A Ea: SO eae 
+ For a trial H in Fig. 17(b), compute, 
P dy = 2 CAr hy Bs Ae. SO ee (73) 

and, 
True H & Se eittel at ee een a aN 


3 
With H computed, the problem becomes statically determinate. This solu- 
tion also includes the effect of side-sway. 

These applications of the traverse method show that it is not restricted 
in its scope as are slope deflection and end-moment distribution. The latter 
methods require the assistance of some other type of analysis to obtain the 

*“Rahmenformeln,” von A. Kleinlogel, Wilhelm Ernst & Son, Berlin. 
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: embers. _ The traverse method is applicable to the entire oe flexure, 
; ich moments only and straight members only are involved. >) 
ow ere axial shortening of members is involved the traverse can be passed — 
through the altered positions of the joints as shown by Fig. 10, and solved. 
The method, therefore, appears to be entitled to the adjective, “general,” io 
rather than “special.” rit. Be. 
es Mr. Paulet feels that the bending-moment diagram should be available 
-as a guide, so that the traverse can be sketched. Proper observance of laws 
_ governing the angles will make it as easy or easier, to draw the traverse 
independently. 

é Note that in all diagrams each obtuse angle which is the supplement of 
each A-angle presents its opening toward the load, and if there is restraint 
_to resist free rotation of an adjoining joint, there will be a A-angle of oppo- a 
e, site direction between the A-angle and the joint. Using this rule the writer zie Be 
- finds it easier to draw the traverse first, the moment diagram being either ee 
2 omitted or drawn as a by-product of the traverse. It was noted with pleasure 


? 


that Mr. Nelidov omitted the moment diagram from his Fig. 9. 

4 Mr. Paulet gives Equations (15) and (16) to express the series for rela- 
, tive moments and slopes, respectively, in Fig. 6 (c), with the letters exchanged 
_ for the series in Fig. 6(d). The fact that a single simple equation—namely, 
Mans = 4M, — Mn+, in which M = either moment or slope—expresses 
the law of progression for all four series is of more interest. 

With reference to Mr. Paulet’s feeling that for single spans the moment- 
area method is preferable to the traverse method, it is suggested that Fig. 18, 
and the supporting text, be examined. The traverse method is simply an 
amplification of the moment-area method by expressing moment-area relation- 
ships as a geometrical traverse, which clarifies the computations, discloses 
“short cuts” which would not be evident except for the traverse diagram, 
provides a graphical check, and also an easily understood method of treating 
settled supports. 

The statement by several discussers that the traverse method compares 
unfavorably in speed with other methods merits a brief investigation. The. 
solution for Fig. 3 was presented in complete detail to illustrate the prin- — 
ciple of the traverse. A computer familiar with the method would not write 
Equation (6), but would analyze the traverse as follows: In Fig 3 where 
two lines intersects to form a A or A-angle, note that the intersection dis- 
plays both an acute angle (the A-angle or the A-angle) and an obtuse 
angle which is the supplement of the A-angle or the A-angle. Designate angles 
with obtuse openings that face the interior of the frame as negative and 
the angles with obtuse openings that face outward, as positive. The sum 
of the negative angles must-equal the sum of the positive angles. Noting 
that each column traverse forms an isosceles triangle, the computer can add 


fy 


Bote rae 


2 
the angles mentally and write, 3A. = a The subsequent steps can also 


be performed by mental arithmetic. 
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a. ‘ioe will also equal the rotation at each bottom corner. From the stiffness ratios 
i the curvature due to the moment in the bottom of the column will be 2 A; as 
pa shown. The. fixed-end moment of 1000. at the end of the top member gives 
a g-angle of 167, as shown. The solution is now given by the following equa- 


ci tions: . Traversing from the bottom to the top of the column, —3A —4A 

; ‘+ A, = 0; and, therefore, A, = 7A. The upper joint rotation, by adding 

“ angles from the bottom to the top of the column, is equal to 4A; the curvya- 
iy 
a ture due to the end moment in the beam, by stiffness ratios, will beta ; a 
‘% 
and the closure of the triangle under the end of the upper beam gives, 
3 7A = 167, from which A = 26.32, and the bottom and top moments in the 7 


column are, respectively, 105.2 and 868.2. This is a direct solution which ~ 
; appears to the writer to be more expeditious than the previous solution cited _ 
and will compare well with solutions by slope deflection or by progressive 
end-moment distribution. 


’ °8 Proceedings, Am. Soc, C. F., December, 1934, p. 1460. 
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The rotation at the left end of the lower beam is — 
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a fore, if P < 3Q, the rotation is positive, as shown, and if P > 8Q, the rota- 
tion is counter-clockwise. + ice 
_. Continuing the solution it is found that if Q > 5P, the left top corner — 


3 of the frame will rotate contra-clockwise. It is observed that slope-deflection 
+ sign rules may be applied to the traverse method of solving problems involy- 
ing unknown directions of joint rotations. 

a Mr. Mitchell showed an application of the traverse method to Fig. 3, 
assuming that the column bases are 50% fixed. The traverse method, how- 


% 


ever, offers a shorter solution than that given by him, as follows: For the 


4 columns fully fixed at the base, A: = 3 A» For 50% fixation, Ai will become s 
_ equal to $ Az and a rotation will occur at the bottom of the column which A 
will be two-thirds of the change in the value of A, = 4 Ac” This angle of wie 


_ ’ rotation will present its obtuse opening outward the same as A. Equating. 2 
_ the sum of the negative angles around the frame to the sum of the positive — Bae 


fs Pl, agreeing : 


2 
angles, one may write at once 34 A. = a from which M = 
t ) 


with Mr. Mitchell’s result, but doing away with his Equations (17), (18), 

and (19), and his traverse equation. 
The writer cannot agree with Mr. 

“holds no advantage over slope deflection’, 


number of equations required. 
- Professor Tsai indicated in his discussion that he had fully grasped ; 


the traverse method. His equations showing how simple moment areas 
and the positions of their centers of gravity may be derived from (although 
not “given in”) the Ruppel tables’ constitute a useful contribution to the 


29 See properties of 1:2:3 diagrams in “Improved Method of Finding Beam Deflec- 
tions,” by Ralph W. Stewart, M. Am. Soc. C. B., Civil Bngineering, February, ; 
6 Transactions, Am. Soe. C. B., Vol. 90 (1927), pp. 167-187. 


Mitchell that the traverse method 
nor with his findings as to the 


that for this problem it is easier to apply ready made formulas 1 


“method” ea the end moment-distribution method. The traverse solutions for 

hi _ single-span frames, including arched frames as given herein, should certainly 
compete with the column analogy, especially for frames of variable moment 
of inertia to which the Ruppel tables are applicable. . 

The writer thanks Professor Drummond, Mr. Reeves, Mr. Kirkbride, 
Professor R. B. Ketchum, and Mr. Floris for their interest and the last 
three particularly for their encouraging comments. He feels that consider- 
ing the versatility and general applicability of the traverse method, it should © | 


o gain recognition. 
a It is superior to slope deflection in the following respects: (1) It pro-— | 
| 
| 
| 


vides an easy graphical verification of the computed analysis; (2) it uses 
By _ basic principles rather than memorized formulas (“the anesthetics of 
me thought”) ; (3) it does not require the assistance of some other method; and 
; (4) the entire elastic deformation of a structure may be pictured from the 
* -golution by the traverse method. 
x , It is superior to end-moment distribution in Items (8) and (4), and also | 
a in its utility in deriving specific formula which may be recorded for future | 
use. : 


oa % Transactions, Am. Soc. C. E., Vol. 96 (1982), p 145. ; 
. " Bulletin No, 215, Eng. Experiment Station, Univ. of Tlinois, Urbana, Tl. 
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THE SILT PROBLEM 


Discussion 


formed a real service to water-supply engineers by making available for ready 
e reference such a mass of quantitative information on the silting of reservoirs _ 
and the suspended matter carried by streams. His statement of the silt prob- 
lem and his discussion of the origin, transportation, and control of silt 
and of sedimentation processes include a summarization of facts and prin- 
ciples, largely devoid of expressions of opinion, that serves to round out a 
< general picture of the silt problem in a manner both interesting and illuminat- 
a ing. In his discussion of the origin of silt the author touches on geology Res 

but, in the writer’s opinion, gives it too little emphasis. As defined by the : 
author; silt is stream-borne material derived from the disintegration of rocks. 
Two processes, weathering (disintegration in place) and corrasion (the tear- 
ing away or placing in motion of disintegrated material) are involved in the 
origin of silt, and the speed of the first of these processes is materially affected 
by the nature of the country rock. Igneous rocks, in general, are hard, 
dense, and crystalline; and they weather slowly. The older sedimentary rocks 
are likewise relatively dense and offer resistance to weathering. The younger 
sedimentary deposits, particularly recently laid alluvium, are especially sus- 
ceptible to rapid weathering. Consequently, other things being equal, regions 
of igneous and precarboniferous rocks are regions of clear-water streams, 
and regions of recent sedimentaries and valleys filled with alluvial deposits 
are most likely to give rise to streams heavily laden with silt. 

Among the younger sedimentary rocks, differences of texture and cemen- — 
tation are reflected by differences in speed of weathering. Loosely cemented 
sandstones and friable sandy shales are readily susceptible to weathering and 
corrasion and are prolific producers of silt. 

The chemical composition of rock mater 
problem. Rocks that contain high percenta 


tris BS LE 
J. C. Stevens, M. Am. Soc. C. E., was published in October, 
ion on this paper has apeared in Proceedings as follows: 


February, 1935, by. Harry G. Nickle, Jun. Am. Soc. C. B.: March, 1935, by Messrs. BE. W. 
Lane, and Frank H. Bonner, and May, 1935, by Messrs. Morrough P. O’Brien, Harry F. 


Blaney, W. W. Waggoner, and Philip R. R. Bisschop. 
48 Chf., Conservation Branch, U. S. Geological Survey, Washington, D. C. 


88a Received by the Secretary July 29, 1935. 
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4 - Herman Srasrer* M. Amu. Soc. CO. E. (by letter).“*—The author has per- 
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ial is likewise a factor in the silt 
ges of aluminum silicate weather 
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a POE mission ee: ieee and other forms ‘Of sedi piicath rom tl 

noment they strike the earth, is to wear the land surface down to_ the le e 

the sea. This is accomplished mainly by weathering and corrasion, by 
solution, and by water transportation. In some regions of the United States — 

more material is carried by streams in solution than as silt. In other sections - 

. ie opposite is true. By and large the silt load of the streams of the United — 

‘States is several times as great as the load carried in solution, but the latter: f 

is by no means an unimportant factor. 4 

_ ‘The relative speed of weathering and of removal and transportation of the 

, entered material, depending primarily on geology, topography, and climate, 

is a factor in the silt problem. In some regions streams run clear because 

little transportable material is available for water to move and only extraordi- 

‘nary stornis and flouds can muddy the waters. In other regions relatively 
gentle rains are sufficient to place in motion a wealth of weathered material 
awaiting only the means of transportation. 

Vegetation, an incident of climate, is an important factor in the eoahed 
ing of rocks and in the availability of weathered material for water transpor- 
tation. By disruptive action of roots and chemical action of organic solvents 

cat vegetation aids weathering and adds to the volume of material suitable for | 

transportation as silt and to the dissolved load of streams. By protection of | 
H 
: 


bo the land surface from action of wind and water it retards or prevents silty 
2a material from corrasion and thus tends to avoid the overloading of streams 
with suspended matter. Cultivation of field crops artificially aids weathering 
and accelerates the speed of the natural process. Road building, the herding 
of stock, lumbering, and a multitude of other activities of Man that scar the 
earth’s surface also speed up the natural process of weathering. Most of these 
activities are essential to human occupation of the land but some of them 
may be modified so as to minimize their effect on weathering and still retain 
their beneficial character. On the other hand, works of Man, such as build- 
. ings, streets, and paved roads, substantially prevent weathering whereas his 
: river-regulating works and consumptive uses of water tend to decrease the 
power and quantity of water that may serve to transport silt and other rock 
Om débris. After all, Man’s activities, although by no means a negligible factor in 

the silt problem, must be considered as incidental. Wisely guided, Man will 

undertake, when practicable, to diminish the removal of material valuable in 
i situ to places where it will be deleterious, and to direct the deposition of trans- 
ported silt at places where it may serve a useful purpose. In attempting to 
abolish the silt problem Man can be but a Don Quixote tilting at the wind- 
mills of Nature. 

In attacking the silt problem, should the energies of Man be directed  ‘ 
toward prevention or cure? Rather obviously good productive soil in a farm- 
ing region should be maintained in place for agricultural use so far as 
practicable. In this case prevention pays, regardless of the silt problem; but 
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Fig. 10.--SuBDIVISIONS OF THD COLORADO River Basin ABOVE THE PRINCIPAL 
GAUGING STATIONS 


ah use, and silt have stoked the writer for many years. Most of the 

extensive surveys of its stream channels, reservoir sites, and dam sites made 
by the U. S. Geological Survey in the two decades, 1915-1935, were under | 
his general direction. The quality-of-water surveys (determination of daily Z 
silt load and dissolved mineral load of the Colorado and of its more important 


ams. Otherwise, costs are likely to exceed the xalues” conserve ih 
erie! Colorado River occupies one of the major stream basins in thee 
United States and has silt eee perhaps in greater — than any other 


tributaries for a period continuous from 1929, which were made largely by 


©. S. Howard and S. K. Love under the direction of W. D. Collins, of the 
Geological Survey) were instigated by him and he has studied their results” 
with interest. These researches, basic to land and water planning, although 
still inadequate to supply full information, afford probably the best available 
basis for the study of the-silt problems of a large river in the United States. 
A brief consideration of some of their results will serve by example to throw 
light on some of the problems discussed by Mr. Stevens. Fig. 10 is a map 
showing the subdivisions above the principal gauging stations of the Colo- 
rado River Basin and Table 13 shows some of the characteristics of these 


TABLE 13.—Darta Retatep To THE Sint PRoBLEM ON Conorapo River 


Axsove GRAND Canyon, ARIZONA. 
WateER-SHED Supprvisions (Srp Fre. 10) 
Description Se Se a ee ee eS Se ee ee 
A B ae D E F G 
Area, insquare miles............... 40 600 | 24 100 | 24 000 | 31 000 | 19 300 | 50 300 | 139 O 00 
Area, on perenulne: of. totaly a? ieee 29.2 17.3 17.3 22.3 13.9 36.2 Hoo 
Surface Geology, Percentage of Area: 
PTe=LPIASsO cas is. cis « svelte wuwiels 14.3 38.5 10.8 30.8 16.0 24.5 21.6 
Triassic-Jurassic............... 3.0 a3 21.5 37.5 55.3 45.0 22.2 
Postal UTasel O:ts)3, 4)4< 5. sfans loemusteca.s 82.7 54.2 67.7 31.7 28.7 30.5 56.2 
MOUAL oa risks Pee iets pba giere 100.0 |} 100.0 100.0 100.0 100. 
Vegetable Cover in Percentage of Area: : oe gar 
Timber 2G spruce, acceeue lati 23.5 70.2 30.0 26.1 11.0 19.8 30.9 
Woodland (juniper, pifion).. st 646.9 17.0 28.6 36.0 26.3 31.8 24.4 
Brush and: Grass... .cvcccuiesscs 59.6 12.8 41.4 37.9 62.7 48.4 44.7 
US AR oe ee ae 100.0 | 100.0 | 100.0 | 1 
Palatable oe animal-unit, years per a ae gi 40.8 
RCMATE TONG. c teins Se cme ance oats 4.8 6.25 4.4 f 
Stock Aree animal units per . a See fi 4-3 
AC UANO AIL gs irnti is, <isualsicieia site seneate 12.9 16.6 12.0 ; 
Area irrigated, in acres per square mile. 13.7 25.8 7.6 5:3 i3 3 10:8 
cb bela ren square mile ../ 2.1 3.3 2.3 1.9 0.6 1.4 2.1 
ainfa! — in inches. ...... 10.30 14. : : : 
oes 1929-1934, in Percentage of S on en Lae Sith eye 
rea: 
More’ than’ 10 in. 66... kes 46.2 59.4 58.8 51.6 13.2 35.3 4 
BOOMS kv edema teas os: 51.1 | 40.6 | 41.2 | 41.6 | 60.7 | 49.7 47.0 
POBR ERAR Os iNs vases ncdw ab a cea c.. 2.7 0.0 0.0 6.8 26.1 15.0 6.4 
RENCE Miers Cie tsiSly. 6: sin tual 100.0 100.0 100.0 100.0 100.0 100.0 
Run-off, 1929-1934, in inches. ...... 2.00 155 : ; ; aig 
Rtuncof 1929-1934, in percentage of 4.46 1.55 0.21 0.62 0.38 1.78 
POURI rarehncnintete b eisinteien pose bylee 33.2 43.9 . 
Silt, 1929-1934, in tons per square mile.| 696 737 2 440 2530 3 770 2" ao 570 
Silt, 1929-1984, in percentage of total.| 12.1 7.6 19.7 29.5 31.1 60.6 100,0 
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Grand Canyon, with special reference to the silt lo: 
divisions are identified as follows: | oe 


= the basin of Green River above Green River, Utah. 

the basin of the Colorado River above Cisco, Utah. 

the basin of San Juan River above Goodridge, Utah. 

= the basin of the Little Colorado, Paria, and other tributaries 
of the Colorado River between Lee’s Ferry and Grand 
Canyon, Ariz. , 

“ the basin of the Colorado River and tributaries between Lee’s 

ae Ferry, Ariz., and Green River, Cisco, and Goodridge, Utah, 

aS including San Rafael, Fremont, and Escalante Rivers. 

a iM the basin of the Colorado River from Grand Canyon, Ariz., 

i to Green River, Cisco, and Goodridge, Utah. 


the basin of the Colorado River above Grand Canyon, Ariz. 
= In Table 18, the information on area, surface geology, palatable cover, 
run-off, and silt are from records of the Geological Survey; that on rain- 
fall, from the records of the Weather Bureau; that on population, stock popu- 
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. lation, and area irrigated from the Bureau of the Census, and that on vege-_ 


if table cover from the vegetation map of the United States by H. L. Shantz, 
_ of the Bureau of Plant Industry. The 5-yr period considered is from October, 
- 1929, to September, 1934, inclusive. Acknowledgment is made of assistance 
_ by J. O. Miller, C. E. Nordeen, Assoc. M. Am. Soc. ©. E., and Depue Falck, 
3 all of the Conservation Branch of the Geological Survey, in the compilation of 
- the information. 
b Inspection of Table 13 shows. that, for the period considered (which 
includes a year of high rainfall and run-off, a year of low rainfall and run-off, 


- and, as a whole, is not far from an average 5-yr period) about 77% of the 


water and 20% of the silt at Grand Canyon came from the basin above Green 
River and Cisco, Utah, which comprises 46.2% of the entire area of the 
basin above Grand Canyon (see Subdivisions A and B, Table 13); and that 
less than 8% of the water and more than 60% of the silt came from 36% of 
z the basin’s area, situated below Green River, Cisco, and Goodridge, Utah, 
and above Grand Canyon, Ariz. (see Subdivision F’). Properly belonging to 
_ this region of low run-off and heavy silt load is about one-half the area 
related to Subdivision (, or approximately that part of the San Juan Basin 
above Goodridge that is situated in Arizona, New Mexico, and Utah. Includ- 
; ing this part of the San Juan Basin it is probable that what is commonly 
known as the plateau region of the Colorado, comprising an area of 60 000 
~ to 65.000 sq miles, contributes less than 10% of the water and more than 
15% of the silt load recorded at the Grand Canyon gauging station. Clearly, 
if silt prevention on the Colorado is a worthy objective, intensive study 
should be made of the origin, and the possibilities of the prevention of silt 
derived from this plateau region of the basin. Referring to Table 13, it is 
evident that this is, relatively, a region of Triassic and Jurassic rocks— 
loosely cemented sandstones and sandy friable shales—and to these rocks 
must be attributed the origin of most of the silt. 
The plateau region has a low rainfall. One-half the area has an annual 
rainfall of 6 to 10 in. per yr, 15% has a rainfall of less than 6 in.; and 
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the average for the 5-yr period considered is 9.72 in. The distribution of the — 


rainfall is of importance, and the records show that, substantially, it all 


oceurs in one to six storms during the year. These desert storms are tor- - 


rential in character, well adapted to corrasion of the rather pace divided 
weathered material that abounds in the region. 
To the arid climate may be attributed the excess of weathered material 


’ 
awaiting transportation and to the torrential character of the storms, the | 
| 
i | 


heavy: though sporadi¢ flows of silt. 


being’ little amber and that is. confined 6 the eas outer rim of | 


the basin. The palatable cover is light, being confined to brush interspersed _ 

with grass and weeds, occurring as single plants, and the stock and human — 
/ 2 . 

population is low. In small areas there is a grassy vegetational aspect but 


in the main the region has the appearance of a desert with vast expanses — 


of bare rock and sand. Probably 90% of the surface of the region is devoid of 
vegetation. The irrigated area is small and confined almost entirely to the 
head-waters of streams. It is a land of stock ranges that will furnish year- 
round feed supply for less than 225.000 cattle, or the equivalent in other 
stock.. The sparseness of vegetation permits the torrential rains to do their 
work with a minimum of hindrance. 

What measures of silt prevention can be undertaken to advantage in such 
a region? Artifical stimulation of vegetable growth on 60000 sq miles of © 
surface now devoid of cover and more than one-half of which has rainfall 
of less than 10 in., is a stupendous, and a hopeless, task. Could it be accom- 


- plished to an appreciable degree there is no doubt that the results would be 


beneficial. There is no lack of fertility, no lack of seeds in the ground. 
After every soaking rain grass and weeds spring up in abundance, wither, 
and die. Moisture alone is needed, and moisture it is impracticable for Man 
to supply. Already 96% of the scant precipitation is devoted to evaporation 
and transpiration. Only by decreasing the former can much additional mois- 
ture be conserved for beneficial use. Stimulation or vegetable cover by regu- 
lation or grazing operations is worthy of consideration. Total exclusion of 
stock from the region for a time would encourage stronger root growth, 
heavier crowns, and a full opportunity for natural reseeding. This would be 
beneficial, but how much of the 190000000 tons of silt per yr contributed 
by the plateau region would it keep from the river and for how long and 
at what cost? The principal cost would be the annual loss of forage valued at 
about $1500000. To this would be added the administrative expense 
of excluding stock. An annual cost of $2.000000 might be justified if 10% of 
the silt could be prevented from reaching Boulder Reservoir and it were 
worth $150 to $200 per yr to maintain an acre-foot of reservoir capacity. 
Total exclusion of stock would be impossible because of the human relations 
involyed. Over considerable parts of the area, in the Navajo Indian reserva- 
tions, for example, stock raising is practically the entire source of livelihood 
for the people. Nevertheless, the principle of reducing stock population is 
worthy of consideration, and it is reasonably assured that, where possible, 


‘some sligh 


Increased 


os cae , would deplete the water supply. 
‘Since 


ennial streams, a second method of preventing silt, would be to cons 

e water supply of the plateau region so-far as practicable for irrigation or — 3 

ther consumptive use within that region. If the value of the water of the _ 
teau region is $2 per acre-ft for irrigation in the Lower Colorado Basin Beck 

25 cents per acre-ft for the development of power in the canyons, a total = 

value of the order of magnitude of ‘$3 000 000 per yr is indicated. Exclusion — 

of this water from the main river, even without beneficial use in the plateau’ o 

egion on the basis of such an assumed value, would cost at the rate of less | 

_ than 2 cents per ton of silt excluded with it—far less than the cost of main- 

§ taining reservoir capacity by operations at the side under any known method. — 

Reservoir capacity could thus be maintained at a cost of only $20 to $30 
per acre-ft per yr. It seems rather clear that the development of bene- — 

ficial consumptive uses of water in the plateau region should be encouraged | 

‘to the utmost for whatever effects in silt prevention may result. 

P A third method of holding silt in the plateau region would be through the 

“4 construction of detention reservoirs and spreading works. There are many 

- good and some very large reservoir sites in the region. Eventually, many 

of them will be developed. They could serve to retard the flow of silt to the 

3 Colorado although its eventual movement could not be thus prevented. 

: The effective life of such reservoirs would be rather short but the cost of 

many could doubtless be justified by their beneficial effects in conservation 

of water and retardation of silt movement. Theoretically, at least, spreading 9g 
works have merit. If 1 in. of torrential rainfall on 1000 acres could be 2 


- spread to 0.01 in. on 100000 acres, serious silt movement could be pre- 3 
— yented. In the process, local vegetable growth would be encouraged and flow : 
of water to the Colorado prevented. Doubtless there are favorable localities ia 


in the plateau region for works of such character. 
Summarizing, regulation of grazing, development of beneficial consump- 
tive use of water, and construction of reservoirs and spreading works in 
the plateau region all seem to the writer to merit consideration as means 
- of withholding silt from the Colorado River. He would urge that before 
expending largs sums of money in any such undertaking economic benefits 
commensurate with the costs be reasonably well assured, full consideration 
being given to the limitations imposed by the value of storage capacity in 
nd the capitalized cost of maintaining such capacity by 
He ventures the opinion, based on personal 
~ knowledge of the physiography of the plateaus herein considered, that no 
practicable means will be found for holding back permanently from the Colo- 
rado more than a minor percentage of the silt load derived from the sand- 
stones and shales of the region. It is his firm belief that after all reason- 
able preventive measures have been exhausted, a now unknown cure for the 


- reservoirs below a 
operations at the reservoirs. 
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THE SPRINGWELLS FILTRATION PLANT, 
; DETROIT, MICHIGAN . 


Discussion 


By EUGENE A. HARDIN, M. Am. Soc. C. E. 


Eucene A. Harpin,” M. Am. Soc. C. E. (by letter).“*—The favorable atti- 
tade of the discussers leaves the writer with very little excuse for a closure, 
except to express his gratitude. To complete the description of this plant it 
appears that a summary of operating results should be included. Therefore, 


Fic. 19—IntTeRIon Virw OF OPERATING FLOOR, SPRINGWELLS FILTRATION PLANT, 
DETROIT, MICHIGAN. 


the data in Table 5, taken from the March, April and May, 1935, reports of 
W. M. Wallace, Superintendent of Filtration, Department of Water Supply, 


—TThe paper by Hugene A. Hardin, M. Am. Soc. C. E., was published in Novem- 
ber 1934, Proceetings. Discussion on the paper has appeared in Proceedings, as follows: 
January, 1935, by Messrs. F. H. Stephenson, and Robert Spurr Weston. 

10 Asst, Engr.-Designer, Black & Veatch, Kansas City, Mo. 


Joa Received by the Secretary July 2, 1935. 


sphenson’s r ference 
he Fig. 19, a view of the Ries vter shoul 
ever, th t the architectural features were pldeed secondary sharon i 

sch denit n features. The architect developed his design from a general lay 
~~ and plant skeleton already planned and given to him for architectural treat- 
¢ o ment. With this handicap, architecturally, he succeeded commendably Be 
4 obtaining a very pleasing appearance from the straight lines and mass of — 
structures with a minimum of embellishment and with the use of platy, 
ud not unduly expensive, materials. 
BS ‘The writer agrees with Mr. Weston that the Wheeler bottom is a very 
Le desirable type of filter under-drain system. It was one of a number of. types 
“considered for the Springwells Plant, but was not adopted on account of its 
ay considerably greater cost than the perforated pipe under-drain system. 


TABLE 5.—Operatine Dara, Derrorr Firrration Prants: Monrniy AveRAGESs 


i okey Con PLE Sas Tis) Marcn, 1935 Aprin, 1935 May, 1935 


ceo weed ayo) Water- Spring- Water- Spring- Water-— Spring- 
OR i : Works wells Works wells Works wells 
Park Plant | Plant | Park Plant | Plant | Park Plant | Plant 


‘ (a) PLant OPERATION ~ 
‘ie - ‘ Plant Output: | 3 ; 
aS) Lae filtered, in million gallons 
SACRE I AR fer oe ' 128.819 97.105 126.787 97.080 128.129 101.654 } 
\ watee pumped to mains, in mil- 4 
; lion gallons daily........... 125.389 95.687 125.068 95.933 125.953 99.290 
by Loss or use in plant, in million : 
gallons daily............... 3.402 1.418 1.719 2.147 2.176 Pe a | 
‘Percentage of loss or use in Tat | 
a, gov es 1a hoe | Genes ino a 2.66 1.46 1.36 2.15 LZ. 2.32 | 
a _ Wash water, in million gallons 
: BURL ieee conc vee eae ate neces S045 utes 200464 Sos. ae PLGYS> 1, Soares ‘ 
; Percentage of wash water. .... 263." eer mee 22629 | Soe eee LASS, Sh ose ae 
a a 2 | 
nen Filters: 
+3 Filters in service............. 66 29 60 29 61 29 ] 
‘ Filters washed............... 28 25 19 24 18 32 4 
t . Average filter run, in hours.... 38.5 28 60.8 29.1 72.3 22.2 4 
ee rate, in million gallons F *| 
PRECRCEEUBENHT Soy cic cnt eh aN 90 Wie uae 119 137 113 139 ex, 
Rate a head loss, increase in pee 
, feet. per hour.............. 0.14 0.20 0.10 0.21 0.09 0.28 
Nig Gallons per square foot per foot 
P Or head lage! eee in vie ee 826 © 643 1 082 666 1 173 611 
| Chemical Treatment; 
nm ' Fayiccs: a, pounds per million 
PRNOUE > niu deattes wines oc 83.5 61.8 107.1 69.0 78.8 : ~ 
: Alum used, in grains per gallon 0.58 0.43 0.75 0.48 0.55 ted 
a eros puters Heed: 
ulfate, in pounds per day....| 469 211.2 460 213.3 471 228.7 
Ns, in pounds per day. ....| 117 52.8 115 53. 
lpr ie pounds per million : 3 +48 aes 
gallons.......... SSR 0.91 0.54 0.91 0.55 0.9 
f NHzg, in parts per million .... 0.11 0.06 0.11 0.07 0. fi OOF 
Chlorine Used: 
In ity rates in pounds per 
million gallons... i... eke 1.96 1.78 1 . 
In td water, in parts per of 1.0Y 1.86 a0 
PEUSLELOOD s wru: b 34 vA scale gla ev wee fee 0.23 0.21 0.22 
In beg ticle. in pounds per onte are eee 
milion gallons....6....008. 1.06 0.99 1.44 
In filtered water, in parts per ae ab Pie 
MENON. oh Aticaet init 6 0.13 0.12 0.17 0.12 0.15 0.12 - 


wee eee 


375 water 

ac, tap water............ Me 

EiGonfitmed B coli Out. R 900 
Tubes: 


. Sad 
0 cu cm raw water.......... 
10 cu cm a) apuled WAtGE. 6.0. 
; aes WatGr. 0, o-,25 
he a 


Rewied soa al ath) ets See Pee 
Filtered water.............-- 
Plankton eho g 
NU PERS Vee WRECE « biacici Sit mins = pie aime 228 | 333 ‘ , 
a PADDING WALEE soo. isie bane ois ~ er 139 | 122 187 | 125 146] 113 128 
Chemical: 
Chlorine Residual, ap 
Million: Ae 
Applied water 0.10 0.07 0.09 0.10 0.10 | 
Filtered water 0.08 0.03 0.06 0.02 0.04 
Weir water. . 0.19 0.10 0.19 0.08 0.19 
SERIE WHUGE snide tase ee pints = 0.15 0.16 0.12 0.13 0.13 
pH-value, raw water.......... 7.9 8.1 7.9 8.2 8.0 
pH-value, tap water........-. 7.4 7.5 7.4 ELE 7.4 
_ Alkalinity, in Parts 4 Million: ; ' 
COs, raw water. .....-.--.-: 2.4 3.6 2.8 
COs, tap water . (5.6. .sl-s)-> ha FR OO 0 
Total, raw water...... pstsges F 84.4 85.0 84.4 
Total; tap water. ......2%.:.. 77.0 at 78.3 80.1. |. .80.2 
SS 
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ANALYSIS OF MULTIPLE ARCHES — 


Discussion 


me By A. A. EREMIN, Assoc. M. Am. Soc. C. E. 


A. A. Eremin,” Assoc. M. Am. Soo. C. E. (by letter).**—A simplified method 
of computing stresses in a system of multiple-arch spans on elastic piers has 
been developed in this paper. The table of moments and thrusts in Fig. 11 
is a useful guide in the practical application of the method. 

In using the alternative method developed by the author in Fig. 11, time 
may be saved by releasing points starting from the end span of a system. 
This is true especially in reference to a series of variable lengths of arch 
spans. The advantage is twofold: First, the error in computing distribution 
factors for moments and thrusts may be traced easily because the effect of 
one additional span on distribution factors is more evident than in the case 
of a series ‘of spans in which the structure is divided into two or more 
series of spans as suggested by the author; and, second, the step of balancing 
distribution factors in far spans may be omitted because they are negligible. 

Mr. Hrennikoff assumed that joints in a system of multiple-arch spans do 
not move vertically. If the support is compressible, forces at the joint are 
balanced after it has rotated through an angle, a, moved horizontally for a 
distance, A, and vertically for a distance, 3. The angle of rotation, and the 
horizontal and vertical displacements of the joint may be determined from 
the equations of equilibrium of the joint: § H,; = 0; 3 V,;, = 0; 3M, = 0. 
Thus, 


haw + had + hed + he (fixed) = 0.....3...-.-.-.-(86) 
vga + UAA + 090 + vy (fixed) = Oc ......5. 205k, (87) 

and, 
max + md + ms + mz (fixed) = 0.............-..(38) 
in which, hs, ms = the horizontal force and moment factors, respectively, 
corresponding to a vertical displacement of a joint; 0, = vertical force fac- 
puna lente es ats : 


Norn.—The paper by Alexander Hrennikoff, Esq., was published in December, 1 
Proceedings. Discussion on the paper has appeared in Proceedings, as follows Mey. 
1935, by Messrs. L. B. Grinter, N. M. Newmark, T. Y. Lin, A. H. Finlay, and A. W. Fischer. 

25 Assoc. Bridge Designing Ener., Bridge Dept.. Di S i i . 
Ree ee { pt., v. of State Highways, Public Works, 


7a Received by the Secretary July 18. 1935. 
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to a verti cement; and vs = 
at a joint from the given loading. 


Mr. Hrennikof’s method, may be extended so that stresses ean ae 


puted ina system of multiple-arch spans on elastic piers with flexible tie-ro 
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yc 


sit bie: 


Fie, 19 —MULTIPLE-SPAN ARCH ON ELASTIC PIERS WITH TIE-RODS 
AT JOINTS. 


in arches, such as that in Fig. 19. Equations for computing stresses in | 
tie-rods may be written according to Maxwell’s theorem: 


DE, Cn + Aa) + DG + Nene = Nereis hee ee o)) 
DE, Ax + BE (As cL As) + OX, DNS = Aoz igluotieiite raltauisiisitat of eitetaay (40) 


and, 


PXGrN, ee Nee Nowe Nox (Weg 2289) "Agel cet ne eee Ce) 


in which, X,, X., X¥; = stresses in the tie-rods (subscripts indicate number 
of the span); An, Ac, Ac: = changes in lengths of spans at the tie-rods of the 
system (Fig. 19 with tie-rods removed and loaded with the given loading) ; 
An, Aw, Ais = changes in span lengths at the tie-rods when the system of arch 
spans is sustaining a unit force, X1 = 1, acting along the axis of the tie-rod 
in the first span; An, Ass, Ase = changes in span lengths at the tie-rods when 
the system of arch spans is sustaining a unit foree, X, = 1, acting along the 


-axis of the tie-rod in the second span; Au, Aw, As: = changes in span lengths 


at the tie-rods when the system of arch spans is sustaining a unit force, 
X, = 1, acting along the axis of the tie-rod in the third span; and Au, As, As = 
changes in the lengths of tie-rods corresponding to a unit load acting along 


the axis of the tie. 


Then: A, = Lh : hel vam ie == Ly , in which Jn, Le, and 
A, Es A, Es 3 Es 
L, = lengths of tie-rods; A, A», and A, = sectional area of tie-rods; and, 
E; = modulus of elasticity in steel. 
Stresses in a system of multiple-arch spans on elastic piers with removed 
ties (Fig. 19), sustaining external loading, or unit forces acting along the 
axes of the ties, may be determined by Mr. Hrennikoff’s method. The hori- 
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_ ANALYSIS OF THICK ARCH DAMS, INCLUDING 


| ABUTMENT YIELD 


Discussion 


By Messrs. I. M. NELIDOV, AND. A. FLORIS 


I. M. Netmoy,’ Assoc. M. Am. Soc. C. E. (by letter).”*—The developmen ; 


of the method by which the effect of abutment deformations on stresses in 
a circular arch ring can be estimated, is offered in this paper. The author 


has plotted curves for Hy = Eo; mr = me = 8; and Lh = 2 — 12, so that 
bibs ue 


the stresses can be determined directly as a function of + and the central | 


“ 
angle, 2¢, The method considers the deformation of the neutral fiber of 
the arch ring decreased by the effect of Poisson’s ratio and by the introduc- 
tion of a finite thickness of the arch. Furthermore, it is based on the 
assumption that plane sections remain plane after deformation. It ignores 


any possible interference from adjacent arch rings. As is known, the method ~ 


is conventional, being derived from the use of thin railroad arches, in which 
the ratios of cross-sections to the length of the axial line are small. Its 
advantage is in the workable form of the equations. On the other hand, 
any attempt to introduce the effect of the firtite thickness and of Poisson’s 
ratio in the form of an interaction between the radial and tangential stresses 
leads to unusually complicated formulas.’ 

When an important structure is under consideration, a method based 
on the general equations of the theory of elasticity should be used preferably. 
For instance, applying such a method as that developed by C. W. Comstock,” 


M. Am. Soc. C. E., to an arch ring with ! — 0.348 and 2¢, = 60° and 
: r 


er by Philip Cravitz, Jun. Am. Soc. C. E., was published in January, 
1935, Proceedings. This discussion is printed in Proceedings in order that the views 
expressed may be prought before all members for further discussion. 
TSenior Engr. of Hydr. Structure Design, State Dept. of Public Works, 
Calif. 
ta Received by the Secretary July 2, 1935. 
8 “Stresses in Thick Arches of Dams”, by B. F. Jakobsen, Transactions, Am. Soc. 


C. B., Vol. 90 (1927), pp. 500-507. 
“9“On the Stresses in Arch Dams”, by C. W. Comstock, N. Y., 1931, pub. by the author. 


Notre.—The pap 


Sacramento, 


oe 


“TABLE 3 ComPutanions BY THE Teretiinroe OF See) 


Taree is Positive) 


SSS 


Unit STREssES, IN POUNDS PER PERCENTAGE OF DIVERGENCE 
ropaam pais 4 ~ 


Caer - Abutment Crown Abutment 
EEE eee _ Ss 
Ex- In- Ex- In- Ex- In- Ex- 
trade eld trados | trados | trados | trados | trados ie trados 
(a) <= 0.345; and, 2 ¢:=60 DrcrREEs 


—2.75 | 1.99 —5.46 Ey a 0 
1.80 3 38 22545" | 60. ee 66 48 


: (b) = = 0.02; and, 2 ¢: =47 Decrees 


Ce ee | —74.5.) —18.1 | —6.1 | —89.2 0 0 0 0 
Gonacode 56 OL PROT 68.4 | —26.5 | 14.1 | —79'8| —s | 46 ey 


It is to be noted that when the elastic properties of the material of the 


-arch are considered there is a notable increase of tensile stresses in arches 


with large ratios of Bee and small central angles. In the development of the 
r 


results in Table 3 it was assumed that the plane between the adjacent arches 
remains a plane and that no shearing stresses occur within this plane. 

Referring to the derivations of the paper, the writer wished to arrive 
at the author’s results by following a different line of reasoning. Consider 
a circular arch ring of uniform thickness, loaded with the uniform radial 
load. Its main, or statically detérminate, system is shown in Fig. 7(a). 
For this system the moments, thrusts, and shear are known, as follows: 
se Patee: P= fetes and V’ = 0. Furthermore, in accordance with 
Fig. 8: 


Equation (44)” indicates that the funicular frame of the loads and the 
neutral line of the arch are spaced at a distance, e’. 

If the left abutment is fixed, as shown in Fig. 7(b), and the right end 
of the arch is supported with a rigid cantilever extending to the elastic 
center, H, three unknown redundant forces, XY, Y, and Z, are thus introduced, 


which will be determined if the deformations which they should counteract 
are made known. 


10 “Strength of Materials’, by S. Timoshenko, Pt. II, p. 532, 1981 Wdition. 
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(c) (d) 


Fie. 7. 


in which, smz is the unit tangential stress due to bending moment; spi is 
the unit tangential stress due to thrust; sy is the unit tangential stress due 
to shear; and, ¢ is the corresponding total deformation. 

The unit stresses in a cross-section and the corresponding deformations 
are shown in Fig. 8. This stress and deformation distribution is modified 
near the supports by the influence of the deformation of abutments, but to 


an unknown extent. 
The unit tangential stress due to water pressure and the force, Y, will be: 


Keys ee + X cos Gn 4 (46) 
Ue oe try 


Fig. 8. 


F. “The unit tangential stress due to the bending moment produced by the ~ 
water pressure and the force, X, will be: 


maine betel fate (1 +2) |+ X ra (cos 6 — S81) ts (fo — 1) (48) 
pa fi Pp as Oath Flag ad oe SA 2 


and the total deformation due to this stress is expressed by: 


ei) — pete al ry 1 ( r )] Pe Te Tn 
Ace ce eee ee eaten ae Pe Te Tn 
eke, 7, om a) ee 


+ Xr, cos p— Str) ta (Hs — 1) or eee en (49 
( on E.In Ty i : 
In accordance with Fig. 7(b), the total stresses in a section are: 

: Me Oo ey 2 ny ee ee (50a). 


Poss PE A Oe Ot - ots aw cs haw oes ee 
and, 


Y= .On+ Xai ott winds weak os. vanes che OUR 


It is noted that as with @¢ = 0, and V = 0, the force, Y = 0, and does 
not enter in Equations (50). 


The relative displacement of the End EF C of the rigid bracket directed 
along the X-axis will be in accordance with Fig. 7(c) and Fig. 7(d): 


Ax = Ax, + Ax = 2Aycos gd: — 2Aysin g, — 2p € et: — + 08 ds ) 
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| and cu 


in which My = — pe ree +Xrn (cos $d; — ut) + Z. Due to the small : 
F 1 
influence of A’, c = c will be assumed in future operations. 
“2 After substitution of Equations (50) and (53) into Equations (51) and 
a (52) and subsequent integration, the following expressions are obtained for 
4 evaluating the unknowns, XY and Ze * 
4 Tn Er Tr Sin dy ae 
Pe Te 4 — Sin Tp + cos os — | nin g — pe cos ¢ | oe 
4 a ia ‘Gare eee ) a 
: 3 d sin? d, or ote 
=e Np a ae beets ($4 : A sin? ¢— mes) 4 ae 4+ — + sin 24) re, 
ee E, \ 2 oe, dy ath 4 A ey 
fag pe = (2-4 sin2 6) +54 nsin 6,| 2 TSS Z (£4 fen 6] " 
Go beeN 2 4 t 6, t t 
Z tT, {sin oy: %, SIN gy r 
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gx 
«aera En ret ian. oe 
ed against Equations (18) to (25) of the pap 


: Gia eae 
uations (58) may be chee 
by noting that Equation (52) is ¥ = ¢: mag Z, and by comparing Equa- 
d * - b- ctn a 


on (50a) with Equation (3) of the paper,' from which it Callous tha 


I 


E or sin , Pons misisie oe 
= ( for —r,). After substituting these expressions into Eque 


k. tions (55), Equations (18) to (25) of the paper are obtained. oF 
ais A question of the relative importance of Poisson’s ratio and of the 
deformation of abutments arises. Table 4 shows a comparison of stresses, in 


{ 


TABLE 4.—Comparison or STRESSES 


(Tension is Positive) 


; 26 4 Unit Stresses Parcenrace or DiveRGeNce 
ly i ety es See a a ee ee eee : 
r ad me=mr| Es t At Crown At Abutment At Crown At Abutment 
86 8s 8e a 8e St 8e 86 8 5 55 
0.01} 150 oo 0 —44.2 | —43.0 | 41.9 | —45.3 | 0.005 0 Ok. '0.4e Orheratnes 
0.01| 150 5 0 —44'2 | 43.0 | 41.9 | —45.3 | 0.005 0 Oh ono 0 
0.01| 150 co | E* | i | —44.2 | —43.0 | 41.9 | —45.3 | 0.005 0 Oliees Oh aan) 0 
0.01} 150 5 E.* | 1 | —44.2 | —43.0 | 41.9 | —45.3 | 0.005 0 O17 0-0 0 
0.01| 150 o | Eo* | 30 | —44.2 | —43.0 | —41.9 | —45.3 | 0.005 0 oO}. Of 0 0 
0.01| 150 5 E.* | 30 | —44.2 | —43.0 | 41.9 | —45.3 | 0.005 0 (Na iaeatnelezt) 0 
0.15| 100 co 0 —4,20 | —1.07 | +0.48 | —6.78 | 0.28 0 OaleOslites) 0 
0.15| 100 5 0 —420 | —1.07 | +0.42 | —6.76 | 0.27 0 0 |—1951 200 0 
0.15| 100 coe oe ee | 4 93 | 1-02, |) —0.11\.] —6 216, |» 0:24 ON ees eet Me Onl eata 
0.15| 100 5 E+ | 1 | —4.22 | —1.14 | 0.14 | 6.17] 0.24 On Sez al tao te 
0.15| 100 x | B* | 30| —4.45 | 0.78 | +0.20 | —6.46| 0.27 +6 | —27 |-58 | —4 0 
0.15) 100 5 E.* | 30 | 4.43 | —0.81 | +0.16 | 6.42 | 0.27 +6 | —24 |—66 | —5 0 
2 0.50} 50 oo o | .. | —0.42 |'+0.17 | +0.36 | —1.13 | 0.38 0 0 |) 0! |.0 0 
> 0.50} 50 5 0 | | —0.51.) 40.19 | 40/26 | —1:08 | 0.37 | +20) +-12 1-28 2 
0.50} 50 o> | B® | i | —0.78 | +0.67 | +0.00 | 0.67] 0.39 | +85 | +286] f+ [41] +1 
me 0.50] 50 5 B* | 1|—0.78| +0.62| 0.00| —0.68| 0.38 | +85 | +259] it |—40 0 
——-«5.50| 50 >» | B* | 30 | —0.86 | +0.84 | 0.04 | 0.54 | 0.40 | +104 | +386) jf |—-52| +3 
; 0.50| 50 5 B.* | 30 | —0.85 | +0.80 | 0.05 | 0.54] 0.40 | +101 | +363 | + |—52 | +2 
a | ee a Ee ——— 
a * Finite. + Reversal of sign. = i 
, a 
pounds per square inch, for pe = 1 lb and for arch rings varying from a 


very thin arch, with a very large central angle, to a very thick arch with 
a very small central angle. Table 4 also gives the percentage of divergence 
of stresses from those computed by a conventional method; that is, with 


small Bt ; with m, = 00; and with rigid abutments. a 
¢ 
The unit stresses in Table 4 were computed by the combined expressions, 
Equations (46) and (48). The bending stress due to water pressure expressed 


by the first part of Equation (48) is numerically equal to +0.20 lb per sq in. 


3 for the arch, with Pegs 0.15, and to +£0.07 Ib per sq in. for the arch, 
fc r 


Peis pene = 0.50. ‘This miecat is OF minor ins ik Bt at. the Aes s , but 
Ba ~ bain 
Beat! ‘of. significance at the spidoete if considered. } ‘ 
_ _‘The data in Table 4 indicate that i in an arch ring of an average thickness * 


Se aa with an average central angle the divergence of stresses is about 

ves 10%, and only the deformation of abutments causes a greater divergence. — 
ae For a thick arch ring with small central angle, the divergence due to 
cs", ) Poisson’s ratio is about 20%, whereas that due to the deformation of abut- 
a 5 ments is increased many times this aaah The divergence of shearing 
“a stresses is negligible. ~~ 
. uh) The limits of the ratio, te = 1 — 30, were assumed, in order to visual- 
#4 t . 
, 


. ize the effect of the extremely large deformations, although actually this 
ae ratio will never be greater than 1. It was originally applied by Fredrik 
Vogt, Assoc. M. Am. Soc. OC. E., for the case of an abutment of constant 
thickness, ¢, and with the forces, M, P, and V, uniformly distributed along 
‘the height, h. In the case of an arch dam this condition is not fulfilled; 
both the thickness, t, and the acting forces vary along the height, A, pro- 
ducing a warped surface of an original abutment plane (which is ordinarily 
a warped surface from the start, due to the requirements of the excavation). 
The problem becomes complex and a practical issue for its solution must 
be found. If the arch rings of a unit height, as well as the areas near the 
abutment were separate one from another, then under the forces acting they __ 
would displace in relation to each other, as the keys of a piano, each one an 
independent element. The relative elasticity and continuity of the founda- 
tion makes the deformation spread to the neighboring units over a certain 
distance. From a list of seventy-four arch dams in California the writer 


found that the ratio, “ = = , at the crest varied from 0.046 to 1.000 


id ad 


(average 0.200) and at the base from 0.009 to 1.000 (average 0.042). If the 
entire height of the dam was assumed as h, the same ratios would be: At 
the crest between 2 and 45 (average, 16); and, at the base, between 1.4 and 19 


(average, 5). As the deformation coefficients, ¢, y, and p, increase with ue 


the selection of larger values of ie will be on the side of safety; but even 
t 


with an increase of five to ten times the ratios referred to, h will be only 
about equal to 1. The ratios in which the entire height of the dam is con- 
sidered do not seem to be justified. 


A. Frorts, Esq.” (by letter)."*—The title of this practical paper and the 
introductory remarks made therein, are rather misleading, because the plotting 
of curves is not an analysis or graphical solution of a problem. The author 
utilizes a known theory by arranging its results in the form of diagrams for 


Dipl. Ing., Los Angeles, Calif. 
lla Received by the Secretary July 22, 1935. 


eae Tis t be heii tn countries in ‘which the. cee 8) . 

e possibl e, however, arrange, the diagrams i in such a - 
make them Sika ‘of this 2 Ree N. Kelen, in his iat D 
0k on arch and multiple-arch dams, for instance, gives erapha the, ve 
ae independently of the pen of units.” . 


Jog te Atta addition. 


Lhe author is to be Saha ee for his 2 eS efforts if fait iy 
- the work of engineers rs engaged i in the design of arch dams. 


; 2 “Die Staumauern,” von N. Kelen, Berlin, 1926. 
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hagsts HYDRAULIC JUMP IN TERMS OF 
DYNAMIC. SIMILARITY 


Discussion 


By Messrs. NOLAN PAGE, ANDREI I. IVANCHENKO, AND 
F. T. MAVIS AND ANDREAS LUKSCH 


* Nowan Pacs,” Jun. Am. Soc. C. E. (by letter).“*—A research problem 
with which the writer recently had contact, involved a study similar in 
general to that presented by the authors, but resulted in values for length 
of jump sufficiently different to be of interest. This problem is one of several 
which have been and are under study by the Corps of Engineers, War Depart- 
ment, at the Hydraulic Laboratory, University of Iowa, Iowa.City, Iowa, in 


connection with construction activities in the Upper Mississippi River and 
Ohio River Divisions. 
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Fig. 11.—Mopen or Srinuine-BAsIn DESIGN. 


To Sump 


As a step in the development of criteria for design for stilling-basins to 
protect navigation dams being built in the Upper Mississippi River, the 
simplified model shown in Fig. 11 was constructed and tested through a range 
a a ge Ls 


Norn.—The paper by Boris A. Bakhmeteff, M. Am. Soc. C. E., and Arth 
Matzke, Jun. Am. Soc. C. H., was published in February, 1935, pruners Dineen 
on this pap per has appeared in Proceedings, as follows: March, 1935, by Hunter Rouse, Bsa. ; 


May, l! by Messrs. Sherman W. Woodward, Robe K 
niieuen le A tele Robert EH. Kennedy, L. Standish Hall, and 


; and August, 1935, b Mes ry: 
aenreSua Ea ae y srs. F. V. A. B. Engels, Baldwin W. Woods, 


“ Asst. Engr., U. S. Engr. Sub-Office, Iowa City, Iowa. 
“co Received by the Secretary July 31, 1935. 


x ren ih a 


h each supporting pier, built one-twentieth the size 


; tes for Mississippi River Dam No. 20, and set up in a glass-sided { nes i 


-TeT- 


TABLE 3.—Test Data 


Measurep Deptus | Meas- Kinetic 
or FLtow, IN Frnt | ured 
PE Se SS Hength “of “Ratio, 

jump, L, L 


2 f 2 . ey ol ae Fe et | ap Ry iy 

oe iste 4h » de } in feet da 5 dy an €1 as €1 & 

e8) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
pile 0.404 1.004 0.031 0.295 1.6 5.42 62.8 0.0381 | 0.294 | 0.203 
13 0.800 1.004 0.058 0.396 2. 5.30 32.6 0.058 | 0.394 | 0.398 
14 0.800 1.004 0.072 0.406 2.1 5.18 25.9 0.072 | 0.404 | 0.398 
16 0.800 0.764 0.066 0.364 2.05 5.63 21.2 0.087 | 0.476 | 0.606 
17 0.800 0.763 0.072 0.376 pe 5.58 19.2 0.094 | 0.493 | 0.607 
19 1.205 1.002 0.084 0.472 2.95 6.25 21.9 0.084 | 0.471 | 0.606 
+20 1.207 1.001 0.090 0.484 mE 5.17 20.2 0.090 | 0.483 | 0.606 

22 1.980 1.000 {| 0.140 0.584 3.1 5.31 12.3 0.140 | 0.584 | 1. 

23 1.988 1.002 0.153. 0.599 3.0 5.01 Ad 0.153 | 0.597 | 1.000 
25 2.000 0.765 0.180 0.529 Bod 3.97 6.50 0.235 | 0.691 | 1.510 
26 1.992 0.768 0.213 0.565 2.5 4.42 5.21 0.278 | 0.736 | 1.498 
28 2.990 1.004 0.227 0.682 2.65 3.89 6.85 0.226 | 0.680 | 1.502 
29 2.996 1.006 0.256 0.719 3.0 4.17 5.86 0.254 | 0.715 | 1.498. 
31 2.010 0.763 0.187 0.549 1.8 3.28 6.17 0.245 | 0.720 | 1.523 
32 2.985 0.827 0.273 0.625 ia 2.72 4.06 | 0.330 |! 0.755 | 2.005 
33 2.985 0.826 0.274 0.632 1.8 2.85 4.03 0.331 | 0.765 | 2.006 
34 2.990 0.827 0.274 0.642 1.85 2.88 4.03 0.331 | 0.776 | 2.008 
36 4.00 - 1.002 0.341 0.778 2.5 3.21 3.88 0.340 | 0.775 | 2.010 
38 5.00 1.162 0.446 0.931 2.4 2.58 3.21 0.384 | 0.800 | 2.015 
40 5.01 1.001 0.463 0.836 1.75 2.09 2.32 0.462 | 0.835 | 2.520 
41 5.00 0.999 70.453 ED iis Sie, oi oars oc 2.41 0,454 | 0.801 | 2.525 
46 4.98 1.163 0.389 0.906 3.1 3.42 3.98 0.334 | 0.778 | 2.010 
49 6.00 1.312 0.422 1.000 3.2 3.20 4.21 0.322 | 0.761 | 2.020 


Instead of the designed sill and stilling-basin, a smooth, level, wooden floor 


was installed, as shown in Fig. 11. Water was supplied to the model through 


a 10-in. line from a constant-head tank two floors above the flume, passed 
through the apparatus as indicated, and discharged into a sump where 
it was re-pumped to the constant-head tank. Model quantities were deter- 
mined by means of a rectangular, suppressed weir. Depths up stream from 
the model gate and down stream from the hydraulic jump were measured by 
vernier hook-gauges in stilling cans connected to the piezometers indicated 
in Fig. 11. Depths below the Tainter gate, up stream from the hydraulic 
jump, were measured with a vernier point gauge. 

Table 3 contains only those data taken in tests on this model which 
are pertinent to this discussion. The lengths of jump tabulated were 
measured directly: in the model, taking the distance from the beginning of 
the jump, which was always sharply defined, to the point of highest rise 
down stream from the jump as observed through the glass sides of the flume. 
The accuracy of measurements has been indicated in the tabulation. It will 
be noted that the values, «, in Table 3 are not computed from d,, but were 
taken as equal to do, the height of pond-water surface above the temporary 
floor of the flume. This was considered more accurate than computing & 
because of the magnified effect of slight inaccuracies in measuring d. 
Theoretically, do should have been corrected for velocity of approach and 


ve ee! ae 


rere ee 


< eects widths at Section 1 and Section 2 known, the lateral expansi 


Values of di and di 


Ree 1c s iy on : € ra de Mv eS ; 
ar atus such as this model, in which the hydraulic jum 
ally as well as vertically, measured values of d’, could not b 


follow the curve for a straight-sided channel given in Fig. 3. Howe 


d be taken into account in theoretical equations for height of jump. 
nis was done for the model and resulted in the curve marked “d’s, computed k 
mode ” in Fig. 12. Except in the range near the maximum height a 


y 
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2° 


ee oO + 
ab Rs Hoa's eae nO 
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Values of —> 


Fie. 12, 


jump, measured values fall very close to the computed curve. As was the 
case in the authors’ tests, the range in which the departure of measured from 
computed values occurs was one of unstable flow. In some tests a standing 
wave condition alternated with a hydraulic jump, the alternations following { 
an approximately equal period for a given model setting. 

Lengths of jump measured in the model have been expressed as ratios of 8 
the depths down stream from the jump, and are plotted in Fig. 12. The 
authors’ curve expressing the relation of this ratio to }, has been transferred 
from Fig. 6 to Fig. 12 in order to compare directly the results obtained 
in the model with those obtained by the authors, It will be noted that, 
although the two curves assume a uniform slope at about A, = 14, the authors’ 
curve is below the curve drawn through model test points and slopes down- 
ward as A; increases, whereas the model curve in this range was drawn with 


zero slope at a value of = = 5.60. It is thought that the jump may have 
2 

been shortened in the authors’ apparatus by friction on the sides of the nar- 

row channel. 


Sa -_ 


g ye elevation corresponding to the quantity. When using the curvy 
of Fig. 12, it is necessary first, to assume a value for d,; from that, compute 
_ an elevation for the stilling-basin floor; and, then, compute by means of ds, an 


_ of d, are not directly indicated by computed values of tail-water elevations. __ 
, A curve that is more directly useful as a means of computing gtillnne eS 
_ basin floor elevations is that shown in Fig. 13. This curve was first called 
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to the writer’s attention by the late Floyd A. Nagler, M. Am. Soc. C. E., 
and was used in the design of Mississippi River Dam No. 2.“ The abscissa te 
in this case involves only the quantity per unit length of stilling-basin, which 
is readily determined from known river quantity and length of spillway, and 
the energy head up stream from the hydraulic jump, which may be taken as 
equal to the height of the pond-water surface above the stilling-basin floor. 
In using this curve, a stilling-basin elevation is assumed to provide data from 
which a basin elevation is computed, this process being continued until both 
elevations agree. Although trial and error methods must be used with this 
curve as with that of Fig. 12, computed and assumed values are directly 
related and converge rapidly. The basis for the abscissa of Fig. 13 is the 


=f equation for energy head at Section 1, 


« “Taboratory Tests on Hydraulic Models of the Hastings Dam,” by Martin B. 
Nelson, Bulletin 2, Univ. of Iowa Studies in Eng. 


; pene pecayse of this. If i caeaia curves ae used, and they are 


usually sufficiently accurate for design Peryones, either Fig. 12 or Fig. 13 


; mall’ give identical results. 


- Acknowledgment is made to the District Engineer, U. S. Engineer Office, ; 
- §t. Paul, Minn., for permission to use the data given in this discussion. The 


. - tests were run by Frederick S. Witzigman, Jun. Am. Soc. C. E P 


 Awprer I. Ivancuenko,“ M. Am. Soo. C. E. (by letter). —The results 
reported in this paper fully confirm the usefulness of the general methods of 


analysis applied. The authors have presented the first comprehensive explana- 
tion of the hydraulic jump, a phenomenon which was once considered quite 


puzzling. 

Little can be added to the clear graphical representation of the relation- 
ships between various elements of the jump, and particularly the longitudinal 
characteristics illustrated in Figs. 6 and 9. The writer is particularly inter- 


ested in the curve of Z plotted against } in Fig. 9. This curve was con- 
a 
sidered by the authors as the basis of their solution of the longitudinal 


elements of the jump. The writer proposes the following exrressions for the 
equation of this curve, obtained by the customary method: 


= = 108 A haa ets 


2 
or, solving for L, 


D = 10.0 dp Oras. Sat babies, eee 


in which dj = d, — dh, or, 


Le 10.0 Ch SE ee eee 


Equation (52) is quite simple,. being derived from experiment; and, 
assuming that it is based on methods of approach that are physically correct, 
it is more justified than any other known formula for the length of the jump. 
In fact, the length of the jump is doubtless the function of the height 


of the jump, (he — hi), and the state of flow, which is represented by A, 
the kinetic flow factor. 


“Civ. Hngr.; Assoc. Prof., Hydraulics and Wate 
Novocherkassk, Union of Socialistic Soviet Rennie nS Oman. dncusteal tae, 


“4a Received by the Secretary August 7, 1935. 
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values of 2 taken from Table 1. The results are given in Tabli 


; TABLE 4.—Cueck Computations or THE Ratio, ie 
— — : J pat 2 F 
at ee rch a , 
» ‘ Values aR Values + 
a — etic og = Kinetic Hai: =a 
I ow y Experi- , i- pages 
. | factor, ment: — ae xe fatter, eae be a i = 
ee A. || By Equa-| curve, — > A |—————_] By Equa-|'curve, as 
4 tion (50) | _ a tiom(50) oe 
oa From | To a a From | To ; HL, in 
; ig. 6 
@!| @ |@ | @ | © © |@m!| @ | @ | @ (6) (6) 
721 | 3.76 | 7.47 | 8.77 | 8.297 8.23 36 | 16.79 | 6.21-| 6.63 | 6.290 0 i: 
30 | 3.94 | 7.55 | 8.79 | 8.225 8.12 1s | 20.80 | 5.83 | 6.32 | 6.046 8.05 
40 | 4.51 | 7.06 |.8.32 | 8.022 7.87 6 | 24.80 | 5.40 | 5.90 | (5.852 5.85 
-43| 4.52 | 7.67| 8.91 | 8.019 7.80 17 | 29.66 | 5.43 | 5.99 | 5.662 3.707 al 
25| 5.35 | 6°51 | 7.65 | 7.773 7.65 39 | 30.53 | 5.32 | 5.67 | 5.630 5.63 
45 | 5.47 | 6.62 | 7.55 | 7.740 7.55 35 | 44.81 | 5.11 | 5.53 | 5.246 5.25 
41 | 6.63 | 7.13 | 7.79 | 7.470 7.35 37 | 53.89 | 4.98] 5.42 | 5.188 510 
24| 6.72 | 7.06 | 7.83 | 7.452 7.30 32 | 62.25 | 4.67 | 5.14] 4.940 4.95 
98.| 8.55 | 6.63 | 7.10 | 7.127 7.13 34 | 68.70 | 4.81 | 5.06 | 4.847 4.85 
261 9.90 | 6.63] 7.03 | 6.936 6.94 33 | 74.56 | 4.65 | 5.17 | 4.774 4.80 
29 | 11.87 | 6.73 | 7.13 | 6.710 6.70 38 | 78.69 | 4.24 | 5.10 | 4.727 4.72 
Beir Petia nny leche best att owe pitt d cuiiieil 1) Nunoiilog le eaeoee, 


The writer is fully aware that the data in Column (6), Table 4, indicate 


certain discrepancies. However, obviously, as the authors state, the tracing 


of the curves was carried out so “that they must pass within the plotted 
ranges, indicating the possible range of the L-values.” For all the curves of 


Fig. 6 the range of L-values from Runs Nos. 25 and 45, as well as for Run 


No. 28, lies below the curves. This may indicate some degree of unreliability 
in the experimental data for these runs, because the proposed formula, 
Equation (50), agrees very well with experimental data in all other cases. 


F. T. Mavis, Assoc. M. Amu. Soc. O. E., and Anpreas Luxsou,” Esq. (by 
letter).“““—The authors’ observations of the length of the hydraulic jump in 
a narrow channel are a welcome addition to the literature of hydraulic 
research. The transfer of data from the model tests to the prototype should 
be made with caution, however, in the light of the authors’ statement (follow- 
ing Equation (15)) that “in analyzing open-flow cases * * * the principle 
[of similitude] in general may be applied to geometrically similar cases only.” 

Referring to Equations (16) and (17 ) presenting the theory of the hydrau- 
lie jump the authors state that these equations “apply equally to a jump. at 
the foot of the Boulder Dam and a small-scale model in a laboratory flume.” 
Tt should not be overlooked, however, that in the derivation of these equa- 
tions the effects of channel friction are wholly neglected. It is not unreason- 


Iowa Inst. of Hydr. Research, and Act- 


a » Di in Charge of Laboratory, 
en oes ok ee The State Univ. of Iowa, Iowa City, 


ing Head of Dept. of Mechanics and Hydraulics, 
Iowa. 

4 Research Asst. Engr., 
City, Iowa. 

46a Received by the Secretary August 26, 1935. 
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‘ire it seems reasonable uit the tow of jump in a narrow oo 
ould be less than in a channel relatively much wider, even if the side walls 4 
both channels were made of glass. Due to the effects of the side walls the 
| ‘ 
R ngth of the hydraulic jump in a narrow channel may be relatively less — 
than in a wide channel. Further tests are needed to show what effect the — 
relative width of a channel may have upon the length of the hydraulic jump > 
3 under otherwise similar conditions. 
: On the basis of tests conducted for the Miami Conservancy District, Riegel 
and Beebe stated that, 
“The length of the jump is approximately five times its height. While the 
beginning of the jump is fairly definite, its lower end is indefinite, and this 
figure represents merely a general estimate. The lower end was taken as the 


‘place where the water surface became and remained sensibly level, a place 
which was variable i in position and difficult to locate.” ; 


In 1927, Safranez® studied the hydraulic jump experimentally at the 
a Technische Hochschule in Berlin, Germany. The tests were conducted in a 
glass channel 19.6 in. wide. Lengths of jump were reported” for eighteen 


tests in which the kinetic flow factor, 4, = “4 varied from 3 to 365, 
@ g % 

lengths of jump from 1.0 to 3.3 ft, and heights of jump from 0.18 to 0.76 ft. 

Safranez proposed a tentative formula for the length of the jump which, using 

the authors’ symbols, reduces to: 
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An analysis of these data“ for kinetic flow factors less than 80 (the range 
¥ reported by the authors in Table 1), showed that the simple formula: 


Lj. = 32 dp. wunteba-~Raet eee. ee 


2 was in better agreement with the observations than Equation (53). The 


‘ L : 
observed ratios, —4, for the nine tests ranged from 4.6 to 5.5 and the mean 
4] 
absolute error in calculated lengths of juhp, using Equation (54), was 6 per 
cent. For kinetic flow factors greater than 80 the length of the jump was 
4.5 times its height, ranging from 4.0 to 5.1 for nine tests, with a mean abso- 
lute error of 5 per cent. The question naturally arises whether this general 


‘ 
‘"“The Hydraulic Jump as a Means of vee as Bnergy,” by Ross M. Riegel and 


arty Beebe, Technical Rept., Pt. III, Section XVI, p. 85, Miami Conservancy District 


48 “Untersuchungen ueber den Wech 7 r ; 
v. 10, Heft 87, p. 38 (1929). selsprung,” von Kurt Safranez, Der Bawingenteur 


* Loc. cit., Table 4. 
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cat lar flume, 26 in. wide, lined with galvanized sheet metal. Water flo 
r a sluice-gate with a rounded lower lip upon a level apron of steel plate 


and s smooth cement mortar. The rates of flow varied from 0.48 to 1.44 cu ft ; 
per sec per ft width of channel, the depths of flow above the jump, from 0. 05 a 
to 0.15 ft, and the heights of jump, from 0.22 to 0.73 ft. In these tests the nee 
posses of jump varied from 1.4 to 4.5 ft. Fig. 14(a) summarizing the data 
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4 Fic. 14—SumMMARY OF OBSERVATIONS, LENGTH OF HYDRAULIC JUMP. os 
shows as ordinates the ratio of length to height of jump and as abscissas, the 
; kinetic flow factor. Fig. 14(b) shows as ordinates the ratio of length of roller 
to height of jump. An analysis of the tests conducted in a flume, 26 in. wide, 

- indicated that the length of jump was approximately equal to 6.2 times 3 

its height and that the length of roller was approximately equal to 5.2 ia 


_ times its height. 


2 V. J. Aravin® has presented a semi-rational analysis for determining the ; 
length of the hydraulic jump in which he followed the suggestion of Professor 

_ M.A. Velikanov “to divide the stream, when computing the losses, into two 4 
= regions: dead zones, where only rotatory motion of the liquid is taking place, 3 


and zones where progressive motion is observed.” Empirical coefficients were 
d from tests by Bakhmeteff, Pietrkowski, Safranez, Einwachter, and 


obtaine 
Aravin. The results obtained by the simplified formula suggested by Aravin, 
are: 
— 3 £2 
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BOE inatio f the Length of the Hydraulic Jump,” by V. J. Arayin, 
eae Belen fic Research Inst. of Hydrotechnics, Vol. 15, pp. 4-57 (1935), 


Transactions, Scienti 
(Translated ‘trom the Russian by Andreas Luksch.) 
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_ STABILIZING CONSTRUCTED MASONRY DAMSE ee 
BY "MEANS OP CEMENT INJECTIONS’ "3 


Discussion aan 


author’s clear, complete, and interesting account of methods, difficulties, and 
results, make this paper unique and, therefore, worthy of something more = = 
ig than casual and perfunctory discussion. a 
‘i In order fully to appreciate the situation which the author’s work was 
intended to, and undoubtedly did, improve, and how such a condition of 
affairs could come about, it is necessary to cast a glance backward over the 
development of hydraulic construction in India. 

Masonry dams exceeding 100 ft in height, and about a mile long, were 
undertaken by the Public Works Department prior to 1870 and pushed to 
é completion within a few years. These were the early years of the theory i 
of dam design developed by Sazilly, Delocre, and Rankine, and not much ee 
had materialized from it except the rule of the middle third. During the 
ensuing third of a century many dams between 100 and 200 ft in height, each 
- comprising more than 500 000 cu yd of masonry, were designed on this prin- 
~ ciple and constructed in various parts of the country. Because the dams 
were long the sections were made as small as possible, and it was considered 
good practice to bring the lines of pressure as near as possible to the limits 
of the middle third. This resulted in sections rather more slender than are 
to-day thought conservative. 

As to materials, they were those which the country afforded. The earlier 
works would not have been possible otherwise. Portland cement, of course, 
was known, but it was manufactured in but few localities, was costly, and was 
little understood. Concrete was practically unheard of. The origin of surkhi 


Norn.—The paper by D. W. Cole, M. Am. Soe. C. B., was published in February, 1935, 
Proceedings. Discussion on the paper has apeared in Proceedings, as follows: August, 
1935, by Messrs. Oren Reed, F. F, Fergusson, and Joseph Wright. 

: 6 Jackson Heights, N. Y. 
6a Received by the Secretary July 17, 1935. 


t deh 12 sates ‘aid 79 tons at 18 months; RE ee 1060 
30 Ib per sq in. Such material cannot be summarily dismissed 
sideration. - 
ee et ~The Indian Sieh are essentially agriculturists. They are not pohianieatom 

ay minded. Elaborate and complicated machinery does not enter into their 
scheme of things, although their skill and patience in many handicrafts are My 
revelations to the westerner. Time is not important to them. They do not 
think in terms of the individual’s span of life. 
_ In these circumstances it is not surprising that the administrators of the 
Indian public works designed and built as they did—with a rudimentary. 
theory, native materials, a minimum of mechanical equipment, and a maxi- 
mum use of manual labor of the kinds to which the people were accustomed. 
After half a century without a major failure, notwithstanding the construc- 
tion of a great number of large and important works, they were justified in 
considering their practice sound, and it is not strange that they were slow 
to adopt the later developments in materials and methods. 

With this background the engineer who designed the Walwhan and 
Shirawta Dams, and who had spent most of his long professional life in India, 
naturally preferred to follow a tried and proven road rather than to pursue 
what might be only a will-o-the-wisp. : 
at The Shirawta Dam contains about 650 000 cu yd of masonry. To construct 

-—~——s such a volume, almost entirely of one-man stone, by hand labor within a few 
4 working seasons required an army of men and women. These literally 
+s 
i 


swarmed all over the work. In fact, seen from an elevation and from a little 

distance the activity resembles nothing so much as a gigantic swarm of bees. 
be The labor turnover is enormous, for these people will work on outside employ- 
a ment only during such times as their personal agricultural operations do not 
require attention. 

Herein lies the weakness of the system. Really good masonry can be 
and has been executed by this somewhat primitive procedure, but with hun- 
dreds of masons distributed over a large area and served by thousands of 
coolies who bring the stone and mortar, adequate supervision and inspection 
are impossible. Here and there a careless or indifferent mason, an occasional ~ 
batch of improperly ground mortar, a few stones which are dirty or cracked— 
the result is a weak spot, a discontinuity, to and through which water under 
pressure may find its way. Time and the never-ending erosion by running 
water do the rest. 

Thus the slowly but steadily increasing leakage through these dams was 
regarded as symptomatic of a condition which should not be allowed to per- 
sist. The loss of water was unimportant, at least in any quantity which 
might be reached for some years. It was not believed that the dams were 
unsafe, but it seemed certain that they would become so if the deterioration 
should be permitted to continue. If that time should come, every one would 
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deprive the ‘City of Bombay and surrounding territory of tramway servi e 
light, fans, and refrigeration. It would shut down about 80 cotton mi Bais 
; and deprive the railroads entering Bombay of power. In short, it would 
_ mean ruin and pestilence to a million and a half people. The Indian Govern: Si 
_ ment does not take a fiendish delight in destroying private capital, but even 
i 4 ‘it were indifferent to such destruction it could not face this possibility. On | 
the other hand, failures of the dams would mean frightful havoc to agricul- 
E tural regions below, and almost certainly much loss of life. Ee: are 
Ca Clearly, the damage already suffered must be repaired and the structures 
» must be rendered permanently immune to future deterioration, without inter- 
4 fering with the use of the reservoirs. The plan adopted after thorough investi- 
_ gation and careful study was endorsed by some of the most experienced engi- 
neers in India and received Government approval. The engineers for the con- 
_ tractor, Francois Cementation Company, Ltd., believed, in the light of their 
__ experience, that it promised success. They had not previously undertaken a 
: precisely similar job, but no major engineering problem is an exact duplicate. 
4 
5 
“ 


_ of any other. It was believed that cavities existing in the masonry, whether 

_ resulting from carelessness during construction or from subsequent washing 

out of mortar, could be filled with good cement grout, and that the old surkhi 

_ mortar in the neighborhood of these cavities might to some extent: be impreg- 

-- nated with cement. It is the execution of this plan which has been so well 
deseribed by the author. 

It has been objected by some persons that the degree of betterment is 
uncertain since there is no visual evidence that all cavities are filled wholly 
— or in part. No visual evidence, it is true, but the volume of cement retained. 
in the structures (in the case of Shirawta 0.4% of the total volume of the 
masonry) is very convincing, especially as the pressures employed were suffi- 
cient to force the grout through minute openings. It is not claimed, of course, - 
that the dams are of as high quality as though they had been originally laid Ag 
up in cement mortar under close inspection, but they are certainly as good as _ 
(probably even better than) any of those built with surkhi mortar in the 
manner prevailing in India, and many of which have served for more than 
fifty years without manifesting any sign of distress. 

It has been said that this work was almost unprecedented. The grouting 
of rock formations underlying dams prior to construction is a well-established 
practice, as is also subsequent grouting of parts of the structure through 
holes left for that purpose. 

Grouting at the Camarasa Dam in Spain, completed in 1931, was not a 
parallel case to the Indian dams. At Camarasa the dam was absolutely tight 
and sound; leakage through the surrounding and underlying rocks was the 
fault to be remedied. This object was accomplished with marked success. 
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middle of the Nineteenth Century. As a result of political Hat hte, the 
undation work was not properly done, cracks and settlements supra 
y almost immediately, it was impossible to maintain the pool level which had ~ 
been intended—the dam was a failure. When Sir Colin Scott-Moncrieff — 

; became head of the Public Works Department in Egypt he took immediate 
steps to remedy the difficulties. The final stage of the repair job was to. 
| aril five vertical holes along the axis of each pier and to inject cement 
- grout. No pressure was used except that due to static head; however, cement 
- injected at one pier frequently appeared in the holes in adjoining piers. One 
pier took as much as 78 tons of cement; a total of 1000 tons was used in 
- 182 piers. This work was completed in 1898 or 1899; the dam has given no 
trouble since. 

One work which in some respects niet the Indian dams is the 
“barrage de ’?Oued Fergoug” in Algeria. This dam, 34 m high, also was 
designed in the early days of the middle-third theory, construction com- 
pleted about the beginning of 1872. Like the Indian dams, it was built of 
on one-man rubble by native workmen using lime mortar, perhaps feebly 
_ hydraulic. From the beginning it leaked profusely at many places, and lime 
deposits on the down-stream face indicated serious destruction of, the mortar. 
In December, 1881, a length of 125 m of the highest part of the structure was 
carried away, most of it clear down to the foundation. Between 1881 and 
1885 the gap was filled by a new structure of slightly greater section than 
the old, but built in the same manner and with the same kinds of materials 
except that some hydraulic lime and some Portland cement were used in places 
regarded as especially critical. The dam continued,to leak but in 1900 it 
withstood an unprecedented flood with no indication of weakness. In — 
mi November, 1927, there was another major failure, carrying away a part of 
wt the old and a part of the later work, and resulting in heavy loss of life and 

property. For some time prior to the failure the authorities had been dis- 
turbed over the persistent leakage and the continual leaching of lime from 
the mortar, and had approved a project for drilling and the injection of 
cement. This work was in actual process near the left bank at the time 
, of the failure, but had not progressed far enough to throw any light on the 
~ condition of the masonry. Whether this method of reinforcement and stabili- — 
zation would have been effective if undertaken a few years earlier can only 
be conjectured. 
The outstanding peculiarity of hydraulic construction in India is the 
general aversion to the use of Portland cement. This was understandable 
in the early days when cement would have had to be imported from Europe 
and would have been very costly. To-day, however, excellent Portland cement, 
meeting all standard specifications, is manufactured at several places in India 
and can be purchased at the factory at prices of about 35 rupees per ton— 
equivalent to about 9 shillings per bbl—so that this reason is no longer 
operative. 
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: “Indian ‘engineers never use “Portland cement in the constraction 
masonry dams, whether of stone or of concrete, except perhaps in wet foun-— 
dations; not merely on account of its high cost—though good cement is now 
_ being made in the country more cheaply—but also because it is considered 
] inferior to good hydraulic-lime mortar for dams. Cement work is too rigid, | 
_ whilst the hydraulic-lime mortar, being somewhat elastic, withstands the 
variations of temperature without cracking. There is therefore not the some- — 
what alarming necessity of providing expansion joints. The range of tempera- 
/ ture is probably not greater in India than in England, but it may occur _ 
- more suddenly. % es: 


_ These are arguments Foie to justify an existing practice rather than — : 
reasons on which that practice is based. The reader may judge for himself 
of their validity. — 
J With increasing education of Indian masons in the properties of Portland _ 
cement, and with the gradual change in personnel of the Public Works © 
- Department with the retirement of the seniors and the advent of younger | 
men, it is probable that Indian practice in dam construction will eventually 
come to conform quite closely with that in other parts of the world. Tradi- — 
a tion is strong in the East, and the inertia of 350 000 000 people is not easily 
overcome. 
5 One of the most valuable features of this paper is that it emphasizes the 
dominant importance of execution, as compared with design, in the construc- 
‘4 tion of masonry dams. Volumes of mathematical formulas and theoretical 
- discussions in many languages have come and gone, each purporting to be 
the last word as to stresses in dams; yet sections have changed little in 
seventy-five years. On the other hand, with few exceptions, all important dam 
failures have been traced to defective materials or faulty execution. A dam 
is made safe, or otherwise, on the job—not in the office. : 
y The author is entitled to great credit for his masterly handling OL 2 4 
difficult job. He is further to be commended for the exhaustive and incisive 3 
~ manner in which he has presented the subject in print. There is nothing 
superfluous; nothing has been omitted. Aspiring authors of engineering 
papers should be referred to this as a model. 


7™“The Tata Hydro-Hlectric Power-Supply_ Works, Bombay,” by the late Robert Bat- 
son Joyner, Minutes of Proceedings, Inst. C. E., Vol. CCVII, p. 55. 
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_ By I. K. SILVERMAN, JUN. AM. Soc. C. E. - 


I. K. Smverman," Jun. Am. Soc. ©. E. (by letter).“*—That thin shells, | 
_ properly supported, serve as very economical structures, has been known for 
some time. Mr. Schorer, for the first time, has introduced into American ~ 
engineering literature the theory of such shells which has been developed — 

largely by German and Swiss engineers. 


Water Surface 
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Fie. 9. Fie. 10. { 

Even with the simplifications introduced in this paper, the application of ; 
the theory requires considerable computation as, for example, when the struc- ; 
ture analyzed is not a complete surface of revolution and the loading is not i 
symmetrical. Such a condition occurs in the analysis of the structure such _ 
as the sector gate shown in Fig. 9. : 
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: sses may occur in the longitudinal supporting memt 
¢ y in the member at ¢ = 0. ‘These large stresses are due to the f 
the inevitable arch action developed in the curved skin-plate causes 
- member to bend about Axis 1-1, Fig. 10. Furthermore, stresses of the typ ieee 
'S, subject the longitudinal supporting members to loads which are axially Me 
eccentric. In combination with the direct water load, these stresses may  __ 
; reach large values, especially in gates of long spans. y 
In outline, the method of attacking such a structure as that presented in 
f Fig. 9 is perfectly straightforward. The author has shown that the loads ona 
structure are carried by a primary system of stresses, the membrane system, 
- and by a self-equilibrating system (secondary), the function of which is to 
preserve geometrical continuity. The stresses of the secondary system are 
; applied along a line common to the shell part of the structure and its 
_ supporting members. _ 
The stresses and deflections of the membrane system are obtained as if the 
_ shell were a complete surface of revolution by means of Equations (101) to 
eS (110). The stresses at ¢ = 0 and ¢ = ¢zx (it is assumed now that inter- 
. mediate longitudinal members have been removed) are applied to the sup- 
porting members as reactions, and the deflections and stresses in these 
members are obtained in the usual manner. In general, there will be 
differences in these stresses and deflections, and it is the purpose of the 
secondary system to annul such geometrical discrepancies. It is the object 
of the author to derive the stresses and deflections in the shell due to this 
secondary system. When these stresses and deflections have been obtained 
““s in the shell for the various conditions of loading (see Fig. 4), the correspond- 
ing stresses and deflections are obtained for the supporting members. Br 
In general, there will be eight elastic equations for the indeterminate 4 
-- quantities in the secondary system. In a symmetrical structure under sym- < 
metrical loading, such as that treated by the author, these equations reduce 
to four. 
Finsterwalder“ has reduced the number of unknowns to six by assuming 
- that @ — 0 at the connection of shell to supporting member. This condition 
is quite possible because of the relative rigidities of shell and boundary 
member. Corresponding to Equations (91) to (94) for symmetrical loading 
there would result three sets of four simultaneous equations each for the 
determination of the constants, A, B, CO, and D, namely: in Equation (92), 
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i = :-N; ='0; and Nz = 0; in Equation (93), T, = 0; 1 ee and 
c 
c 


os 4 ; 
T. = 0; in Equation (94), os == (0s Coes O: and = = -— ; (and in Equa- 


3 = 


ox dx Ou G 
tion (98), 0 = 0; 6 = 0; and 6 = 0. 
In all cases, M, (Equation (91)), would be different from zero. With 


the constants known the method of procedure follows that of the author 


3 i ee 

Edad ie der zylindrischen Schalengewdlbe, System Zeiss-Dywidag,” von Dr. 
Ing. Dine seo walder, International ‘Assoc. for Bridge and Structural Eng., Zurich, 
1932. 
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in which the subscripts, s and b, refer to shell and beam, respectively, 
hae ae yee stress. When the loading is Speeds or ' the eee mem 


s ‘ ara aan of shell due to membrane system at Ss = 0; wv, = aaa ‘ara 


shell due to membrane system at ¢ = ¢x; w. = deflection of beam due to 


f membrane system at ¢ = 0; w’, = deflection of beam due to membrane sys- 
tem at 6 = dx; ws = deflection of shell due to secondary system at ¢ = 0; 
w’, = deflection of shell due to secondary system at 6 = $x; ws = deflection 
of beam due to secondary system at ¢ = 0; and, w’, = deflection of beam due 
to secondary system at ¢ = ¢x. These deflections are broken up so that a 
symmetrical and contra-symmetrical system of eager stresses may be 
applied. 

For the symmetrical system: = 


Ary = Oat wl) __ (ot t's) on he's) (ou a ') _ 47) 
2 2 2 2 ; 

and, for the contra-symmetrical system: 
paps Ma W's). Ws ws) 4 me) eo! yea 
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Corresponding equations are obtained for v and g. Unsymmetrical loadings 


tend to cause large secondary stresses, and it is to be expected, in the case — 


of the sector gate, that these high moments and stresses occur at ¢ = 0. 
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_ PHOTO-ELASTIC DETERMINATION OF 
2 ‘SHRINKAGE STRESSES 


Discussion 


oe By MEssrs. THOMAS H. EVANS, AND I. K. SILVERMAN 
“4 


« THomas H. Evans,° Juy. Am. Soc. C. E. (by letter).°—The problem of 
Pe determining, with some exactness, the stresses in a concrete dam due to shrink- 
age seems well adapted to the photo-elastic method. Although the writer disa- 
" grees with the author’s contention that the stress problem is only two-dimen- 
sional, it will not greatly affect the “philosophical picture” to consider it as 
such. The principal inaccuracy in the determination is not in the method 
itself but, as Mr. Smits points out, in the physical inconsistency between the 
prototype and the necessarily idealized model. However, the method does 
provide a rational first approximation to an extremely complex stress dis- 
tribution. It should be a distinct aid in supplementing the engineering 
experience and judgment of the designer. 5 

As indicated by the author, the “membrane analogy” as a supplement to 
_ photo-elastic analysis appears to have decided advantages over the method of © ; 
graphical integration for the determination of the principal stress ~ 4 
sum. There was also another excellent paper by Supper and McGivern’ that ie 
exhibited in great detail the use of thd membrane analogy for the determina- 
tion of the quantity, p + 4. The accuracy of this method decreases, how- 
~ ever, as the slope of the membrane increases. This occurs in regions of high 
stress concentration, such as the heel and toe of the dam illustrated in the 
_ paper. Since a designer might have more interest in such regions than in 
others, it would be desirable to avoid the effects of this limitation of the 
membrane analogy. 
Another supplementary procedure which seems to be adaptable to all cases 
is that proposed® by Professor V. Tesar. This is a purely optical method . 
requiring even less additional apparatus than the membrane analogy, and 
its accuracy is independent of the concentration of stress in the model. Its 


Norp.—The paper by Howard C. Smits, Esq., was published in May, 1935, Proceed- 
ings. This discussion is printed in Proceedings in order to bring the views expressed 
before all members for further discussion. 

® Asst. Prof. of Civ. Eng., Dept. of Eng., Univ. of Virginia, Charlottesville, Va. 

6a Received by the Secretary June 25, 1935. 

7 Journal, Franklin Inst., April, 1934, 

8 Revue @Optique, March, 1932. 
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a a and, in Figs. 10 and 11, presents curves for the resulting stress distribu 


i deties Patera in the bakelite siodl be a ‘aninealae neat resting on a cont 
~ ous elastic foundation and subjected to shrinkage? In this paper, Mr. 
has described a conception of this phenomenon. His description of how tl 
hy _ stresses might vary seems logical. Based on this preconceived picture, t 
proceeds to load his model so as to obtain this possible distribution of st 


near the base. He has failed to show, however, that this is the only way this | 
stress pattern could have been produced and, for. this reason, and for others 
cited herein, no great value can be attached to the numeri- 
cal results shown in Fig. 10. a 
The author states that the form of the test specimer 
was chosen so that the stress pattern would be a gauge of fe 
the degree of axial loading. It is evident that in choosing» | 
this type of specimen, he has clouded the stress distribution — 
in the very region that he is examining. The stress dis- 
tribution in any model is dependent to a high degree on i 
geometry. From Fig. 8 the width of the strap is of the 
same order of dimension as that of the base of the dam. — 
For this reason the stress distribution of Fig. 7 may be far — 
different from that in a model which is of the form shown 
in Fig. 12. | 
By means-of Fig. 9 Mr. Smits has shown the accuracy 
obtained in measuring the stresses shown in Fig. 7. It 
would be interesting if, in his closing discussion, he would 
present curves similar to those of Fig. 9 for the resultant stress distribution © 
shown in Fig. 10. Evidently, the total areas under these curves should be zero. 


* With U. S. Bureau of Reclamation, Pare Colo. 
®a Received by the Secretary August 8, 5. 
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